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I. Introduction 


For years man had desired the capability to remotely monitor 
various' phenomenon in his environment. For example, the measure- 
ment of ocean wave and wind conditions are of vital interest to 
many marine industries and government agencies. Furthermore, the 
knowledge of ocean-surface temperature on a global, all-weather, 
and day-night basis is also of importance to the fishery and marine 
transport industries, as well as the oceanographers and marine 
meterorologists, and weather forcasters. In recent years, the 
microwave radiometer has proven itself to be a feasible remote 
sensing device. To monitor the environment on an all-weather basis 
microwave sensing has the immediate advantage of being affected 
less by fog and rain than infrared. In addition, microwave radio- 
meters have been designed [1] that can measure the incident 

O 

radiation to an accuracy of j^O.l K and remain calibrated, 
unattended, for a year or more. In fact, small, lightweight, auto- 
mated radiometer systems have recently been flown by NASA on the 
Nimbus satellite program [2], There is presently a great deal of 
research and development being conducted in the area of microwave 
remote sensing from satellites. 

In order to make precise measurements of the radiometric 
brightness temperature of a target (and thereby infer certain physi- 
cal parameters) one must be able to mathematically model the inter- 
action between the electromagnetic radiation properties of the 
antenna and the incident radiation from the environment. This. 



interaction can be described by Fredholm integral .equations of the 
first kind which are extremely unstable. This instability has been 
studied in considerable detail by investigators in many fields, 

Twomey [3] and Phillips [4] have devised matrix filtering techniques 
to stabilize the solution. Although these matrix methods are not 
without merit, Bracewell and Roberts [5] have demonstrated the value 
of a successive substitution solution. Assuming that the intensity 
of the emitted radiation of the environment can be represented in 
scalar form, they have shown that the antenna is only capable of 
responding to those frequency components of the function representing 
the environment below a cut-off determined by the antenna aperture. 
The high frequency components of the emission function are invisible 
to the antenna. The low frequency components are accepted but their 
relative magnitude is altered according to the system (antenna) fre- 
quency characteristics. Inversion through the method of successive 
restorations leads to the principal solution [5], in which frequency 
components accepted by the antenna have been restored to their 
original values, but the rejected components are not represented in 
the solution. The work done by Bracewell and Roberts was, however, 
more applicable to astronomical observations than to general micro- 
wave radiometric measurements. They assumed that the antenna was 
very efficient and that the sidelobes and backlobes could be 
neglected, which is not always the case. They also used a scalar 
representation of the interaction between the antenna radiation 
characteristics and the emission by the target. 



The interaction between the emitted radiation from a water 
surface and the radiation characteristics of an antenna is a 
vector relationship. As the sidelobe and backlobe levels of the 
system (antenna) weighting function become more intense, and the 
major lobe beamwidth more wide, the vector model interaction 
becomes more important. Classen and Fung [6] have vectorially 
modeled the viewing of the ocean using matrix techniques. Their 
representation, however, assumes that the observed environment is 
circularly symmetric and infinite in extent. It should also be 
pointed out that computer time using the matrix modeling would be 
quite extensive as compared to some other types of numerical 
techniques.. 

To study the radiometric signature of a controlled water sur- 
face, a wave tank system has been constructed at NASA Langley 
Research Center, Hampton, Virginia. For the wave tank geometry, 
the environment is of finite extent and is no longer circularly 
symmetric. The response of the radiometer for this system was first 
modeled by Fisher [7], using a two-dimensional scalar approximation. 
This approximation works well for high efficiency antennas. The 
direct inversion used by Fisher [7] was, however, sensitive to errors. 
Holmes [8], by applying the iteration techniques of Bracewell and 
Roberts [5] to this problem, was able to restore, with acceptable 
accuracy, the brightness temperature (scalar emission function) of 
the water from measurements that contained error. Both Holmes [8] 
and Fisher [7] used the Fast Fourier Transform techniques, with an 



algorithm reported by Fisher [9], to perform their computations. 

The first three-dimensional modeling of the NASA wave tank 
was done by Beck [10]. He was able to formulate and calculate the 
antenna response of the system given the emission characteristics 
of the surroundings. Beck's formulation requires numerical inte- 
gration for direct computation of the antenna' temperature and is not 
convenient for inversion processes nor can it be modified con- 
veniently for efficient and economic restoration computations. 

The classical design consideration of a radiometer antenna is 
the compromise between resolution and system design constraints 
(size, frequency, etc.). Only after the inversion process has been 
studied is the true resolution of the antenna known and can the 
design for the particular application be made. The study of this 
inversion for- the wave tank geometry is the subject of this 
dissertation. A three-dimensional inversion scheme is described 
which takes into account the interaction between the radiation 
characteristics of the antenna and emitted radiation from the wave 
tank (vector representation), and computations are performed using 
the efficient and economical Fast Fourier Transform algorithm. The 
inherent instabilities of the inversion are overcome by the 
adoption of the filtering properties of the restoration method. 
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II. Theory 

A. Brightness Temperature 

All matter above absolute zero temperature emits electro- 
magnetic radiation due to the thermal motion of its atoms/mole - 
cules. The brightness temperature of a given substance is a 
standard measure of the intensity of this radiation. By definition, 
the brightness temperature of a perfect black body radiator is equal 
to its molecular temperature. For the perfect black body radiator, 
none of the electromagnetic radiation generated from within the body 
is reflected back at the interface between the radiating surface 
and the surrounding transmission media. For all passive physical 
objects, however, the transmission coefficients for the radiating 
surface are less than unity. Consequently, the brightness 
temperature will be less than the molecular temperature. The trans- 
mission coefficient is often called the emissivity, and the 
brightness and molecular temperatures are related by 

■^b = 0;) 

where is the brightness temperature, T„i the molecular temperature, 
and £ the emissivity or transmission coefficient. 

For a flat semi-infinite radiating surface the emissivity can 
be found from the complex dielectric properties of the radiator. 

The emissivity is also a function of both the incidence angle at 
which the interface is viewed and the polarization of the emitted 
wave. Stogryn [11] and Holmes [8] have shown how the emissivity is 
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related to the complex permittivity of the radiator. For the 
perpendicular (horizontal) polarization the E-field is perpendi- 
cular to the plane of incidence and the emissivity is given as [8] 



II 

6 = incidence angle 
p = \TF cosy 
q =\/r" siny 



e = Re[e^] 

H 


e^Hcomplex dielectric constant of 
the radiator 

For the parallel (vertical) polarization the E-field is parallel to 
the plane of incidence and the emissivity is expressed as 
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Through (2) and (4), the radiation Is being described as two 
orthogonal,- linearly polarized waves. The radiation from the sur- 
face of the water can be described in this manner if the dielectric 
properties of the water are known. Stogryn [11] has concluded that 
the dielectric constant for sea water may be adequately repre- 
?;pntprf hv thp following equation of the Debye form 


^ . jo- ■ 

. ... ‘ ^ 

1-j 2it T f 2-^ e* f 



where and are, respectively, the static and high frequency 

dielectric constants of the solvent, x the relaxation time, e* the 

_12 

permittivity of free space (= 8.854x10" farads/m), ct the ionic 
conductivity of the dissolved salt in mhos/m, and f the electro- 
magnetic frequency. 

In order to evaluate (5), the variations of x, and a as 
functions of salinity, frequency, and temperature need to be known. 
■By using resonant cavity techniques, Stogryn [11] reported, through 
numerous measurements, empirical equations to evaluate the variables. 

The high frequency dielectric constant is considered to be 
a constant (=0.48), The low frequency dielectric constant and the 
relaxation time x are expressed as 


£o(T,N) = CqCT.O) a(N) (6) 

2itx(T, 1^) = 2-irx(T,0) b(T,N) (7) 

0 

where T is the water temperature in C and N is the normality of 



8 


the solution. The series expansions used to evaluate (6) and (7) 
for 0 < T 5 40 C and 0 ^ N < 3 are 

a(N) = 1.0 - 0.2551 N + 5.151 x 10"^N^ (8) 

- 6.889 X lO'" n' 

b{N,T) = 1.463 X 10“^ N T + 1.0 - 0.04896 N (9) 

- 0.02967 + 5.644 x 10“^ 

£q{T,0) = 87.74 - 0.40008 T (10) 

+ 9.398 X 10“'' + 1.410 x 10“^ t' 

2ttt(T,0) = 1.1109 X 10-^“ - 3.824 X 10"^\ (11) 

_14 2 16 3 

+ 6.938 X 10 T - 5.096 x 10" T 

Given the salinity in parts per thousand, the normality can fae found 
as 

N = S (1.707 X 10"^ + 1.205 x 10"^ S + 4.058 x lO"^ S^) (12) 

The series is valid for 0 < S < 260. The expression reported for 
the conductivity a of sea water is 

o(T,S) = a(25,S) e"^^ (13) 


where A = 25 - T and 


C = 2.033 X 10" 
-S (1.849.x 10"^ 


2 


4 


+ 1.266 X 10“ A + 2.464 x 10" 
- 2,551 X 10"’ A + 2.551 x 10" 



C14) 
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a(25,S) = 5(0.182521 - 1,46192 x 10"^ S 

5 z 7 3 V,l3J 

+ 2.09324 X 10“ S - 1.28205 x 10" S ) 
in the range 0 s S < 40. 

Using (6) - (15) we can obtain e and e from (5). With e' 

tl |l 

and e, £|^(0 ) and £^(0 ) are found by the use of (2), (3), and (4). 
As seen by (1), the horizontal and vertical brightness tempera*- 
tures of the polarized radiation emitted by the water are 


° 06a) 

06b) 


Equations (16a) and (16b) yield the intensities of two linearly 
polarized, orthogonal waves that are needed to describe the radi- 
ation emitted from the water. 

Brightness temperatures of the earth and sky were also part of 
this investigation. These brightness temperatures have been found 
(experimentally) to be nearly randomly polarized and therefore 
related to molecular temperature by (1). For the sky, Peake [12] 
expressed the brightness temperature as 

“T-sec0_ 

TbsC^s) = TeffO - ^ 07) 

where 


08 ) 
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To = (19) 

The angle 6^ is the angle measured from zenith. For the lack of a 
more accurate brightness temperature model, the earth emissions are 
usually assumed to be constant and unpolarized. If a more accurate 
polarized brightness temperature model were known, it could be 
utilized in the analysis and computations. 

In addition to its own generated radiation, the water surface 
reflects the incident sky radiation and directs it toward the 
receiving antenna. To account for this reflection, (16a) and (16b) 
can -be modified as 

^bwh(»") = + ('-e^)T|, 3 (e,=e'') (20a). 

’'bwv<8") = + (l-ev)Tt,s(es=8“) (20b) 

In Figure 1, we have plotted and Tj^^^ as functions of inci- 
dence angle for T^^, = 284°K, S = 0 Q/oo, and f = 10.69 GHz. The 
shape of the plots are basically the same for any temperature, 
salinity, and frequency. The peak in the Tj^^^ curve occurs when 
the water is viewed at the Brewster angle (e^ =1). At this angle, 
^bwv equal to the molecular temperature of the water. 

B. Antenna Temperature 

The antenna temperature measured by a radiometer is the 
brightness temperature of the observed environment weighted by the 
power pattern of the antenna. We shall define Tj^(e.,0) as the 




Fig. 1. V/ater brightness temperatures. 
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unpolarized brightness temperature of the environment, T the 

a 

measured antenna' temperature,- and G(e,0) the antenna power pattern 
which has been normalized so that the integral of G(e,0) over the 
entire solid angle is equal to unity. The variables are related by 
the following relationship 


2ir t'R 


T 


a 


T|j(e,0) G(e,0) sine de d0 


( 21 ) 


J J 

0 0 

If G(e,0) were a delta function 6(e-6„,0-0_) ,T, would then equal 
^b^®o’^o^’ P’^sctical antennas, however, do not have such convenient 
radiation characteristics, and is generally not equal to the 
T|^ at boresight. 

If a significant fraction of the emitted radiation from the 
observed environment (brightness temperature) is polarized, such as 
that emitted by the water surface, (21) is no longer a valid expression 
to be used to calculate the antenna temperature. To explain the 
coupling between the radiation properties of the antenna and emitted 
polarized radiation from the environment as well as the concept of 
partial and total antenna temperatures, let us assume that the radio- 
meter system is over the ocean in clear atmospheric surroundings. 

Since the observed environment is the water and sky, the total antenna 
temperature T. is equal to the contributions from the water T„, and 
sky T^g. Assuming that the radiation from the sky is unpolarized, T^^- 
is expressed as 
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J J G( 0 , 0 ) sine de d0 

over . 
skv 


( 22 ) 


Since the radiation emitted from the water is polarized, to 
find its antenna temperature contribution, the weight of the gain 
function 6(9,0) needs to be found at each integration point in 
directions perpendicular and parallel to the plane of incidence. To 
do this, we form the unit vectors h(0,0) and v(6,0) within the water 
integration limits. The vector h(6,0) is perpendicular to the plane 
of incidence formed at the integration point on the water surface, 
and v(0,0) is orthogonal to h(0,0) and r(0,0), where r(9,0) is 
the radial unit vector. For a given antenna, the normalized 
electric field intensities in the 6(0,0) and 0(0,0) directions. 

Eg (0,0) and E0(9,0), can also be found. In turn, the power 
intensities, at each integration point, for the horizontal and 
vertical polarizations, 6^ and g'^, are then formulated as 

6^0,0) = [h . '0 Eg + h . 0 Ejj]" (23) 

G'^(e,0) = [v • 9 Eg + V . 0 E^]^ (24) 

The antenna temperature contribution from the water can then be 
expressed as 
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Tb^f^(0i0) G^{0,0) s1n0 de d0 

over 

water 



r 


j 


over 

water 


g'^( 9,0) sine d0 d0 


The total temperature measured by the radiometer is 


(25) 


T = T 
a as 


+ T. 


aw 


(26) 


Equations (22) and (25) define the relationships between the 
power pattern of the antenna, the brightness temperature functions of 
the observed environment, and the measured antenna temperature for 
both unpolarized and polarized emissions. 


C. Wave Tank Geometry and Theory 

In order to obtain the microwave emission signature of a 
water surface in a controlled environment, a wave tank system has 
been constructed at NASA Langley Research Center, Hampton, Virginia. 
The model, as illustrated in Figure 2, consists of a fourteen foot 
square tank with the antenna and radiometer placed at the end of a 
boom over the tank. The antenna and radiometer can move along a 
circular arc above the tank and can be scanned, at each position, 

O 

through a complete 360 in a plane which bisects the wave tank. The 
angle g is the scanning angle and a describes the position of the 


boom. 



Earth 


Fig. 2. Radiometer and Finite Wa 
NASA Langley Research Cei 
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tenna 



e Tank Configuration at 
ter, Hampton, Virginia.' 
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For the wave tank measurements, the total antenna 
temperature is composed of three partial antenna tonperatures; 
namely that of the water, earth, and sky. The partial antenna 
temperatures are of the same form as (22) and (25) and are given by 




T^^^I^(9J) G*^(0,0) sine .d9 d0 


over 

water 


■ ■ 

■ s'! (9.0) sine de d0 


over 

water 


(27) 




Tt,g(9,0) G(0,0) sine de d0 


over 

earth 


(28) 



J J 
over 

sky 


T[^g(0,0) G(e,0) sine de d0 


(29) 


h w 

where G and G are defined by (23) and (24), respectively. The 

total antenna temperature T, is then equal to T + T + T . 

a ^ aw ae as 


1 . Z-Axis Normal to Radiometer Antenna Aperture 

Patterns from directional antennas used in radiometry, 
such as horns, are nearly circular symmetric about the boresight. It 
would therefore be convenient to express the gain functions in a 



coordinate system v/hich uses the z-ax1s as the boresight. The 
problem was .origitially formulated 1n this manner by Beck [10], The 
coordinate systems used are illustrated in Figures 3 and 4. The 
origin of the x,y,z coordinate system is the phase center of the 

I I I 

antenna and the origin of the x, y, z coordinate system is the 
center of the wave tank. Although the radiometer antenna may be of 
any type, let us assume one with an aperture E-field polarized in 

* •^ I t 

the X direction (aperture E-field parallel to the x-z or x - z 
plane of Figure 3). This type of antenna has a strongly polarized 
pattern. As the antenna is scanned parallel to the x-z plane, as 
shown in Figure 3, the antenna will principally see the vertically 
polarized emissions from the water. When scanned parallel to the 
y-z plane, as shown in Figure 4, the horizontal polarization will 
predominate. In the system implementation, if the boom is allowed 
to. move only along one plane, the system polarization can Be changed 
by simply rotating the antenna aperture 90° about the z-axis. The 
system at NASA has, however, the capacity to move along either plane. 
To describe both rotations, shown in Figures 3 and 4, with one set 
of functions, a fictitious constant 0^^ is introduced which is set 
equal to zero when the rotation is as shown in Figure 3 (vertical 
polarization scan) and equal to -rr/2 for the scanning displayed in 
Figure 4 (horizontal polarization scan). 

The two variable brightness temperature profiles in (27), 
^bwh^®’*^^ and ’^’ 5 y,v( 0 > 0 ) > can be expressed as a function of a single 
incidence angle variable 0. To find 0, one begins with the relation- 
ship between the primed and unprimed rectangular unit vectors of 
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A 


y 



Fig. 4. Z-axis Coordinate System Orientation for 
Horizontal Polarization. 



Figures 3 and 4 which are written as 

X = X (cosa cos0j^ + sinP^) + z sinct 
y = y (cosa sinP^ + cos0^) + z sina 

A A j A I ^ j 

z = - X sina cos0q - y sina sin0j^ + z cosa 
The angle 9" can be expressed as 

ccse" = R . £' / |R| |z'| 

where 

^ A /s A 

R = R r = R(x sine cos0 + y sin9 sin0 + z cose) 

Substituting (30a), (30b), and (30c) into (32), one obtains 
R = R(R^ x' + Ry y' + z') 

R = R {x' [sine COS0 (cosa cos0^ + sin0^) - sina cos0^ cose] 
+ y [sine sin0 (cosa sin0Q + cos0^) - sina sin0Q cose] 

+ z [sine COS0 sina cos0q + sine s1n0 sina sin0Q 


(30a) 

(30fa) 

(30c) 

.(31) 

(32) 


+ cose cosa]} 


(33) 
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Since the magnitude of the unit vector z' is unity and R • z' is 
known from (33), cose" in (31) can now be expressed as 

COS0 = sine COS0 sina cos0p + sine sin0 sina sin0^ 

+ cose cosa (34) 

This allows the evaluation of the incidence angle e" as a function 
of e, 0, a, and polarization (0 q). 

In order to evaluate G^(e,0) and 6^(e,0) for the z-axis 

0^ A 

jeometry, h(e,0,a,0p) and v(ej0,a,0^) must be found. The vectors 

V A ’ 

I and h are defined by the vector relationships 


, h • r = 0 C35) 

A A A A . 

h • z = h • z =0 (36) 

AAA 

V = h X r (37) 

and can be written as 

h = + W^z (38a) 

V = V^y’ + V^z' (38b) 


Equation (36) implies that = 0 and (35) can then be expanded to 
yield 

H, R, + H,, Ry = 0 (39) 

A 

Since h is a unit vector. 
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2 2 

H + H =1 
X y 

(40) 

Solving (39) and (40) simultaneously, one finds that 



(41a) 


(41fa) 

Expanding (37)' yields 


; = Hy R^J' [H, Ry - R^] 2' 

(42) 

Therefore , 



(43a) 

\ = -«x«x 

(43b), 

''z = »x«y - «y«x 

(43c) 


/\ /\ 

The vectors v and h have now been broken into their 

/N 

primed rectangular coordinate components. The unit vectors 0 and 
6 can also be expressed this way to allow the dot products to be 
taken. One must begin with the vectors in the un'primed coordinates 

A /s ^ 

0 = X cos9 COS0 + y cos6 sin0 - z sin0 C^4a) 

/V A. /\ 

0 = -X sin0 + y cos0 


C44b) 



Using (30a), and (30b), and (30c), we can write (44a) and (44b) as 


< X 

X 

1— 

II 

< CD 

1 

+ V 

(45a) 

0 = P,x 

+ PyV' + 

(45b) 


where 


= cose COS0 (cosa cos0^ + sin0^)‘+ sina sin6 cos0^ (46a) 

Ty = cose sin0 (cosa sin0^ + cos0^^) + sina sine sin0^ (46b) 

= sina cose (cos0 cos0^ + sin0 sin0^) - sine cosa (46c) 

= sin0 (cosa cos0^ + sin0^) (46d) 

Py = COS0 (cosa sin0Q + cos0^) ‘ (46e) 

= -sina sin0Q cos0 + sina cos0^ sin0 (46f) 


The dot products in (23) and (24) can now be evaluated using 
(45a)-(46f) as 


e ■ h = 

A j\ 

0 ' h = P^ 
e • V = T 

X 

/N /N 

0 • V = P^ 




'x " P, 


V H- T V + T 
X .V y 'z 


Vx Py V, P, 


(47a) 

(47b) 

C47c) 

(47d) 



To evaluate tfie variables and using (27), 

(28), and (29), respectively, the only parameters still not known 
are the limits of integration and the angle 6^('0,0,a) (the angle 
measured from zenith) needed to evaluate (8^). 

Since the z-axis always passes through the center of the 
wave tank, then for any value of -0, as 0 is, varied from 0 to ir, we 
will always be integrating first over the water, the earth, and ther 

the sky. What is needed then are the values of 0 as a function of 

1 

0 at which the water-earth boundary and the horizon (earth-sky 
boundary) occur. We will define as the water-earth 

boundary and 9g^(0,a,0p) as the horizon. Due to the symmetry of the 
problem, the integration limits of 0 can be made 0 to ir for the 

vertical polarization case and ir/2 to 3ir/2 for the horizontal polari 
zation. 

We will first outline the procedure in determining 0^^. 
Referring to Figure 3, we can write that 


where ' 


cos6 


we 


G ■ Rq 

!g| IR^I 


(48) 


^o "" '^ 0 ^”^ " y sina sin0 + 


cosa) 


(49a) 



(49b) 


The vector G is found by defining the position of each end of the 
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vector referenced to the primed coordinate system. The coordinates 

t 4 

of the point A are x = sina cos0j^, y = Rq sina sin0^, and 
z' = -Rq cosa. The coordinates of the points along the edge of the 
tank are shown in Figure 5, which is a view of the wave tank looking 
straight down, and are given by 


= W, 

1 t 

y = W tan0 

1> 

1 

0 

> 

ir 

(50a) 



4 



4 


= W cot0 , 

y' = W 

■ ^> 
4 

1 

0 

> 

2L 

4 

(50b) 

= -w. 

1 * 
y = -W tan0 


1 

0 

> 

Stt 

(50c) 



4 



4 


= -W cpt0 , 

y' = .-W 

4 

0 ' 

> 

5tt 

4 

(50d) 


For all cases z = 0. This gives four difference expressions for 
G. To eliminate the redundancy of showing the derivations for all 
four cases, we will show the details of finding 0^^ for case 1 
when ^>0' > - ^ and then list for the other values of 0'. 

For case 1 > 0* ^ vector G can be expressed 


as 


G = (W >- Rq sina cos0p) x' + (W tan0* ^ Rq sina sin^g) y' 

+ Rp cosa z' (51) 

By defining as the ratio of W/R^, we can write’ the dot product 
in (48) as 




5. Overhead View of Wave 
Z-axis Geometry. 







G ■ Rq " [(-sina cosP^) (W^ - sina cos0q) 

I 2 2 

+ (- sina sin0Q) (W^ tan(5 - sina sin0Q) + cos a]R^ (52) 
Equation (52) can be simplified into the form 


G • Rq'= [1 - sina {cos0q + sin0^ tan0 )]R^ (53) 

The magnitude of the vector G is given by 

- z t 2 

|G1 = R„[(U„ - si na COS0) + (W^ tan0 - sina sin0^) 

+ cos^a]^ (54) 


which can be reduced to 

|G| = Rq[Wj^ " ^ '^n cos0^ + 1 + tan 0 - 2 tan0 

1 

sina sin0Q]^ (55) 

Substituting (49fa), (53), and (55) into (48) and solving for 6^^ of 
case 1 yields 

1 I 2 

®we ” (cos0Q sin0^ tan0 )]/[W^ - 2W^ sina cos0 

1 


+ 1 + 

•For case 2, 


2 2 1 • 
tan 0 •" 2 tan0 sina sin0g] 

— V. ffl' ^ X 9.,« "is found by 
4 - > 0 > we . 




( 56 ) 
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®we "• (cotp' COS 0 Q + sin 0 ^}]/[W^^ - 2 W^ sina sin 0 ^ 

2 i- 

+ 1 + cot 0 * - 2 cot 0 sina cos 0 q]^} 1[5?) 

Stt ^ 3’t 

When — >0 > — , we have case 3 and for this 
4 4 

_1 j 2 - 

®we " ''^n (cos0^ + sin0Q tan0 )]/[W^ + 2W^-, sina cos0q 

+ 1 + tan^0* + 2 tan0* sina sin0^]^} (58) 

For ” > 0 ' > ~ j e,,Q is given by 
4 - 4 we ^ 

_ 1 I 2 

®we ~ ^ ^n (cot0 COS0Q + sin0^)]/[W^ + 2W^ sina sin0^ 

i_ 

+ 1 + cot 0 + 2 cot0 sina cos0j^]^} (59) 

I 

Given a and 0 one can now find 0 ^^. However, the integration will 
be performed in the unprimed coordinate system, so a relationship 

I 

between 0 and 0 is needed. This can be found from the relation- 
ship, referring to Figures 3 and 4, R = + r'. The vector R has 

already been expressed in the primed coordinate system by (33), Rq 
by (49a), and r' is given by 

r' = r' (x' COS 0 ' '+ y' sin 0 *) (60) 

I *ir » 

since 6 = water sur-^ace. Using (33), (49a), and (60), the 

three vector components of the equation r' = R - Rq yield 
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I I 

R COS0 = R[sin6 cos0 (cosa cos0^ + s1n0^} 

- 3ina cose cos0q] + sina cos0q - (61) 

I t 

R sin0 = R[sin0 sin0 (cosc; s1n0p + cos0q) 

- sina cose sin0^] + sina sin0Q (62) 

0 = R [sine sina (sin0 s1n0j^ + eos0 cos0^) 

+ cose cosa] - R^ cosa (63) 

Substituting (63) into (61) and (62), one can write 

(^) sin0 = sine sin0 cos a sin0^ + sine sin0 cos0^ cosa + sin^a sin 
sine COS0 cos0^ + sin^a sin 0^ sine sin0 (64) 

R ' * 2 

(•j^) COS0 = sine COS0 cos a cos0^ + sine cos0 cosa sin0^ + sine cos0 
.2 . 2 

sin a COS0Q + sine sin0 sin a sin0^ cos0^ (65) 

Since 0^ is either 0 or | , sin0^ = sin^0^, cos0^ = cos'^0q, and 
sin0Q cos0^'= 0, (64) and (65) can be reduced considerably to 

R* - ' 

(^) sin0 = sin0 sine (cos0q cosa + sin0^) [66) 

o' I 

(^) COS0 = CO30 sine (sin0Q cosa + cos0q) 


( 67 ) 
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Dividing (66) by (67), we get the desired relationship between 
0 and 0* to be 

, sin0 (cos0 cosa + sin0^) 

• tan0 = 5 (68) 

COS0 (sin0Q cosa + cos0^ 

With the above relationship between 0*, 0, 0^, and a, (56), (57), 
(58), and (59) can be used to evaluate 0 (0, a, 0 ). 

Finding an expression for 9^^ is considerably easier. 
As the observation point moves farther away from the wave tank, the 
vectors R and r‘ become nearly equal. In the limit, as the obser- 
vation point approaches infinity, R = r'. With this approximation 
for the horizon, equating the z components of R and R yield 


cose = sina cos0^ cos0 sine + sina 


sin0Q sin0 sine + cosa cos0 


(69) 


Solving (69) for 9 , which occurs when e' = ^, results in 

2 


1 cosa 

0 _ r 

es ‘--sina cos0 cos0 - sina sin0 sin0 


rar ] (70) 


The angle 0^, measured from zenith, can also be found 
from (69), since 


0„ = ir - 0 


(71) 
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Using (71) and (69) yields 

6g = cos t-sina cosP^ cos0 sin0 

- sina sin0Q sinj3 sin0 - cosa cos6 (72) 

Now. given T|,3(0j, Tj,^. T^„,,(e“), EgtS.®). 

Epj(0,0) and the tank dimensions, we can now find T ,,, T and I ^ 
as function of a for both polarizations. The scan angle 3 would be 
equal to zero in these calculations. To calculate the T, 's as a 

a 

function of 3 (for a given a) requires a transformation of variables. 
Referring to Figures 6a and 6b, the vertical scanning involves 
coordinate system rotation about the y-axis and for the horizontal 
scanning an x-axis rotation. It can be seen that the antenna gain 
functions will be known in the Xp yp z^ coordinate system or as 
functions of 0-| and 0p To integrate in the unprimed coordinate 
system, 0-| and 0^ need to be expressed as functions of 6 and 0, 
Appendix I contains a derivation of these transformations. For the 
vertical scanning, pictured in Figure 6a, 



( 75 ) 



Fig. 6a. Coordinate System Tra 
Vertical Scanning (6 
Geometry. 








(76) 


0-j = cos" [-sin3 sin0 sine + cosB cos6] 

By using GV(6p0^), G^(0^,0^), and G(0^,0^) in (27), 
(28), and (29), and can be found as functions of a 

and 8 for both polarizations. However, using this geometry, (27), 
(28), and (29) must be evaluated by numerical integration for each 
value. of a and 8. One would also have to solve the transcendental 
equations relating 6-| and 0^ to 0 and 0 at each integration point. 
This would require considerable computer time and can be avoided 
if the scanning of the antenna is described by a rotation about 
the z-axis instead of a rotation about the x- or y- axis as given 
by (73)-(76). 

2 . X-Axis -Normal to Radiometer Antenna Aperture 

An alternate coordinate system that avoids the trans- 
cendental equations describing the scanning is illustrated in 
Figure 7. In this case the x-axis is used as the boresight of the 
antenna for both horizontal and vertical scans. These scans are now 
mathematically described by rotations about the z-axis, and the 
transformations of coordinates during the scan (as shown in Appendix 
1) leave the 0 variable unaffected and change 0 by a constant value. 
The elimination of the transcendental equations is not the only 
advantage of rotating about the z-axis, It will be shown that by 
utilizing this geometry, the integration with respect to 0 and the 
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functional variation with respect to 3 of (27), (28), and (29) can 
be established in a correlation form and evaluated conveniently and 
efficiently by Fourier transform techniques. It was for this 
reason that this system was adopted. 

Since the antenna system is now restricted to rota- 
tions about the z-axis, it is not going to be scanned in two ortho- 
gonal planes' to establish the two different polarizations. Instead, 
the scanning will be restricted in one plane but the antenna orien- 
tation (aperture field) will be changed to accomplish this. To 
receive primarily the vertical polarization, the aperture field is 
assumed to be oriented in the y direction. If the horizontal 
polarization is desired, the Rg and fields are those calculated 
with the aperture field in the z' direction. We shall use the sub- 
script p to denote a function that depends upon polarization. The 
subscript p will represent h for horizontal or v for vertical polari- 
zation. 

To evaluate (27), (28), and (29) with the new geometry, 
we again need to find the dot products that represent the degree of 
alignment between the 0 ' and vectors of the antenna's coordinate 
system and the horizontal and vertical unit vectors. The incidence 

II 

angle 0 and the various limits of integration also need to be 
known. The dot products and incidence angle will be found utili- 
zing the geometry as represented in Figure 8. 

The planes defined by constant values of 0 form on the 
water surface straight lines N which are parallel to the z-axis. 
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Along each one of these lines, 0 is defined as the angle between 
the projection of the radial line into the x-y plane and the 
x-axis. The radial line is R and its projection into the x-y plane 
is M. Therefore, 0 is the angle between H and M. The angle 
between the z-axis and R is 0. Since the line N and the z~axis are 
parallel, the lines R,M,N, and the z-axis all lie in the same plane. 
The z-axis and M intersect at right angles, so the angle between M 
and R is ^ - 6, defined positive in the direction shown, 

A 

The horizontal vector h always lies in the plane that th 
surface of the water defines. Therefore, if one can find as a 

A 

function of G and 0 then h can be found from The plane the water 
surface defines is parallel to the y-z plane and h is given by 

A A A 

h = y sini|> - z cosip (77) 

Line segment length L is found from 

L = H tan0 (78) 

and M by 


M =v = H sec 0 

To find N we use the relation 

^ = tan(|- • 6) = cote 


(79) 


(80) 
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Substituting (79) into (80) we can find N to be 


= H sec 0 cot 0 


( 81 ) 


Using (78) and (81), we find that 


. /N\ ^ -^r cos e T 

4) = tan (j-) - tan 


(82) 


The vector h is now known from (77) and we can get v from the 
relationship 


A A. 

V = r X h 


(83) 


The unit vectors 0, 0, and r are expressed in (44a), (44b) and (32), 
respectively. Using (32) and (77) in (83) yields 


A A 

V = X (“ cos4< sin0 sin© - sin4) cos0) 


A A 

+ y cosijj COS0 s1n0 + z sin^; cos0 sin0 


(84) 


Knowing the vectors h, v, 0, and'0, the various dot products needed 
to evaluate (23) and (24) can be expressed as 


A A 

0 ’ h = sin0 COS0 sinil; + s1n0 cosif; 


A A 

0 • h = COS0 simf; 


(85a) 


(85b) 
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6 • V = COS0 cos6 (- cost|j sin0 sine - sini|j cos9) 

. « 2 
+ sin0 COS0 cos^i COS0 sine - sini() cos0 sin 0 (85c) 

0 • V = sin0 (cosijj sin0 sin9 + sinti) cos0) 

2 

+ cosij) cos 0 sine .(85d) 

Expanding and simplifying (85c) and (85d), we find that 

0 • V = - 0 • h = -sin^) COS0 (86a) 

/N ✓v 

0 • V = 0 • h = sim(; cos6 sin0 + cosi{j sine (86b) 

flie now have the needed dot products in the unprimed coordinate system. 
The incidence angle e" is found by the relationship 

tan (J - e") = H/P (87) 

and P from 

P =Jl + = Hy cot 0(1 + sec^0) (88) 

Substituting (88) into (87) yields 
6 = tan" [/ cot^.0 ( 1 + s^c " 0 ) 


( 89 ) 
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Now that the incidence angle is known, let us next 
find the relationship between the primed coordinate system of the 
antenna and the unprimed coordinate system representing the water. 
Referring to Figure 7, the primed coordinate system is rotated 
about the z-axis through the angle a + 3. From Appendix I, we know 
the transformation to be 


e = e' (90) 

0 = 0' + a + 3 (91) 

We can now express (27), (28)., and (29) as functions of a and 3. 



(a ,6) 


j j 

over 

water 


[h(0,0) • 0 Egp(e,0 - (a + 3)) 


+ h(e,0) • 0 Epp(e»0 - (a + 3))] T|^^^f,[0"(0>0)] sin0 de d 


[v(0,0) • 0 Egp(0,0 - (a + 3)) + V (0;0) 


over 

water 



(«>e) 



i J 
over 

earth 


(0,0 - (a + 3))]^ Tj^^y[0"(0,0)] sine de d0 
T|^g(es0) Gp(0,0 - (a + 3)) sine d0 d0 


(92) 

(93) 



(« 53 ) 



T(,,(e.0) 


over 

sky 


Gp(9,0 - (a + 3)) sine d9 d0 


(94) 



where Gp has been normalized so that its value over the entire 
solid angle is unity. 

In order to evaluate (92), (93), and (94), we need to 
find the limits of integration. The edges of the tank which are 
parallel to the z-axis, as illustrated in Figure 9, lie in con- 
stant 0 plane. The values of 0 which define these edges are indi- 
cated in Figure 10 as 0-j and 0^. The boom length is p and 


H = p cosa 

(95) 

The distances C and D are knov'n as 


C = K/2 - p sina 

(96a) 

D = W/2 + p sina 

(96b) 


The relationship between 0-|, C, and D are 


1 


0^ = tan" (D/H) 

(97a) 

02 = tan" (C/H) 

(97b) 


Substituting (95), (96a), and (96b) into (97a) and (97b), we find 
the limits of integration for 0 between the water and earth to be 

_i 

0-| = tan [(W/2 + p sina)/p cosa] (98a) 

02 = tan" [(W/2 - p sina)/p cosa] 


(98b) 
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Fig. 9. Wave Tank Configuration to Determine the 6 Limits of 
Integration for the X-axis Geometry. 


W/2 


lave Tank Scanning P‘ 
f Integration for ti 
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Referring to Figure 9, for each value of 0 between 
-02 and 0^ we need to find the value of 0^^ that defines the other 
edge of the tank. First, we find B by 

B = H tan0 (gg.) 

We express E as 

E =/h + (100) 

By combining (96), (99). and {100). we find E as a function of p, a, 
and 0 as 


E = p cosa sec 0 (101) 

It can be seen from Figure 9 that 

tan (| - e^„) = F/E . cot(e^„) (102) 

and for this square tank 

^ ^ (103) 

The angle is found from (101), (102), and (103) as 

®ew "" [{W/2)/o cosa sec 0] ( 104 ] 

For values of 0 not between -02 and 0, we will define 9 as 

c I GW 2 * 

The sky is integrated for values of 0 between ^ and for all 
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values of e. We can now rewrite (92), (93), and (94), showing 
the limits of integration, as 


f 

^awp " ] 


0^ tt/2 


[same as in (92)] sin0 d0 d0 


■®2 %(a,p,U,0) 


0, ./2 


[same as in (92)] sin0 d0 d0 


(105) 


~h ®ew(a,p,W,0) 


V2 0 


'aep 


(a,6) = 


ew(a,p,W,0) 

[same as in (93)] sin0 d0 d0 


(106) 


-rr/2 0 


3tt/2 tt/2 




tt/2 D 


[same as in (94)] sin0 d0.d0 


(107) 


where 6p has been normalized so that its value over the entire solid 
angle is 2. 

We need only integrate 0 from 0 to it/ 2 because the 

geometry is symmetrical about the x-y plane. It should be noted 

here that the geometry does not have to be synmetrical to obtain a 

solution to the problem. The tank need not be square but of any 

shape. We need only to know F as a function of B to find 0 as a 

ew 

function of 0. If the tank is not symmetrical about the x-y plane, 

then we need to integrate 0 from 0 to it. The value of 0 that 

ew 

describes the back edge of the tank can be found in the same manner 



as the for 6 between 0 and tt/2. 

By utilizing the rotation about the z-axis to 
represent the rotation of the support arm through the angle a 
and the scanning angle B, we obtain a much more powerful repre- 
sentation of the problem than is possible with the earlier geo- 
metry which utilized the z-axis perpendicular to the aperture of the 
antenna. Previously we found from (56), (57), {58}, and (59) 
and from (70). Equations (98a), (98b) and (104) are much 
simpler. The dot products for the first coordinate system (z-axis 
normal to antenna aperture) are functions of e, 0, a, and polari- 
zation (0g). By rotating about the z-axis, the dot products are 
only functions of 0 and 0 because the transformation between 6, 0 

I I 

and 0 , 0 is performed merely by the addition of a constant to 0. 
Figure 11a shows the vectors 0 and 0 on the surface of the water when 
the z-axis is directed straight down into the water surface. The 
vector 0 is equal to v and 0 equal to h and all of the dot products 
are either 1 or C. When we rotate this coordinate system about 

A. j 

either the x or y axis, the 0 vector no longer lies in the plane 
of the water surface and the dot pro^cts become functions of 9*, 

t 

and the amount of angular rotation. The dot products would, 
therefore, need to be calculated for each rotation angle. This 
is a consequence of the transformation of variables between 0, 0 

I I 

and 0, 0 shown in transcendental form in Appendix I. Figure 11b 
shows how the vectors 0 and 0 align on the surface for the second 
system (x-axis normal to antenna aperture). If this coordinate 
system is rotated about the z-axis and the vector e' and 0* were 

















plotted they would look exactly the same all the time. Rotation 
about the z-axis does not' affect the alignment of the vectors on 
the surface and consequently does not alter the dot products. 

In addition to all the above advantages, the rotation 
about the z-axis has another tremendous advantage. The integration 
with respect to 0 is in the form of a correlation, which can be 
evaluated fay the use of Fourier transform techniques. The inte- 
gration with respect to 0 will be executed by numerical integration. 

In order to reduce the computation time, we make the 
gain functions independent of a. This way the spectrum of the gain 
need only be found once and used for all values of a.' To do this 
we add a to 0 in the integrands of (105), (106), and (107) and then 
subtract a from the limits of integration creating 0^ = 0 - a. 

This yields 



(a.B) 


I [h{9>0o+ a) * 0 Egp(6,0o- B) + h(0,0o+ a) • 

- 02 “Ct ®ew(a,p,W,0o+ a) 


^0.| -a ir/2 



[v(0,?L+ a) . 0 Egp(0,!2!,- B) + vCe,0o+ a)- 
0ew(«»P»W,0o+ a) 


0 EjgpCo.Po- T^^yC0"C0,0o+ a)J s.tn0 do d0„. (108) 



(109) 



(a, 3) 


yjZ a 

J Gp(e,0„- 6) sine d8 d0„ 

-Tr/2-a Q 



(a .6) 


3Tr/2-a it/2 

’ 

Tbs(0.0o+ ct) Gp(0>0o“ B) sine de d0 

7r/2-a 0 


The following functions will now be defined: 


^bwb 


( 0 , 00 + a) 


■J*bwh^®"^’ fo»" - 02 * ^ 0o^ 01-01 

and ^ 0 i tt/2 


0 el siiwhere 



(e,0^+ a) 


"^bwv^®"^ for - 02 - a ^ 0,s 0^-a 
and < 0 < ir/2 


0 el sewhere 


Tbe(e,0o+ a) 


ot) for -tt/2 - a ^ 0^5 ir/Z-a 
and 0 < 0 < 0g^ 


0 elsewhere 


o 010) 


(llOa) 


(110b) 


(110c) 
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Tbg(8.0o+ a) for tt/ 2 - a < 0g< 3ir/2 - a 
and al 1 0 

Tbs(9.0o+ a) E (nOd) 

0 elsewhere 

Using the primed functions we can rewrite (108) •* (110) as 

2 tt 7r/2 
r 

" Uh(e,0o+ a) • 0 Egp(0,0„- g) + h(0,0o+ a) • 

o'. 0 ^ 

0 E0p(e,0o- 3)] [v(6>0o+ «) • 

0 Egp(0,0,- 3) + v(0,0o+ a) • 0 E0p(6,0<,- 3)]"^ 

a) } s-'n0 de d0„ (111) 

2 it - 7t/2 

■’'aep^'^’®^ = I &) sin0 d0 d0o (112) 

0 0"^ 

2ir fr/2 

Tasp(a,3) = t'^(0,0o+ «) Gp(0,0„- 3) sin0 d0 d0„ (Tl3) 

0 0”^ 

Expanding (111) and dropping the arguments for convenience, yields 
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2ir ir/2 


T,„p(a.B) - J J [(h . S) E,p . (h . 0)’ 


0 0 


/s ^ 


* 2(h ■ _e) (h • 0) Eep E(,p 


" ■'b,™ ^0p' t;„ 


+ 2(v ■ 5) (5 . 0) Egp E^p sine de d06 (114) 


Equations (112), (113), and (114) are all of the form 


2ir ir/2 


Tg(a, 3 ) = 


Tjj(0,0^+ 0 .) G(0,!3„- 6) sine d0 d0<, (115) 


0 0 


By performing the integration with respect to 0 as a summation, (115) 
can be written as 


■ Z-n 


Tg(a,3) = S W.sin0^ J Tj^(6i,0„+ a) G(0^,0„- 3) d0^ (116) 

1=1 


The variable 0 is sampled N times while 0-j = 0 and 0^^ = ir/2. The 
function W^- is the weighting factor for the particular numerical 
integration technique used. 

For any constant value of a, (116) can be evaluated as 

^ J * 

Tg(aQ, - 3) W.sine. F"'{Tf^(f). ^.} • 


(117) 
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where 

= the periodic Fourier transform of Tj^{6^,0o+ a^) 

* 

G(f). = the complex conjugate of the transform of G(0^. ,0^) 

_i - - 1 _ - 

Since F [A] + F [B] is equivalent to F" [A + B], computation 

time can be reduced by evaluating (117) as 
1 N * 

Tj^{f). G{f). W. sine.} (118) 



(119) 



■ 1 N . 

= F {S L , G„, W, sine,} 


,aep' 0 - 


^'bei '"pi "i ^"'^i' 


i=l 


( 120 ) 


^aspK’“ ~ "^bsi %i ^i sine^J 


021 ) 


.Given a, equations (119), (120), and (121) can be used to find 
^awp’ ^aep’ ^asp’ ’"^spectively, for all values of 3. The 

5 * 

transforms of the gain functions for each polarization, E^p- , 

^0pi’ ^ 0 p ^0pi’ computed only once, since the gain 

functions are never needed in the time domain. The function 6p. is 
equal to Egp.| + E^p.j and need not be computed separately. 

The problem has new been formulated in two coordinate 
system configurations. The first used the z-axis as the antenna 
boresight to conform with the circularity of the antenna^s power 
pattern. The second required the x-axis as the boresight. This 
allows the use of Fourier transforms to perform the integration and 
reduce the computation time when the T^'s for all values of B are 

a 

computed. The necessity of knowing the T^*s for all values of 3 

a 

will become apparent in the inversion process. 

D. The Gain Functions 

In order to use either the z-axis or x-axis analysis, one 
must know the radiation characteristics of the antenna. Figure 12a 
shows the antenna geometry when the z-axis is taken to be perpen- 
dicular to the antenna aperture and Figure 12b shows the geometry 
with the x-axis perpendicular to the aperture. For either 


REFROfiUCHM^ cyp W- 
ORIGINAL'PA@E 'iS'POQR 
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Fig.- 12a. Z-axis Antenna Geometry. 
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y 



Fig. 12b. X-axis Antenna Geometry, 
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geometry. 

the three-dimensional power pattern 

can be obtained by 

G(6,0) 

= G[e (6,0), 0 = 0]cos 0 (9,0) 



+ G[0 (6,0), 0 = |-]sin 0 (0,0) 

(122) 

where 




1 _i 

0 = cos {cos0 sin6) 

(123a) 


1 i 

0 = tan {-sin0 tan6} 

(123b) 


Equations (122) and (12a) - (123b) can be used to construe 
the three-dimensional gain patterns from principal plane measure- 
ments whenever they can not be determined analytically. 

E, Cross-polarization 

Most practical antennas possess what is usually referred 
to as a cross-polarization pattern. For example, referring to 
Figure 12a, if there is no cross-polarization pattern, then in the 

I 

E-plane (0 =0) there would only be an Eq component and in the 

I 7T 

H-plane (0 = ^ ) only an E^ component. Any radiated component 
which is orthogonal to the principal polarization is usually 
referred to as cross-polarization. Having no cross-polarization in 
the principal planes does not insure no cross-polarization in any 
other plane, as is demonstrated for a reflector system by 
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Silver [13], In this investigation, the cross-polarized 

pattern will be assumed to have the same shape as the principal 

pattern. With this assumption, the antenna temperatures with 

cross-polarization and can be related to the antenna 

temperatures with no cross-polarization (T and T .) by 

sv 3n 

"^av = ‘av^^'' + ^ah * CROSS/0 + CROSS) (124a) 

‘‘‘ah " "^ah/^l + ^av ’ ^ROSS/O + CROSS) (124b) 

where CROSS refers to the fraction of power in the cross-polarized 
pattern. For example, if there is a -20 db cross-polarized 
pattern, then CROSS is 0.01. Equations (124a) and (124b) are exact 
if the shape of the cross-polarized pattern is the same as the 
principal pattern. 
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MI. Inversion 

A. Two-dimensional Approximation 

So far, we have only concerned ourselves with the direct 
problem, that of finding the antenna temperature T,(a,6) from the 
brightness temperature Tj^(e,0), Let us now approach the inverse 
problem of finding Tjj(e,0) from Tg(a,B). This inversion problem 
for the wave tank geometry was first approached by Fisher [7], 
Referring to Figure 7, Fisher used a two-dimensional approximation 
to represent the wave tank system. Assuming that the antenna maxi- 
mum (boresight) is directed only along the O = J plane, the water is 
scanned along the line Lp. The two-dimensional approximation assu- 
mes that most of the energy of the antenna is within the major lobe 
and only integration along the 0 = plane is necessary. With this 
approximation (1.11), (112), and (113) can be "reduced considerably. 
The dot products can be expressed as 


8 ’ h(|,0o+ a) = 0 • v(|,0o+ a) = 1 025a) 

0 • h(f,0o+ a) = 0 • v(|-,0o+ a) = 0 (125b) 

Renormalizing the gain functions so that 


0 


Gp(|,0)d0,= 1 


(126) 
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•we can write (111), (112), and (113) as 

2ir 

* 

“ J Eeplf-H.- e) “) 

0 

Zu 

* j B) a) it 

G 

2if 

Taep<“*e) “ ] Bp(f.0<r 6) t;^(| ,!9,+ a) dft 
2it 

Tasp(“-«) “ 1 ®p(f-M) 0.) it 

0 


(127) 

(128) 
(129) 


The limits of integration can be made 0 to 2ir since the primed 
brightness temperature functions are equal to zero in their res- 
pective regions to avoid overlapping. Assuming no cross-polari- 

2 

zation in the principal planes, E^p is zero for the horizontal scan 
and Egp is zero for the vertical scan. We can then write (127), 
(128) and (129) as 
2w 

Tawpt“-6) “ I «) d0„ 030) 

0 

2‘ir 

Taep{“-8) = J Sp(H.-B) “) dJl. (131) 

0 
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,2it 


Tjsp{a,8) = 


0 


Gp(0o-6) 




+ a) d0„ 


(132) 


which are of the same form as reported by Fisher [7] and Holmes [8]. 
With this two-dimensional approximation we can find the total 
antenna temperature by summing T T and T or T „ can be 

dWfJ dcp aSp ap 

found from one integral. To calculate directly we use the con- 
tinuous Tj^p{ 2 > 05 + a) which is 

Tbp(«°+ “) = T^„p(|,|9.+ a) + T'g(|,U„+ o) + T^3(|-.0„+ a) (133) 


We can therefore find T 
2n 


ap 


(a, 3) as 


Tap(a,e) 



Sp(0,-B) T|^p(0.t<») 


D34) 


0 

Equations (130), (131), (132), and (134) are all now in correlation 
form and can be evaluated, using Fourier transforms, as 


^awp^“o’~ 3) - F 
Taep'%-- 6) = r' 
^asp^“o’" S) = F 


«p • 


«P ■ TbeJ 


«P • 


(135) 

1136) 


(137) 
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• T^p} 038) 

Using (138) as an example, T^^p can be found from T^p by a simple 
division in the frequency domain followed by an inverse transform, 
or 


_ _ ^ 

Although (139) is a valid expression for the T^p 

expressed in (134), the inversion technique is extremely sensitive to 

error in T^p. Small errors in T^p can cause large oscillations in the 

inverted function for Tpp[8]. The instability of the equation can be 

explained in both the spatial and frequency domains. In the frequency 

domain, we note that, for the type of functions used for Gp, the 

high frequency components of its spectrum are small compared to the 

low frequency components. Due to the division by Gp, relatively 

small errors in the high frequency components of T7T can cause 

ap 

large errors in the corresponding components of T^, These errors 

cause high frequency oscillations in the spatial domain solution of 

Tpp in (139). This instability can also be explained by observing 

(134). If a high frequency sinusoid is added to T^p, the function 

T^„ will be nearly unaffected [3,4]. Therefore, (134) does not 
ap 

uniquely define T|^p for a T^p known with a moderate accuracy. 

The inversion of the Fredholm integral equation, of 
which (134) is one type, has been encountered and studied in the 
fields of aerosol detection, astronomical measurement 
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Interpretation, and spectral analysis where cause and effect 
situations are of Interest. Twomey [3] and Phillips [4] investi- 
gated the Fredholm equation of the first kind and were able to 
stabilize its inversion by employing matrix filtering techniques. 
Bracewell and Roberts [5] have reported an iterative restoration 
process which is particularly adaptable to our needs. The process 
introduced by Bracewell and Roberts and applied to a two-dimensional 

modelling of the wave tank geometry by Holmes stabilizes the 

* 

inversion by avoiding the direct division by G in (139). Writing 

— * * " ■ 

l/6p of (139) as 1/[1 - (1-Gp)], and then performing a series 

expansion [5] results in 


Tbp(f * »o> = F fT,„p tl Ml - Gp) 

* * 

+ 0 - Gp)" + (1 - Gp)* + ]) (140) 

* 

The infinite series expansion of 1/G converges provided that 

|l-Gpj < 1. For most antennas used in radiometry, their gain 

patterns are symmetrical, smooth varying functions which insure that 
— * 

G is always real and positive in the dominant frequencies. The 
^ — * 

maximum values of these Gp‘s will be the average value of the 
spatial domain functions. Since Gp is normalized by (126), the 
average value is 1 / 2 tt. The necessary conditions to insure the 
convergence of the series are therefore met. As the series converges, 
(140) becomes equal to (139) and the presence of error in T., will 
again cause oscillations. Fortunately, these unwanted oscillations 
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mainly arise from the higher order terms in the series expansion 

— * 

of l/6p. By properly truncating the series we can obtain the 
smooth principle solution of the inversion. This inversion process 
can also be performed in the space domain. The first term of the 
series expansion assumes that T is the first approximation of 
Tjjp and each addition term represents a new approximation of T^^p, 
This restoring process can be interpreted as letting the values of 
Tpp be equal, respectively, to 


T - T 

bpo ap 

(141a) 

^bpl " ^bpo ^ ^\p “ ^bpo^ 

(141b) 

'*'bp2 " ''’bpl ^"^ap " *^p * ’’’bpl i 

(141c) 

^bpn " ^bp(n-l) ^ ^^ap “ ^p ^bp(n-l)^ 

(141d) 


where * implies correlation. The altered inversion procedure 
reduces to an iterative method, as indicated by (141a) - (141d)-, and 
will be referred to as restoration [5]. The second term in (141a)- 
(141d) is a correction factor which is added to the values of the 
previous restored brightness temperature to obtain the newly 
created function. 

Holmes [8] has tested the restoration process with the 
two-dimensional simulation of the wave tank problem. He has calcu- 
lated Tgp(aQ,6) from the semi -empirical brightness temperature 
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models, added errors to represent measuring inaccuracies, and 

— * 

restored the data using a truncated series to represent 1/Gp in 
(140). Errors that caused large oscillation by direct inversion 
(139) did not create a stability problem in the restoration method. 

B. Three-dimensional Inversion 

Let us now investigate the problem of inverting the 
data using the three-dimensional model. As seen in (114), is 

» t 

dependent on both and and the same is true of T^^^, The 
equations are coupled and inversion is more complicated than in the ' 
two-dimensional case. At this point, we should review the goals of 
our inversion and what we will be given to obtain the necessary 
information. The desired results will be to obtain Tj^^^hC©”) and 
Tbwv(9")> where 0 " is the incidence angle, given Taf,(%»B) and 
Tj^v(%»3). To accomplish this, we will first need to make esti- 
mates of and (0,0) to calculate T^ep^“o’^^ ''^asp^%’^^ 

The estimated T,„^ and T^^^ will then be subtracted from the total 
aGp asp 

antenna temperature to find and In turn, the 

T .u(a.,,B) and T,,,,,(a„,3) functions will be inverted to restore 
awn 0 awv 0 

T^v,n(0") and 

At this point, we must make some assumptions about either 
the environment or the antenna in order to estimate T^g and T^g. 

The sky brightness temperature (T^g) is not extremely critical. 

In the calculations the brightness temperature of the 
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sky was assumed to be only a function of 0. The functional 
variation of along the 6 - plane as a function of 0 can be 
approximated by either the empirical sky model or by the T 

ap 

measured through the sky. Using these approximations of T^g(0) for 

all values of e is a good representation of the hemispherical 

brightness temperature profile except for the values of 0 near 0 

and TT. Since these areas are only seen by the sidelobes of the 

antenna, this error causes negligible error in T* . 

asp 

The brightness temperature of the earth is a little more 
critical since it borders the water. Since, for the earth, there is 
no functional relationship between the brightness temperature in 
the e = J plane and the other e planes, we assume that T^g is only 
a function of 0. T^p can be used as an approximation of T^^ except 
for the values of 0 which put the boresight near the wave tank. 
Referring to Figure 10, we define a new function tJJ( 0) as 


T,{fl) - T^p(e = 0 , - a 4 Ij) 

0, 4fo>0> 0 

(142a) 

fb(0) = T,p(S - 0 - a) 

— > 0 > 0 + — 
2 ^ ^1 90 

(142b) 

= 0 

0 > 2: 
2 2 

(142c) 

Tb(0) ■ Tgpfc - 0 - a) 

2tr - 02 r ;> 0 I. |2L 

(142d) 

Tb<0) - T,p(B - -02 - a - H_) 

27T 0 > 2TT •• 02 ^ 

(142e) 
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Using 6^^ as expressed in (104) and t|^( 0), we can form T^g(050) as 


Tbe(s>0) = t”(0) 

ew 

(143a) 

T^g(e,n) = 0 

Tf 6 > 

ew 

(143b) 


Equation (120) can now be used to find T (a ,3), 

aep 0 

Referring to Figures 8 and 10, the x‘-axis (antenna 
boresight) intersects the plane of the water surface at incidence 
angles between 0 and 0-| if 02 is positive or between *^2 0.j if 

02 is negative. Knowing the for these incidence angles is 

sufficient to compute the l^^^^p's over nearly the entire water sur- 
face of the tank, since the incidence angle can be computed exactly 
over the entire water surface with the use of (89). The only points 
which have incidence angles not in these ranges are the corners of 
the tank that are the farthest from the antenna. The brightness 
temperature functions of the water can be interpolated for these 
points. Given Tj^^^(| ,0) and brightness temperature 

all over the water surface can be found and, from (119), Tg^j^(a^,3) 
and T^y^,v{%^3) can be calculated. 

The restoration process used by Holmes [8] will now 
be applied to the three-dimensional problem. The first approxi- 
mation of and Tj^^^y(|-»0o) > Tor the particular range of 

incidence angles needed, will be Tg^^(3) and Tg^^(3), respectively. 
These first approximations will be called T^^^^(0J and Tj^^^^(0J, 
and from them we can find Tg^p^^( 3 ) and T^^y^(3) through (119). The 
difference between Tg^f^^(3) and T^^j^( 3 ) will be defined as ER^^( 3 ). 



REPR0!)U<5mUllf W' 'fM 
ORIGINAL'' PA©E IS POOS 

Similarly, the difference between and will be 

called ERy^(3). 

A second approximation of and 
^bwh2’ >^espectively) needs to be found from ^ and 

ER^2* Since the brightness temperatures and antenna temperatures 
are coupled, the following algorithm will be used to restore 


^bwh ^bwv' 

’’’bwhZ " "^bwhl ^^'"l ’ ^‘^hl '^'^2 ■ ^'^vl 

^bwv2 " '•'bwvl ' ^'^hl ^ ^^4 ‘ ^^1 ’ 

or in general 

"^bwh(n + 1) " "^bwh(n) '^^1 ' ^‘^h(n) ‘^''2 ' ^^(n) 

^bwv(n + 1) " \wv(n) '^''3 ’ ^*^h(n) ^^^4 ‘ ^*^v(n) (145b) 


The terms WF-j, WF2, WF^, and WF^ are weighting functions. 

To determine the weighting functions for any value- of 
6, we find what percentage of the power of the antenna incident on 
the water surface aligns v/ith the horizontal vectors and what per- 
centage aligns with the vertical vectors. If x% of the antenna's 
power picks up Tj^^^ while measuring Iay^jf^(3Q)> then x% of the error 
in the approximation of Tj^^^ at that particular incidence angle can 
be corrected by ERh^(3o). The functions WFp WF2, WF3, and WF^ are 
defined to be the percentage alignment of Gj^, G^, and G^ in the 



h, V, h, and v directions, respectively, for each value of 3. 

For example, the weighting function WF, would be found as 

» P * 

[h • 0 Egh + h • 0 sine de d0„ 

oyer 

WF^(3) - (146) 

sine de d0o 

over 

water 

All the functions in (145a and b) have now been defined and the 
restoration can be performed. 

The three-dimens"'Onal restoration acknowledges the 
influence of the following: 

1. the vector misalignment off the principal axes of the antenna, 

2. any cross-polarization in the antenna pattern, 3. the entire 
three dimensional environment, and 4. the true two variable power 
pattern of the antenna. The previously mentioned stability charac- 
teristics of the restoration procedure are also retained. 



IV. Computations and Results 

A. Finite Wave Tank 

The modeling of the interaction between the waye tank 
environment and the radiometer antenna has now been described in 
several different ways. In order to establish the validity of the 
various methods, computations were made with each of the methods 
while viewing identical environments with the same antenna. The 
results of these computations will establish the accuracy of the 
three-dimensional vector formulations and also the shortcomings' of 
the two-dimensional scalar method. 

Given the brightness temperature characteristics of the wave 
tank environment, the antenna temperature for the horizontal scan 
Tg^ and the vertical scan T^^ can now be predicted with four 
different computer programs developed during the course of this 
investigation.. The first of these is a computer program based on 
the three-dimensional analysis when the z-axis is normal to the 
radiometer antenna aperture, similar to the one developed by 
Beck [9], of the wave tank system. The integration with respect 
to the spherical coordinates e and 0 is numerically performed with 
the trapezoidal rule. This program is designed to make direct 
calculations of Tgj^(a,B=0) and Tg^(a,B=0) given the brightness 
temperature profiles of the water, ^bwv^®’^^’ 

earth Tjjg{0,0), and of the sky Tjjg(6,0). Since B==0 (no scanning) 
for all cases with this program, the antenna- is always assumed to be 
viewing the center of the tank at the incidence angle a. Secondly, 
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another three-dimensional computer program has been developed 
that takes the x-axis perpendicular to the radiometer antenna 
aperture and performs the calculations again with numerical 
integration. The integration with respect to 0 is done using a 
256 point midpoint rule and with respect to e using a 32 point 
Gaussian quadrature method. This program was written to yield 
^ah ^av various a s and 6=0, as was the first program. 

The third and most important program uses the same coordinate 
system alignment as the second program (x-axis perpendicular to 
the aperture) but has the capability to handle scanning in the 0 
direction. It uses a 32-point Gaussian quadrature method to 
integrate with respect to 0. However, the scanning of the antenna 
through the entire 360 range of the angle B and the integration 
with respect to 0 is handled simultaneously in the transform 
domain via the correlation form and a 256 sample point fast Fourier 
transform technique. The results of this program can be compared 
with those of the two previous programs for 6=0. It is the use of 
this program that enables the inversion (restoration) of the data. 
The last program that predicts the radiometer response is the two- 
dimensional approximation used by Fisher [7] and Holmes [8]. This 
formulation also uses fast Fourier transform techniques to carry 
out the integration for the scanning of the radiometer. 

Computed data of the same geometry using these four programs 
will be used to verify the analyses and computer programming of the 
equations. Several different comoarisons between them will he m^He 
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to validate the- techniques. If the results from the 2 -axis and 
,x-axis programs agree, the validity of both formulations as well 
as the accuracy of the programming and integration techniques will 
be established. When the results of the x-axis numerical inte- 
gration program match the output of the x-axis fast Fourier ■ 
transform program, they insure that the transform techniques is 
accurately performing the necessary integration. Finally, compari- 
sons of the two-dimensional approximation with the three-dimensional 
data provides criteria as to when it is mandatory to use the three- 
dimensional program to obtHn the desired accuracy. 

To test these various methods, the radiation characteristics 
of two different radiometer antennas will be used. Both of the 
antennas are pyramidal, corrugated horns with square apertures. The 
first horn has an aperture width of 12X and a total flare angle of 

O Q 

13 . The half- power beamwidth of this horn is approximately 6 . 

This is the antenna that is being used by NASA to take measurements 
in the wave tank system. In addition to the 12A aperture antenna, 
the response of the system using an 8X corrugated horn with a half- 
power beamwidth of 10 and a total flare angle of 19° will also be 
examined. The flare angles of these horns were designed so there 

O 

would be a 120 phase lag in the wave at the edges of the aperture 
as compared to the wave at the center. This particular phase taper 
in the aperture field creates a far-field pattern with no appreciable 
sidelobes or backlobes, which is desirable for radiometric measure- 
ments. Shown in Figures 13' and 14 are the principal plane patterns 







of the 8X and 12X horns, respectively. These will be used to 
carry out the necessary computations of the finite wave tank 
system. 

Before attempting any inversions, it will be desirable to 
perform some direct computations, calculating the antenna 
temperatures given the brightness temperatures, to validate the 
analyses and computer programming. Upon a successful evaluation 
of the formulation and programming, inversion (restoration^ of 
data will then be examined. 

1 • Direct Computations of Antenna Temperatures 

Physically, the wave tank is 14 feet square and the length of 
the boom that supports the antenna can be either 13 feet or 26 feet. 
In Table I, the predicted results for the 12X horn and the 13 foot 
boom are shown. The variables and T^^^ are the horizontal 
and vertical antenna temperature contributions from the water's 
surface, and T^^^ is the combined antenna temperature from the 
earth and sky. All three of the three-dimensional programs yield 
nearly identical results. With this antenna and boom length combi- 
nation, the two-dimensional formulation also yields fairly accurate 
results. 

In Table II, the predicted results for the 8X horn and the 
13 foot boom are shown. It is quite evident that the data from the 
two x-axis programs are almost identical and agreement with the 
z-axis results is very good. With the wider 8X horn, agreement 
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TABLE I 

Computed Antenna Temperatures for Finite Wave Tank System 

f 

(p = 13 feet. Antenna = 12X horn) 

{f = 10.69 GHz, T^ = 284^, S = 0 °/oo) 



0 

a=0 

a=20° 

a=40 

a=60 

a=80 

^awh 

109.07 

104.12 

89.17 

64.47- 

34.31 

^awv 

109.07 

114.48 

133.37 

176.71 

231.84 

^aes 

0.08 

0.09 

0.19 

1.28 

36,34 

^awh 

109.07 

104.13 

89,18 

64.48 

34.46 

T 

awv 

109.07 . 

114,46 

133.37 

176,77 

232,91 

^aes 

0.08 

0.10 

0,18 

1.20 

35.08 

\wh 

109.07 

104.12 

89.18 

64.48 

34,46 

^awv 

109.07 

114.46 

133.37 

176.77 

232.91 

^aes 

0.08 

0.09 

0.18 

1.20 

35.06 

T 

awh 

108.95 

104.00 

89.05 

64.46 

34.24 

T 

awv 

109.23 

114.59 

133.40 

176.74 

233,14 

^aes 

0.03 

0.05 

0.14 

1.15 

34,79 
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TABLE II 

Computed Antenna Temperatures for Finite Wave Tank System 
(p = 13 feet. Antenna = 8X horn) 

{f = 10.69 GHz, = 284^, S = 0 %o) 



O 

a=0 

a=20° 

O 

a=40 

a=60 

00 

o 

o 

3-D 

Z-AXIS N.I. 

^awh 

108.96 

103.98 

88.90 

63.52 

30,84 

^awv 

108.96 

114.40 

133.17 

173.03 

c. 

201.28 

^aes 

0.41 

0.46 

1.12 

7.19 

63.18 

3-D 

X-AXIS FFT 

T 

awh 

108.97 

104.01 

88.96 

63.58 

31.14 

T 

awv 

108.97 

114.38 

133.22 

173.23 

202.79 

^aes 

0.38 

0,46 

1.00 

6.90 

61.06 

3-D 

X-AXIS N.I. 

T 

awh 

108.96 

104,00 

88.95 

63.58 

31.14 

^awv 

108.97 

114.37 

133.21 

173.22 

202.79 

T 

aes 

0.38 

0.46 

1.00 

6.90 

60.98 

Q h- 
t Lt. 
CM U- 

^awh 

108.74 

103.77 

88.74 

63.54 

31.16 

^awv 

109.35 

114.72 

133.41 

173,30 

204.12 

T 

aes 

.16 

.12 

.79 

6.61 

59.42 

































































between the three- and two-dimensional computations is not as 
good. 

The results obtained using the 12X horn and the 26 foot boom 
are shown in Table III. Again the three-dimensional modelings yiel( 
results in agreement; however, the two-dimensional programming is 
not as accurate even. though the more efficient 12X horn was used. 
With the 26 foot boom, the angular limits of the wave tank are 
appreciably smaller and some of the main beam power spills over intc 
the earth. Since the earth is very warm, as compared to the water, 
an accurate modeling is needed at directions off the antenna bore- 
sight. The two-dimensional modeling can not provide this accuracy. 

To demonstrate how the greater accuracy of the three-dimen- 
sional formulation becomes imperative for antennas with wider main 
beams, the computed responses of the wave tank system with the 8X 
horn antenna and the 26 foot boom are shown in Table IV. By 
examining the data it is clear that, although all of the three- 
dimensional modelings agree well, the two-dimensional approximation 
is no longer an accurate method for predicting the radiometer 
response. Having now established the accuracy and necessity of 
the three-dimensional formulation, an examination of the inversion 
(restoration) procedure, as applied to the wave tank system, will 
be undertaken. 

2. Inversion (Restoration) Techniques for Antenna Brightness 
Temperature 


With the 26 foot boom and the 14 foot wave tank-, the 
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table III 

Computed Antenna Temperatures for Finite Wave Tank System 
(p = 26 feet. Antenna = 12A horn) 

(f = 10.69 GHz, = 284°K, S = 0 °/oo) 



■ 

n 


!■ 

ra 

m 

« 

»*— 1 

■31 

Q </) 
1 »-« 
ro X 
C 

1 

M 


■ 


88.58 

64.40 

25.09 


■ 

114.06 

132.32 

169,73 

173.30 

T 

aes 

.94 

1.16 

2.34 

11.95 

104.03 

3-D 

X-AXIS FFT 

iISfll 

108.76 

103.81 

88.66 

62.65 

25.04 

T 

awv 

108.77 

114.10 

132,49 

170.37 

173.31 

"^aes 

.93 

1.03 

2.02 

10.85 

104.12 

3-D 

(-AXIS N.I. 

^awh 

108.77 

103.81 

88.66 

62.65 

25.04 

T 

awv 

108.76 

114.10 

132.49 

170.37 

173.31 


^aes 

0.93 

1 .03 

2,02 

10.85 

104.10 

2-D 

FFT 

T 1 

T 

awh 

108.82 

103.83 

88.71 

62.83 

27.30 

T- 

awv 

109.09 

114.38 

132.78 

170.84 

188.43 


m 

.41 

m 

■ 

9.90 

85,90 

































TABLE IV 

Computed Antenna Temperatures for Finite Wave Tank System 

(p = 26 feet. Antenna - 8;;^ horn) 

(f = 10.69 GHz, = 284°K, S = 0 °/oo) 

a=0 a=20° a=40° a=60° a=80° 

106.94 101.72 85.60 57.48 18.66 


SS ''’awv 111-78 127.53 154.66 128.18 


2 12.90 36.46 150.02 


ira 


85.82 

58.04 

18.42 

106.92 

112.03 

128.00 

156.01 

126.85 

5.99 

6.45 

12.03 

34.07 

151.43 

106.93 

101.96 

85.82 

58.04 

18.42 

106.92 

112.03 

128.00 

156.01 

126,85 

5.98 

6.44 

12.03 

34.05 

151.36 

107.88 

102.78 

86.76 

59.23 

21,58 

108.42 

113.44 

129.77 

158.99 

147,34 

2.61 

3.33 

8.15 

28.60 

127,10 
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boresight of the antenna will intersect the edges of the wave 
tank when 3 = f 15 for a-0 . With the 13 foot boom, the angular 
limits of the wave tank are g = ± 28. s'* for a=0*’. For all other 
values of a, the angular space for viewing the wave tank is reduced. 
Since the antenna only possesses finite resolution, the range of 
incidence angles from which the brightness temperatures can be 
restored will be less than the angular limits of the wave tank. 
Therefore, in order to get continuous restored brightness 
temperature profiles, it will be necessary to combine the data from 

o o 

several values of a. The values of a that were chosen are 5, 10, 

2oI 3oI 4oI 5oI 55 °, 60 °, 65^ 70^ 75 °, and 80 °. 

In Figure 15, the results of a restoration process are shown. 
For the above mentioned values of' a, antenna temperature profiles 
were calculated, using the 12A horn antenna and the 13 foot antenna 
supporting boom, from the semi -empirical brightness temperature 
models of Stogryn [11]. Each a value only yields a limited range 
of incidence angles. The resulting antenna temperature profiles 
T^, the original brightness temperature profiles T^^, and the 
restored brightness temperatures using three iterations, are 

shown in the figure for vertical and horizontal polarizations. By 
combining the profiles for the different a s, nearly continuous 
curves have been formed. To improve the accuracy of the resulting 
curves, one could draw smooth curves through the 3=0 points of the 
data. These points yield fairly accurate results since they 
represent the system while the boresight is viewing the center of 
the wave tank. This has been done and the antenna, the restored 



TEMPERATURE ( K) 



Fig. 15, Continuous Incidence Angle Restoration Results 
for the Finite Wave Tank (Antenna = 12X horn, 
p = 13 feet, three iterations). 
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brightness and the empirical brightness temperatures for the 12A 
horn, 13 foot boom, and 3=0 data are shown in Figure 16. It can 
be seen from the plotted data in Figures 15 and 16 that the hori- 
zontal polarization results are improved by the restoration process 
for incidence angles greater than 60, and the vertical polarization 
inversion results are less accurate than the antenna temperatures 
for the larger angles. However, little can be inferred about the 
accuracy of the results for the incidence angles less than 60° for 
both polarizations. 

In order to get a more detailed look at the accuracy of the 
restoration process, some of the data from Figures 15 and 16 is 
listed in Table V. This table includes the total antenna tempera- 
ture T^, the restored brightness temperature the difference 

between T^ and the original brightness temperature Tj^, and the 
difference between T|jj,gg and Tj^. As can be seen, the restored 
brightness temperatures are always a better approximation of the 
true brightness temperature than the antenna temperatures for the 
horizontal polarization. Improvement is obtained in the vertical 
polarization case for all incidence angles up to and including 60. 
Although the differences between the antenna temperatures and the 
brightness temperatures were very small, the restoration process 
was still able to improve the results. 

Any instability in the solution will become evident as more 
restorations are taken. To show that the computations with three 
iterations are indeed convergent and are an improvement from those 
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(Antenna = 


Incidence 
Angle 


TABLE V 

Restored Antenna Temperatures for Finite Wave Tank 
with Three Restorations 

12a horn, p= 13' feet, f = 10.69 GHz, T^^^ = 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Angle 


VERTICAL POLARIZATION 



a 


109.18 


110.46 


114.56 


121 .92 


133.55 


151.25 


177.97 


219.07 


267.99 


109.12 


110.41 


114.46 


121.75 


133.17 


150.12 


176.15 


220.28 


283.05 


0.08 


0.07 


0.12 


0.19 


0.36 


0.68 


1.10 


2.11 


- 1.71 


o 

0 

109.14 

109.08 

0.04 


o 

10 . 

107.92 

107.88 

0.06 

0.01 

o 

■ 20 

104.22 

104.17 

0.06 

0.01 

0 

30 

‘98.02 

97.98 

0.06 

0.01 

40° 

89.36 

. 89.25 

0.10 

- 0.01 

O 

50 

78.44 

77.97 

0.34 

- 0.13 

O 

O 

65.69 

64.56 

0.98 

- 0.14 

o 

70 

58.15 

50.42 

8.20 

0.48 

o 

80 

69.54 

39.90 

31.57 

1.93 



bres " ' 


0.02 


0.01 


0.01 


0.02 


- 0.02 


- .45 


- .73 


3.32 
























































obtained with fewer restorations, the results of the inversion 
process for the 12X horn and the 13 foot boom with one restoration 
are listed in Table VI. By comparing Table V with Table VI, one 

o 

can see that, with the exception of the 80 incidence angle, more 
accurate results are obtained with the three-restoration process. 

O 

When a=80 , the angular limits of the wave tank (see Figure 7) are 

o o 

at 3 = + 3.5 and 3 = -11.3 . This is too small of an angular 
sector to expect accurate results. The restoration process is 
shown to be' convergent and to yield improved results. 

In Figure 17, the computed antenna, the restored brightness 
{three restorations), and the empirical brightness temperatures 
for the continous incidence angle data, utilizing the 8X horn and 

O 

the 13 foot boom, are shown. In Figure 18, the smoothed 3=0 curves 
for the same case are shown. As with the 12X horn and 13 foot data, 
an improvement can be seen in the horizontal polarization with a 
slight instability for the vertical polarization at incidence angles 

O 

greater than 60 . For a more accurate analysis of the data, 

Table VII has been included. From the table, one can see that the 
restored data yields a more accurate approximation than the 

O 

antenna temperatures for all angles listed, except 70 for the 
vertical polarization. Table VIII lists the one restoration 
results for the 8X horn and the 13 foot boom. Again one can con- 
clude that the restoration process has proven to be convergent and 
util itarian. 

The restoration process has also been applied to computed 
antenna temperatures for the 12X horn and the 26 foot boom. For 
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TABLE VI 


Restored Antenna Temperatures for Finite Wave Tank 
with One Restoration 

(Antenna = 12X horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S =• 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


VERTICAL POLARIZATION 



109.18 

109.12 

110.46 

110.41 

114.56 

114.46 

121.92 

121.77 

133.55 

133.30 

151.25 

151.17 

177.97 

179.15 

219.07 

224.60 

267.99 

271.27 


'a “ 'b 




0.07 


0.12 


0.19 


0.36 


0.68 


2.11 


- 1.71 


Incidence 

Angle 

. Ta 

T 

bres 

— I 
1 

cr 

Ves - Tt 

0 

0 

109.14 - 

109.08 

0.04 

- 0.02 

0 

10 

107.92 

107.87 

0.06 

0.01 

o 

20 

104.22 

104.17 

0.06 

0.01 

o 

30 

98.02 

97.98 

0.06 

0.02 

o 

40 

89.36 

89.30 

0.10 

0.04 ■ 

0 

50 

78.44 

78.27 

0.34 

0.17 

O 

60 

65.69 

65.16 

0.98 

0.46 

o 

70 

58.15 

51.08 

8.20 

- 1.14 

00 

o 

o 

69.54 

38.90 

31.57 

O'. 93 


^bres ’ ^ 
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TABLE VII 


Restored Antenna Temperatures for Finite Wave Tank 
with Three Restorations 

(Antenna = 8 a horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo} 


Incidence 

Angle 



HORIZONTAL POLARIZATION . 


"*^bres "^a “ 


bres " b 


.03 


0.01 


- 0.40 


- 0.16 


- 0.23 


- 0.19 


0.13 


o 

0 

109.34 

109.12 

0.24 

, o 
10 

. 108.14 

107.87 

0.27 

0 

20 

T04.47 

104.12 

0.31 

0 

30 

98.41 

97.80 

0.44 

o 

O 


89.03 

0.70 

50° 

80.31 

77.90 . 

2.22 

0 

60 

70.48 

64.84 

5.78 

0 

O 

71.81 

51.62 

21.86 

0 

80 

92.20 

38.86 

54.23 


VERTICAL POLARIZATION 



Incidence 

Angle 


T 

bres 

a b 

T - T 

bres b 

o 

0 

109.39 

109.17 

0.29 

0.07 

o 

10 

110.68 . 

110.40 

0.28 

0.00 

o 

20 

114.84 

114.37 

0.39 

- 0.08 

o 

30 

122.36 

121.39 

0.64 

- 0.34 

0 

o 

134.22 

132.57 

1.02 

- 0.62 

o 

50 

152.63 

149.97 

2,06 

- 0.60 

o 

60 

180,13 

177.74 

3.25 

. 0.87 

70° 

221.39 

■m 

4,43 

13.02 

O 

80 

263.84 

272.49 

- 5.86 

2.79 



















































































TABLE VIII 


Restored Antenna Temperatures for Finite Wave Tank 
with One Restoration 

(Antenna = 8X horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S=0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 
Angle 


Incidence 

Angle 


■^a-^b 


■a 


109.39 


110.68 


114.84 


122.36 


134.22 


152.63 


180.13 


221.39 


263.84 


VERTICAL POLARIZATION 


T. . 


bres 


109.16 


110.51 


114.68 


122.19 


134.30 


153.83 


182.87 


224.71 


252.18 


^a-^b 


^bres ” 


109.34 

109.12 

0.24 

0.02 

108.14 

107.96 

0.27 

0.10 

104.47 

104.31 

0.31 

0.15 

98.41 

98.19 

0.44 

0.23 

89.96 

89.68 

0.70 

0.42 

80.31 

79.02 

2.22 

0.92 

70.48 

65.87 

. 5.78 

1.16 

71.81 

51.22 

21.86 

1.28 

92.20 

37.38 

54.23 

- 0.59 




0.06 


0.11 


0.24 


0.46 


1.11 


3.26 


6.00 







































































three restorations, the continuous Incidence angle data is shown 
in Figure 19 and the smoothed g=0 data is shown in Figure 20. The 
x’s in Figure 20 indicate antenna temperatures that did not fit on 
the smooth curve. Improvement can be seen in the results except 
at the larger incidence angles for the vertical polarization. 

Table IX lists some of the data used in Figure 19 and 20. An 
improvement with the restored data can be seen at all angles for 

O * ^ 

the horizontal polarization and up to 60 incidence for the vertical. 
It should be noted that with the 26 foot boom the angular limits of 
the wave tank are 3 = + 4.2° and 3 = -17.0° for a=7oI and 3 = + 2.1° 

o o 

and 3 = -3.6 for a=80. To again demonstrate the convergence of the 
restoration process, the one-restoration computations are listed in 
Table X. Comparing the data in Tables IX and X again verifies the 
convergence and the need for the restoration process. 

The fourth and final antenna and boom combination is the 8X 
horn with the 26 foot boom. For three restorations, the continuous 
incidence angle data is shown in Figure 21 and the smoothed 3=0 
data is shown in Figure 22. For this case, the largest difference 
between the antenna and brightness temperatures is realized and 
the restoration process is needed the most. As can be seen in 
the figures, the restoration process works very well and yields 
an improved result for both polarizations at all incidence angles. 
Tabulating the data shown in Figures 21 and 22, one obtains 
Table XI and can see that the inversion process does recover the 
original water brightness temperatures and T^^^^ with good 
accuracy up to an incidence angle of- about 60°. To again check 
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Fig. 19. Continuous Incidence Angle Restoration 

Results for the Finite Wave Tank (Antenna = 
12X horn, p = 26 feet, three iterations). 



TEMPERATURE ( K) 



TABLE IX ' 


restored Antenna Temperatures for Finite Wave Tank 
with Three Restorations 

1 Antenna = 12x horn, p = 26 feet, f = 10.69 GHz, T^, = 284°K, S- 0 °/oo) 

HORIZONTAL POLARIZATION 


T_ - T. 



Incidence 



•yCRTICAL POLARIZATIO 




a 


110.18 


110.94 


115.13 


122.65 


134.52 


153.70 


181.23 


227.85 


277.43 



0 

0 

.110.14 

109.04 

1.04 


o 

10 

108.41 

107.74 

.54 

- .13 

. o 

20 . 

104.84 

103.88 

.68 

- .28 

o 

30 

98.90 

97.56 

.94 

- .40 

0 

O 

90.68 

88.80 

1.42 

- .46 

o 

o 

82.69 

77.83 

4.60 

- .27 

o 

60 

73,50 

65.17 

8.79 

.46 

O 

o 

88.07 

52.86 

38.12 

2.92 

o 

O 

CO 

129.16 

40.39 


2.42 


'*'bres 

■'a -Tb 

T - T 

bres b 

109.08 

1.08 

- .02 ■ 

110.25 

.54 

- .15 

114.08 

.68 

- .37 

121.12 

.93 

- .60 

132.42 

1.33 

- .77 

150.24 

3.13 

- .33 

179.19 

4.35 

2.32 

231.94 

10.89 

14.98 

• 231.69 

7.73 

11.99 









































































TABLE X 


Restored Antenna Temperatures for Finite Wave Tank, 
with One Restoration 


Incidence 
ng 


{Antenna = 12; 

V horn, p = 26 feet, f = 10.69 GHz, T^^ = 
HORIZONTAL POLARIZATION 

284 K, S = 0 ° 


Ta 

^bres 

■'a "T, 

BH 




0° 

110.14 

108.98 

1.04 

- 0.12 

o 

10 

108.41 

108.16 

0.54 

0.29 

o 

20 

104.84 

104.58 

0.68 

0.42 

o 

30 

98.90 

98.64 

0.94- 

0.68 

o 

0 

90.68 

90.30 

1.42 

1.04 

Q 

50 

82.69 

79.97 

4.60 

1.87 

O 

60 

73.50 

66.96 

8.79 

2.26 

0 

o 

r-.. 

88.07 

51.88 

38.12 

1.94 

o 

80 

129.16 

34.99 

91.19 

- 2,98 


a 


110.18 


110.94 


115.13 


122.65 


.134.52. 


153.70 


181.23 


227.85 


.43 


VERTICAL POLARIZATION 




bres 


109.02 


110.72 


114.97 


122.71 


135.09 


155.02 


184.29 


223.46 


242.77 




3.13 


4.35 


10.89 


7.73 




- 0.08 


0.32 


0.52 


0.99 


1.90 


4.45 


7.42 


6.50 


- 26.92 
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Fig. 21. Continuous Incidence Angle Restoration 

Results for the Finite Wave Tank {Antenna = 8X 
horn, p = 26 feet, three iterations). 
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TABLE XI 

Restored Antenna Temperatures for Finite Wave Tank 
with Three Restorations 


(Antenna = 8x horn, p= 26 feet, f = 10.69 GHz, T^ = 284° K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Anole 

L 

^bres 

\ - Th 

cL b 

^bres “ ^b 

O 

0 

115.43 

106.80 

6.33 

- 2.30 

o 

10 

111.36 

107.24 

3.49 

- .63 

o 

20 

108.41 

103.60 

4.25 

- .57 

o 

30 

103.90 

97.67 

5.94 

- .29 

o 

40 

97.85 

89.47 

8.60 

.21 

O 

50 


80.30 

17.43 

2.21 

O 

60 

92.11 

67.61 

27.40 

2.90 

O 

70 

124.80 

53.63 

74.85 

3.68 

O 

80 

169.85 

36.87 

131.88 

- 1.08 



VERTICAL POLAR] 

'ZATION. 


Incidence 

Angle 

. T, 

^bres 

"a -^b 

^bres ~ 

. O 
0 

115,45 


6.35 

- 2.28 

0 

10 

113.85. 

109.80 

3.45 

- .59 

20° 

118.48 

113.84 

4.03 

- .61 

o 

30 

126.98 

121 .43 

5,26 

- .30 

0 

40 

140.03 

133.84 

6.84 

.64 

o 

50 

161.31 

156.76 

11.25 

6,19 

o 

60 

190.08 ‘ 

188.16 

13.20 

11.29‘ 

O 

70 

238.86 

229.40 

21.89 

12.44 

O 

80 

278.28 

254.82 

8.58 

- 14.88 








































the convergence of the restoration process, the single restoration 
data for the 8X horn and the 26 foot boom is listed in Table XII. 
Comparing Tables XI and XII, the data shows that the process is 
convergent with the exception of the larger incidence angles. 

Now that the restoration process has been investigated with 
error-free data, antenna temperatures with added errors will be 
examined to determine their effect on the inversion. To include 
error in the data, the antenna temperature functions will be 

- o 

sampled every 5.6 and then interpolated between these points to 
obtain the 256 needed data points (sampling every 1.4 ). The first 
interpolation method to be used is a routine called SPLINE which 
was provided by Squire [14]. The SPLINE program uses a polynomial 
to represent the function between the sample points. The poly- 
nomials are formed so that the derivative of the interpolated curve 
is continuous at the sample points. In addition to the SPLINE 
interpolation method, linear interpolation was also used. The 
results obtained using these two interpolation methods are listed 
in Tables XIII through XXVIII. These tables show the data obtained 
with (a) 12X antenna, 13 foot boom; (b) 8X antenna, 13 foot boom; 
(c) 12X antenna, 26 foot boom; and (d) 8X antenna and 26 foot boom. 
For each antenna and boom combination and particular interpolation 
method used, a graph is included of the results obtained with the 
optimum number of restorations using the 8=0 data. These graphs 
comprise Figures 23 through 30. 

For the 12X horn and the 13 foot boom, examination of the 
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TABLE XII 

Restored Antenna Temperatures for Finite Wave Tank 
with One Restoration 

(Antenna = 8X horn, p = 26 feet, f = 10.69 GHz, T^ = 284\, S = 0 %o) 
HORIZONTAL POLARIZATION 


Incidence 

Angle 


T 

bres 

Ta -Tb 

^bres “ ^b 

O 

0 

115.43 

108.69 

6.33 

- 0.41 

o 

10 

111.36 

109.81 

3.49 

1.94 

o 

20 

108.41 

106.56 

4.25 

2.40 

o 

30 

mmm 

100.86 

5.94 

2.90 

o 

40 

97.85 

92.55 

8.60 

3.29 

o 

50 

95.52 

81.88 

17.43 

3.79 

o 

60 

92.11 

67.88 

27.40 

3.17 

O 

70 

124.80 

48.69 

74.85 

- 1 .25 

O 

80 

169.85 

28.66 

131.88 

- '9,31 


VERTICAL POLARIZATION 


Incidence 

Angle 

T 

a 

^bres 

"■a -Tb 

^bres " ^b 

0 

0 

115.45 

108.71 

6.35 

- 0.39 

10° 

113.85 

112.48 

3.45 

2.08 

o 

20 

118.48 

117.35 

4,03 

2.91 

o 

30 

126.98 

125.76 

5.26 

4.03 

o 

O 

140.03 

138.75 

6.84 

5.56 

o 

50 

161.81 

158.61 

11.25 

8.04 

0 

O 

190.08 


13.20 

8,52 


238.86 

203.71 

21 .89 

- 13i25 

O 

80 

278.28 

196.29 

8.58 

- 73.41 
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TABLE XIII 


Restored SPLINE Interpolated Antenna Temperatures 

for Finite Wave Tank with One Restoration 
(Antenna = 12 a horn, p = 13 feet, f = 10.69 GHz, T^ = 284 ” k , S = 0 °/oo) 


Incidence 
gle 


Incidence 

Angle 


a 


109.18 


HORIZONTAL POLARIZATION 


T,. ... 


bres 


109 . U 


T. “ T. 
a b 



a 


109.21 


110.46 


114.56 


121 .92 


133.55 


151.25 


177.97 


219.07 


267.99 


VERTICAL POLARIZATION 


T. 



'b 


109.17 


110.40 


114.48 


121.81 


133.23 


151.16 


179.46 


224.09 


270.61 



o 

10 

107.92 

107.86 

0.06 

. - 0.01 

o 

20 

104.22 

104.20 

0.06 

0.04 

o 

30 

98.02 

98.05 

0.06 

0.08 

0 

40 

89.36 

89.18 

0.10 

- 0.08 

o 

50 

78.44 

78.23 

0.34 

0.13 

o 

60 

65.69 

66.19 

0.98 

1 .48 

o 

70 

58.15 

48.87 

8.20 

- 1,08 

o 

80 

69.54 

38.34 

31.57 

0,37 


1 — 

^bres “ ^b 

0.12 

0.08 

0.07 

0.00 

0.12 

0.04 

0.20 

0.08 

0.36 

0.04 

0.68 

0.59 

1.10 

2.58 

2.11 

7.13 

- 1.71 

0.91 
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. TABLE XIV 

Restored SPLINE Interpolated Antenna Temperatures for 
Finite Wave Tank with Three Restorations 

(Antenna = 12X horn, p = 13 feet, f = 10.69 GHz, T^ .= 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Anqle 


• T 

bres 

T, - T. 
a b 

"^bres ■ ^b 

O 

0 

109.18 

109.17 

0.08 

0.07 

10° 

107.92 

107.80 

0.06 

- 0.06 

o 

20 

104.22 

104.32 

0.06 

0.15 

o 

30 

98.02 

98.26 

0.06 

0.29 

o 

40 

89.36 

88.72 

0.10 

- 0.54 

o 

50 

78.44 

77.88 

0.34 

- 0,21 

o 

60 

65.69 

68.83 

0.98 

4.12 

70° 

58.15 

43.76 

8.20 

- 6.19 

0 

80 

. 69.54 

37.77 

• 31.57 

- 0.20 


VERTICAL POLARIZATION 


Incidence 

Angle 


^bres 

Ta -Tb 

T - T 

bres b 

0° 

109.21 

109.21 

0.12 

0.11 

10° 

110.46 

110.35 

0.07 

- 0.05 

O 

20 

114.56 

114.57 

0.12 

0.13 

0 

30 

121.92 

121.94 

0.20 

0.21 

O 

40 

133.55 

132.87 

0.36 

- 0.32 

50° 

151.25 

150.14 

0.68 

- 0.43 

60° 

177.97 

177.44 

1.10 

0.56 

0 

70 

219.07 

218.79 

2.11 

1.83 

O 

80 

267.99 

280.72 

- 1.71 

11.02 
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TABLE XV 

Restored Linearly Interpolated Antenna Temperatures 

for Finite Wave Tank with One Restoration 
(Antenna = 12 a horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


^bres 

^a-Tb 

■ 

^bres ~ ^b 

0 

0 

109.08 

109.02 

- 0.02 

- 0.08 

0 

10 

1G7..92 

107.91 

0.06 

0.04 

o 

20 

104.22 

104.21 

0.06 

0.04 

0 

30 

98.02 

98.01 

0.06 

0.04 

0 

40 

89.36 

89.23 

0.10 

- 0.03 

0 

50 

78.44 

77.62 

0.35 

- 0.48 

0 

60 

65.69 

63.48 

0.98 

- 1.23 

0 

70 

58.15 

43.21 

B.20 

- 6.73 

0 

80 

69.54 

32.39 

31.57 

- 5.58 



VE 

kTICAL POLARIZATION 

Incidence 
• Angle 

Ta 

^bres 

Ta-Tb 

^bres “ ^b 

0 

0 

109.25 

109.19 

0.15 

0.09 

0 

10 

110.46 

110.36 

0.07 

- 0.04 

o 

20 

114.56 

114.37 

0.12 

- 0.07 

0 

30 

121.92 

121.64 

0.20 

- 0.08 

0 

40 

133.55 

133.09 

0.36 

- 0.10 

0 

50 

151.25 

150.65 

0.68 

0.08 

0 

60 

177.97 

178.38 

1.10 

1.50 

O 

70 

219.07 

223.20 

2.11 

6.24 

0 

80 

267.99 

272.14 

- 1.71 

2.44 
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TABLE XVI 

Restored Linearly Interpolated Antenna Temperatures for 
Finite Wave Tank with Three Restorations 

(Antenna = 12A horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 %o) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

■ ’ - 1 , i’. „ 

T 

bres 

\ 

a D 

^bres “ ^b 

0° 

109.08 

109.02 

- 0.02 

- 0.08 

0 

O 

107.92 

108.00 

0.06 

0.13 

0 

20 

104.22 

104.32 

0.06 

0.15 

o 

30 

98.02 

98.14 

0.06 

0.18 

o 

40 

89.36 

89.35 

0.10 

0.09 

o 

50 

78.44 

77.09 

0.35 

- 1.01 

60° 

65.69 

61.16 

0.98 

~ 3,55 

O 

70 

58.15 

25.12 

8.20 

- 24.83 

O 

80 

•69.54 

19.15 

31.57 

- 18.82 


VERTICAL POLARIZATION 


Incidence 
, Angle 

Ta 

^bres 

Ta -Tb 

T - T 

bres b 

O 

0 

109.25 

109.20 

0.15 

0.10 

0 

10 

110.46 

110.26 

0.07 

“ 0.14 

ro 

o 

o 

114.56 

114.22 

0.12 

- 0.23 

0 

30 

121.92 

121.44 

0.20 

- 0.28 

40° 

133.55 

1.32.70 

0.36 

- 0.49 

O 

50 

151.25 

149.02 

0.68 

- 1.55 

O 

60 

177.97 

174.21 

1 .10 

- 2.67 

O 

70 

219.07 

215.87 

2,11 

- 1.09 

o 

O 

CO 

267.99 

287.06 

- 1.71 

17.36 







TABLE XVII 


Restored SPLINE Interpolated Antenna Temperatures for 
Finite Wave Tank with One Restoration 
(Antenna = 8A horn, p = 13 feet, f = 10.69 GHz, T = ZBA^K, S = 0°/oo) 


Incidence 

Angle 


^bres 

a b 

Ves - Tb 

0 

0 

109.38 

109.17 

0.28 

0.07 

o 

10 

108.14 

107.96 

0.27 

0.09 

0 

20 

104.47 

104,33 

0.31 

0.17 

o 

30 

98.41 

98.20 

0.'44 

0.24 

o 

40 

89.96 

89.63 

0.70 

0.37 

cn 

o 

o 

80.31 

79.15 

2.22 

1.05 

o 

60 

70.48 

65.99 

5.78 

1.28 

o 

70 

71.81 

50.74 

21,86 

0.80 

o 

80 

• 92.20 

37.11 

54.23 

- 0.86 


VERTICAL POLARIZATION 

Incidence 

Angle 

Ta 

^bres 

1 — 

1 

rd 

1 — 

mm 

O 

0 

1 

1 

109.21 

0.32 

0.11 

o 

10 

110.68 

110.51 

0.28 

0.11 

o 

20 

114.84 

HQEBli 

0.39 

0.25 

o 

30 

122.36 

122.19 

0.64 

0,47 

0 

o 

134.22 

134.28 

1.02 

1.08 

50" 

152.63 

153.89 

2.06 

3.32 . 

o 

60 

180.13 

182.90 

3.25 

6.02 

0 

70 

221,39 

224.61 

4.43 

7.65 

0 

80 

263.84 

252.07 

- 5.86 

- 17.63 
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TABLE XVIII 

: Interpolated Antenna Temperatures 

(lave Tank with Three Restorations 
Feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 



^bres 


B1 


109.20 

0.28 

0,10 


107.85 

0.27 

- 0.02 


104.23 

0.31 

0,06 


97.83 

0.44 

- 0.14 


88.73 

0-70 

- 0.53 


78.60 

2.22 

0,50 


65.70 

5.78 

0.99 


49.96 

21.86 

0.02 


38.06 

54,23 

0.09 


VERTICAL POLARIZATION 



T 

bres 

L “ T. 
a b 

IHBI 


109.24 

0.32 

0.14 


110.38 

0.28 

- 0.02 


114.45 

0.39 

0.01 


121.40 

0.64 

- 0.32 


132.40 

1.02 

- 0.79 


150. -30 

2.06 

- 0.27 






177.96 

229.62 

272,15 


3.25 

4.43 

5.86 


1.08 

12.66 

2.45 
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TABLE XIX 

Restored Linearly Interpolated Antenna Temperatures 

for Finite Wave Tank with One Restoration 
(Antenna = 3 a horn, p = 13 feet f = 10.69 GHz, T^^ = 284°K, S = 0 ^/oo) 


Incidence 
gle 



80 


Incidence 


HORIZONTAL POLARIZATION 

1 

T. 


109.30 

109.06 

0.20 

- 0.04 

108.14 

107.93 

0,27 

0.06 

104.47 

104.25 

0.31 

0.09 

98.41 

98.02 

0.44 

0.06 

89.96 

89.24 

0.70 

- 0.01 

80.31 

77.82 

2.22 

- 0.28 

70.48 

63,45 

5.78 

- 1.26 

71.81 

46.68 

21 .86 

- 3.26 

92.19 

33.41 

54.22 

- 4.56 


VERTICAL POLARIZATION 


T. - T,. 



109.47 

109.24 

0.38 

110.68 

MMm 

0.28 

114,84 

114.54 

0.39 

122.36 

121.94 

0.64 

.134.22 

133.90 

1 .02 

152.63 

153,11 

2.06 

18C.13 

181.97 

3.25 

221.39 

223.95 

4.43 

263.84 

253.20 

- 5.86 


T - T 
bres b 


0.14 


0.03 


0.09 


0.22 


0.71 


2.54 


5.09 


6.99 


- 16.50 
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TABLE XX 

Restored Linearly Interpolated Antenna Temperatures for 
Finite Wave Tank with Three Restorations 

(Antenna = 8X horn, p = 13 feet, f = 10.69 GHz, T^^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


^bres 


^bres " ^b 

O 

0 

109.30 

109.11 

0.20 

0,02 

o 

10 

108.14 

107.93 

0.27 

0.06 

o 

O 

CSJ 

104.47 

104.25 

0.31 

0.09 

0 

30 

98.41 

97.76 

0.44 

“ 0.20 

45. 

O 

o 

89.96 

88.44 

0.70 

- 0.82 

0 

50 

80.31 

75.28 

2.22 

- 2.82 

O 

60 

70.48 

57.69 

5.78 

- 7.02 

o 

0 

71.81 

35.77 

21.87 

- 14.18 

80° 

92.19 

26.16 

54.23 

- 11.81 

VERTI 

CAL POLARIZATION 


Incidence 

Angle 


"'^bres 

Ta -Tb 

^bres ■ ^b 

O 

0 

109.47 

109.30 

0.38 

0.20 

o 

10 

110.68 

110.28 

0.28 

- 0.12 

20° 

114.84 

114.14 

0.39 

- 0.30 

o 

30 

122.36 

120.96 

0.64 

- 0.77 

o 

40 

134.22 

131.74 

1.02 

- 1.46 

o 

50 

152.63 

148.18 

2.06 

- 2.38 

o 

60 

180.13 

174.95 

3.25 

- 1.93 

o 

O 

221.39 

227.67 

4.43 

10.71 

o 

80 

263.84 

277.49 

5.86 

7.79 






TABLE XXI 


Restored SPLINE Interpolated Antenna Temperatures for 
Finite Wave Tank with One Restoration 

(Antenna = 12X horn, p = 26 feet, f = 10.69 GHz, T^^ = 284 ^, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Angle 


Incidence 

Angle 




o 

0 

111.33 

110.76 


1.66 

o 

10 

108.41 

107.86 

0.54 

- 0.01 

o 

20 

104.84 

104.42 

0.68 

0.26 

O 

30 

98.90 

98.88 

0.94 

0.92 

o 

0 

90.68 

91.22 

1.42 

1.96 

O 

50 

82.69 

80.11 

4.60 

i 

2.01 

o 

60 

73.50 

65.73 

8.79 

1.02 

o 

o 

88. 07 

49.69 

38.13 

- 0.26 

o 

80 

129.16 

34.05 

91.19 

-- 3.92 



VERTICAL POLARIZATION 




^a-Tb 




111.36 

110.79 

2.26 


110.94 

110.43 

0.54 

0.03 

115.13 

114.82 

0.68 

0.38 

122.65 

122.88 

0.93 

1.16 

134.52 

135.69 

1.33 

2.50 

153.70 

155.12 

3,14 

4.56 

181.23 

183.98 

4.35 

7.11 

227.85 

222.52 

10.89 

5.55 

277.43 . 

242.35 

• 7,73 

- 27.35 
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TABLE XXII 

Restored SPLINE Interpolated Antenna Temperatures 
for Finite Wave Tank with Three Restorations 
{Antenna = 12x horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S 


Incidence 
Angl e 


O 


HORIZONTAL POLARIZATION 


■^a 

T 

bres 

T, - T. 
a b 

111.33' 

m.58 

2.23 

108.41 

106.20 

0.54 

104.84 

103.22 

0.68 

98.90 

98.81 

0.94 

90.68 

92.72 

1.42 

82.69 

78.75 

4.60 

73.50 

62.30 

8.79 

88.07 

45.60 

38.13 

129.16 

37.28 

91.19 




Incidence 
gle 


VERTICAL POLARIZATION 


2.49 


- 1 .66 


- 0.94 


0.84 


3.46 


0.65 


- 2.40 


- 4.35 



"^bres " ’’’b 


111.36 

111.61 

2.26 

110.94 

108.79 

0.54 

115.13 

113.44 

0.68 

■ 122.65 

121.98 

0.93 

■ 134.52 

135.00 

1.33 

153.70 

150.96 

3.13 

131.23 

178.88 

4.35 

227.85 

228.98 

10.89 • 

277.43 

260.02 

7.73 



1.00 


0.26 


1.81 


0.40 


2.00 


12.02 


10.32 
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TABLE XXIII 

Restored Linearly Interpolated Antenna Temperatures 


for Finite. Wave Tank with One Restoration 

(Antenna = 12 a horn, p = 26 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo 

m 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

^a 

^bres 

T. - T. 
a D 

IIR8EE9i 

o 

0 

110.51 ■ 

109.46 

1.41 

0.36 

o 

10 

108.41 

107.80 

0.54 

- 0.07 

o 

20 

104.84 

■ 103.89 

0.68 

- 0.27 

0 

30 

98.90 

97.40 

0.94 

- 0.56 

o 

40 

90.68 

88.27 

1.42' 

- 0.98 

0 

50 

82.69 

74.52 

4.60 

^ 3.58 

o 

60 

73.50 

58.18 

8.79 

- 6.53 

O 

o 

88.07 

38.92 

38.13 

- 11.02 

o 

80 

129.16 

26.52 

91.19 

- 11.45 


VERTICAL POLARIZATION 


Incidence 
• Angle 

L 

T 

bres 

Ta 

T - T 

bres b ■ 

o'" 

110.67 

109.63 

1.57 

0.53 

10° 

110.94 

110.30 

0.54 

- 0.10 

o 

20 

115.13 

114.25 

0.68 ■ 


o 

30 

122.65 

- 121.58 

0.93 

- 0.14 

0 

40 

134.52 

133.48 

1.33 

0.29 

o 

50 

153.70 

151.64 

3.13 

1.08 

o 

60 

181.23 

180.33 


3.46 

o 

70 

227.85 

218.98 

10.89 

2.02 

• O 
80 

277.43 

241.96 

7.73 

■- 27.74 
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TABLE XXIV 


Restored Linearly Interpolated Antenna Temperatures for 
Finite Wave Tank with Three Restorations 

(Antenna = 12X horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 


Incidence 

Angle 


HORIZONTAL POLARIZATION 



T 

bres 

Ta -^b 

^bres “ ^b 

110.51 

109.94 

1.41 

0.84 

108.41 

107.49 

0.54 

- 0.38 

104.84 

103.10 

.0.68 

- 1.07 

98.90 

95.77 

0.94 

- 2.19 

90.68 

84.98 

1.42 

- 4.28 

82.69 

61.76 

4.60 

- 16.33 

73.50 

38.34 

8.79 

- 26.37 

88.07 

8.88 

38.13 

- 41,07 

129.16 

8.35 

91.19 

- 29.62 



VERTICAL POLARIZATION 


Incidence 

Angle 

L 

T 

bres 

^a-Tb 

O 

0 

110.67 

110.11 

1.57 

o 

10 

110.94 

109.76 

0.54 

o 

20 

115.13 

113.04 

0.68 

o 

30 

122.65 

119.28 

0.93 

o 

40 

134.52 

129.22 

1.33 

o 

50 

153.70 

140.55 

3.13 

o 

60 

181 -.23 

167.55 

4.35 

O 

0 

227.85 

217.12 

10.89 

o 

o 

CO 

277.43 

278.83 

7.73 


^bres ~ ^ 


1.01 


- 0.64 


- 1.41 


- 2.44 


- 3.97 


- 10.02 


- 9.33 
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TABLE XXV 

Restored SPLINE Intei'polated Antenna Temperatures 
for Finite Wave Tank with One Restoration 

(Antenna = 8X horn, p = 26 feet, f = 10.69 GHz, T^^^ = 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Anale 

Ta 

T 

bres 

— 1 
Q) 

1 

^bres - Tb 

O 

0 

115.30 

108.51 

6.20 

- 0.59 

o 

10 

111.36 

109.83 

3.49 

1.96 

o 

20 

108.41 

106.66 

4.25 

2.50 

o 

30 

103.90 

101.04 

5.94 

3.07 

O 

40 

97.85 

92.75 

8.60 

3.49 

O 

50 

95.52 

81.81 

17.43 

3.72 

O 

O 

92.11 

67.52 

27.40 

2.81 

o 

0 

124.80 

47.96 

74.85 

- 1,98 

o 

80 

169.85 

28.07 

131.88 

- 9.907 


VERTICAL POLARIZATION 


Incidence 

Angle 

L 

^bres 

- T. 
a b 

T - T 

bres b 

O 

0 

115.32 

108.53 

6.23 

- 0.56 

o 

10 

113.85 

112.49 

3.45 • 

2.10 

o 

20 

118.48 

117.43 

4.03 

2.99 

O 

30 

126.98 

125.90 

5.26 

4.17 

O 

40 

140.03 

138.90 

6.84 

5.71 

0 

50 

161.81 

158.59 

11.25 

8.02 

O 

60 

190.08 

185.27 

13.20 

8.40 

O 

o 

238.86 

203.44 

21.90 

- 13.52 

o 

80 

278.28 

196.18 

8.58 

- 73.52 
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TABLE XXVI 


Restored SPLINE Interpolated Antenna Temperatures 

for Finite Wave Tank with Three Restorations 

8a horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Angle 


Incidence 

Angle 

L 

T 

bres 

T. - T. 
a b 

^bres " ^b 

o 

0 

115.30 

■ 106.66 

6.20 

- 2.43 

o 

10 

111.36 

107.30 

3.49 

- 0.56 

o 

20 

108.41 

104.10 

4.25 

“ 0.06 

o 

30 

103.90 

98.72 

5.94 

0.76 

o 

40 

1 

97.85 

90.76 

8.60 

1.50 

o 

50 

95.52 

80.46 

17.43 

2.36 

o 

60 

92.11 

66.87 

27.40 

2.16 

O 

70 

. 124.80 

51.28 

74.86 

1.34 

o 

80 

169.85 

34.62 

131.88 

- 3,35 


VERTICAL POLARIZATION 


L 

T 

bres 

T. - T, 
a D 

■^bres - "^b 

115.32 

106.69 

6.23 

- 2.41 

113.85 

109.86 

3.45 

- 0.54 

118.48 

114.27 

4,03 

- 0.18 

126.98 

122.26 

5.26 

0.53 

140.03 

134.78 

6.84 

1.59 

161.81 

156.99 

11.25 

6.42 

190.08 

188.03 

13.20 

11.15 

238.86 

228.58 

21.90 

11,62 

278.28 

254.42 

8.58 

- 15.28 
















































































TABLE XXVII 


Restored Linearly Interpolated Antenna Temperatures 
for Finite Wave Tank with One Restoration 
(Antenna = 8 a horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 


HORIZONTAL POLARIZATION 


Incidence 

Angle 

^a 

T 

bres 

^a-Tb 

T - T 

bres b 

O 

0 

116.30 

109.53 

7.20 

0.44 

o 

1=0 

111.36 

108.64 

3.49 

0.78 

o 

20 

108.41 

105.08 

4.25 

0.91 

o 

30 


98.74 

5.94 

0.77 

o 

40 

97.85 

89.55 

8.60 

0.29 

o 

50 

95.52 

77.21 

17.43 

- 0.89 

o 

60 

92,11 

61 .88 

27,40 

~ 2.83 

o 

70 

124.80 

41.61 

74.86 

- 8.34 

O 

80 

169.85 

25.03 

131.88 

- 12.94 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

T 

bres 

- T. 
a b 

"^bres “ 

o 

0 

116.44 

109.67 

7.34 

0.57 

0 

10 

113.85 

111.30 

3.45 

0.91 

o 

20 

118.48 

115.92 

4.03 

1.48 

o 

30 

126.98 

123.89 

5.26 

2.17 

0 

40 

140.03 

136.44 

6.84 

3.24 

o 

50 

161.81 

155.61 

11,25 

5.04 

O 

60 

190.08 

182.42 

13.20 

5.54 

O 

70 

238.86 

201.23 

21.90 

- 15.68 

O 

80 

278.38 

196.41 

8.58 

- 73.29 
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TABLE XXVIII 

Restored Linearly Interpolated Antenna Temperatures for 
Finite Wave Tank with Three Restorations 


(Antenna = 8A horn, p = 26 feet, f = 10,69 GHz, T^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 
Angle • 

Ta 

^bres 

^a-Tb 

^bres “ "'"b 

0 

0 

116.30 

108.70 

7.20 

- 0.40 

0 

10 

111.36 

104.85 

3.49 

- 3.02 

o 

20 

108.41 

100.40 

4.25 

- 3.76 

o 

30 

103.90 

92.46 

5.94 

- 5.50 

40° 

97.85 

81.06 

8.60 

- 8.20 

O 

50 

95.52 

64.89 

17.43 

- 13.21 

O 

60 

92.11 

47.16 

27.40 

“ 17.55 

0 

70 

124.80 . 

27.05 

74.86 

- 22,90 

0 

80 

169.85 

21.46 

131.88 

- 16.51 


VERTICAL POLARIZATION 


Incidence 

Angle 

■^a 

T 

bres 

^a-^b 

T - T 

'bres 'b 

0 

0 

116.44 

108.83 

7.34 

- 0.27 

o 

10 

113.85 

107.35 

3.45 

- 3.04 

o 

20 

118.48 

110.68 

4.03 

- 3.77 

0 

30 

126.98 

116.81 

5.26 

- 4.92 

o 

40 

140.03 

127.37 

6.84 

- 5.82 

o 

50 

161.81 

147.13 

11.25 

- 3.44 

0 

60 

190.08 

178.35 

13.20 

1 .47 

0 

70 

238.86 

- 220.41 

21.90 

3.45 

c 

80 

278.28 

255.46 

8.58 

- 15.24 
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Fig. 23. Restoration of SPLINE Interpolated Data for 
the Finite Wave Tank {Antenna = 12X horn, 
p =13 feet, one iteration). 
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Fig. 24. Restoration of Linearly Interpolated Data for 
the Finite Wave Tank (Antenna = 12X horn, 
p = 13 feet, one iteration). 
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Fig. 25. Restoration of SPLINE Interpolated Data for the 
Finite Wave Tank (Antenna = 8X horn, p = 13 feet, 
three iterations). 
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Fig. 26. Restoration of Linearly Interpolated Data for 
the Finite Wave Tank (Antenna = 8X horn, 
p = 13 feet, one iteration). 
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Fig. 27. Restoration of SPLINE Interpolated Data for the 

Finite Wave Tank (Antenna = 12X horn, p = 26 feet, 
three iterations). 
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Fig. 29. Restoration of SPLINE Interpolated Data for 
the Finite Wave Tank (Antenna = 8X horn, 
p = 26 feet, three iterations). 
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Fig. 30. Restoration of Linearly Interpolated Data 

for the Finite Wave Tank (Antenna = 8X horn, 
p = 26 feet, one iteration). 



SPLINE interpolated data shows that the single iteration results 
are more accurate than those obtained with three restorations. 

With the interpolation error, three restorations can no longer be 
taken and still achieve accurate results. With one restoration, 

O 

improvement is achieved up to an incidence angle of about 60 . 

With linear interpolation, again the single iteration data is 
better than with three restorations and it is an improvement over 
the antenna temperatures. Due to the high frequency error in the 
spectrum of the linearly interpolated antenna temperatures, the 
inversion results, using three restorations, are less accurate than 
for the SPLINE routine. The accuracy of the results yielded with 
one restoration is about the same with either interpolation method. 

With the 8X horn antenna and the 13 foot support boom, the 
tabulation of the SPLINE interpolated data shows that three resto- 
rations yield more accurate results than one iteration, and the 
recovered brightness temperatures are almost always better approxi- 
mations than the interpolated ^tenna temperatures. Examination of 
the linearly interpolated data shows that additional restorations 
are not useful due to the previously mentioned high frequency error 
in the spectrum of the antenna temperatures. However, the single 
iteration case does yield improved results over the original 
antenna temperatures at the lower incidence angles for vertical 
polarization and at all incidence angles for horizontal polarization. 
Comparing the best case data for the two different interpolation 
techniques, one can see that with the 8X horn and 13 foot boom 
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i-umumrtcion, iKniNt interpolation Is superior. 

Next, the effect of interpolation will be examined for the 
12X horn and 26 foot boom computations. With SPLINE interpolation, 
the accuracy of the one- and three-iteration results appears to 
be very nearly the same but with a slight overall superiority for 
the three restorations. With this antenna and boom combination, 
the tank looks very narrow and the antenna, having a very narrow 
beam, creates a rapidly varying function as the tank is scanned. 

O 

The 5.6 sampling is not rapid enough in this case to let the 
SPLINE routine fit an accurate curve. With the rapidly varying 
functions involved with this antenna-boom combination, the high 
frequency error in the linear interpolation is most pronounced 
making the three restoration results much inferior than those of 
one iteration. The one restoration case does, however, yield 
improved results over the antenna temperatures. Due to the rapidly 
varying functions, the linear interpolation yields better results 
than does the SPLINE routine. 

Finally, the antenna temperatures for the 8X horn and the 
26 foot boom vary less rapidly than those yielded with the narrow 
12A horn and 26 foot support boom. Consequently, the restoration 
of the SPLINE interpolated data is more convergent and the three- 
iteration case does yield better results than the single restora- 
tion. With three restorations, results are yielded that are a 
considerable improvement over the interpolated antenna temperatures. 
With the linear interpolation, multiple iterations are not 
desirable and the best results are obtained with one restoration 



where the recovery process definitely yields improved results 
and should be used. For the 8X horn and the 26 foot boom, a 
slightly more accurate approximation of the water brightness 

temperature can be obtained with linear rather than SPLINE inter- 
polation. 

The restoration process has up to now been investigated using 
error-free and interpolated antenna temperatures. The process will 
now be examined with a random error added to the antenna 
temperatures. A Gaussian error with a mean of 0° and a standard 
deviation or 1 was added to the antenna temperature profiles for 
each value of a. The maximum error that was added to the antenna 
temperatures was approximately f 2.s". These profiles were then 
smoothed through the use of a Fortran subroutine named ICSSMU of 
the IBM IMSL library. The function with error and the standard 
deviation of the error is supplied to the subroutine which places 
a smooth cubic spline along the given set of data points. The 
subroutine can also interpolate between the data points. In this 
investigation, the subroutine was used to smooth the antenna 
temperatures that were known every 1.4^ after the random error 
had been added. Also, to show the combined effect of both the random 
error and interpolation, data with the random error was supplied to 
the subroutine at 5.6 intervals and the program was used to both 
smooth the antenna temperatures and interpolate to provide the 

O 

needed 1.4 sampling. Tables XXIX through XXXVI show the results 
obtained by smoothing the antenna temperatures with added random 
error and not interpolating. All four combinations of the 12X, 
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TABLE XXIX 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and One Restoration 

(Antenna = 12A horn, p = 13 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


^bres 

^a-Tb 

^bres " ^b 

0° 

108.94 

108.63 

- 0.16 

- 0.47 

10° 

108.48 

108.81 

0.61 

0.95 

ro 

o 

o 

104.78 

105.12 

0.62 

0.96 

30° 

•98.59 

98.96 

0.63 

0.99 

o 

O 

89.92 

90.26 

0.66 

1.00 

o 

50 . 

78.98 

79.25 

0,88 

1.15 

O 

60 

66.15 

66.11 

1.44 

1.40 

70° 

57.66 

49.73 

7.71 

- 0.21 

CO 

o 

o 

72.66 

43.72 

34.69 

5.75 


VERTICAL POLARIZATION 


Incidence 

Angle 


^bres 


T - T 

bres b 

o 

0 

108.98 

108.67 

- 0.12 

- 0.43 

10° 

111.01 

111,33 

0.61 

0.93 

20° 

115.11 

115.39 

0.66 

0.94 

O 

30 

122.47 

122.70 

0.75 

0,98 

0 

40 

134.09 

134.22 

0.90 

1.03 

O 

50 

151.77 

152.09 

1.20 

1.52 

o 

60 

178.44 

180.00 

1.57 

3.12 

70° 

219.46 

225.10 

2.50 

8.14 

80° 

268,90 

272.60 

- 0.80 

2.90 
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TABLE XXX 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and Three Restorations 

(Antenna = 12A horn, p = 13 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

''a 

T 

bres 

Ta-fb 

^bres " 

0°' 

108.94 

108.33 

- 0.16 

- 0.76 

10° 

108.48 

109.38 

0.61 

1.52 

O 

20 

104.78 

105.69 

0.62 

1.52 

o 

30 

98.59 

99.52 

0.63 

1 .56 

0 

o 

89.92 

90.76 

0.66 

1.50 

o 

o 

78.98 

79.58 

0.88 

1.48 

o 

O 

66.15 

66.58 

1.44 

1 .88 

70“ 

57.66 

47.60 

7.71 

- 2.34 

00 

o 

0 

72.66 

46.79 

34.69 

8.82 

VE 

RTICAL POLARIZATION 

Incidence 

Angle 

Ta 

T 

bres 

jD 
1 — ■ 

^d 
1 — 

T - T 

bres b 

0 

0 

108.98 

108.37 

- 0.12 

- 0.72 

o 

10 

111.01 

111.88 

0.61 

1.48 

' 0 
20 

T15.il 

115.93 

0.66 

1.48 

o 

30 

122.47 

123.23 

0.75 

1.51 

40° 

134.09 

134.62 

0.90 

1.43 

50° 

151.77 

151.62 

1.20 

1.05 

O 

O 

178.44 

177.66 

1.57 

0.78 

O 

o 

219.46 

220.85 

2.50 

3.88 

o 

O 

CO 

268.90 

- ■ ■ 

284.78 

- 0.80 

15.07 
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TABLE XXXI 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and One Restoration 

(Antenna = 8X horn, p = 13 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

Ta 

"^bres 

T, “ T. 
a b 

T - T 

' bres ' b 

0° 

109.21 

108.74 

0.11 

- 0.35 

10° 

108.51 

108.62 

0.64 

0.75 

20° 

104.85 

105.01 

0.69 

0.85 

30° 

98.76 

98.86 

0.80 

0.89 

0 

o 

90.30 

90.37 

1.04 

1-11 

0 

50 

80.46 

79.40 

2.37 

1.30 

O 

60 

70.04 

64.92 

5.33 

0.21 

o 

0 

73.17 

53.15 

23.23 

3.21 

o 

o 

CO 

95.16 

42.20 

57.19 

4.24 


.VERTICAL POLARIZAT 

ION 


Incidence 

Angle 


T 

bres 

Ta-Tb 

"*^bres ~ 

0° 

109.26 

108.79 

0.16- 

- 0.31 

10° 

111.04 

111.15 

0.65 

0.75 

20° 

115.21 

115.36 

0.76 

0.91 

30° 

122.71 

122.83 

0.99 

1.11 

40° 

134.56 

134.96 

1.37 

1.77 

0 

O 

152.86 

154.28 

2.29 

3.71 

60° 

180.22 

182.95 

3.34 

6.07 

o 

o 

222.10 

225.88 

5.14 

8.92 

00 

o 

0 

264.92 

253.94 

- 4.78 

- 15.76 
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TABLE XXXII 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and Three Restorations 


(Antenna = 8A horn, p = 13 feet, f = 10.69 GHz, T^^^ = 284°K, S = 0 o/oo) 

HORIZONTAL POLARIZATION 


-Incidence 

Angle 

Ta 

T 

fares 

Ta-^b 

^fares “ ^fa 

0° 

109.21 

108.28 

0.11 

- 0.82 

io“ 

108.51 

109.03 

0.64 

1.16 

ro 

o 

o 

_J 

104.85 

105.36 

0.69 

1.20 

0 

O 

CO 

• 98.76 

98.99 

0.80 

1.03 

45 :. 

o 

o 

90.30 

90.38 

1.04 

1.12 

o 

50 

80.46 

79.07 

2.37 

0.97 

0 

O 

70.04 

63.29 

5.33 

- 1.42 

70“ 

73.17 

54.00 

23.23 

4.06 

CO 

o 

o 

95.16 

46.49 

57.19 

8.52 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

T 

fares 

- T. 

a b 

"^^fares " 

0“ 

109.26 

108.32 

0.16 

- 0.78 

10“ 

111.04 

111.53 

0.65 

1.13 

20° 

115.21 

115.56 

0.76 

1.11 

30° 

122.71 

122.54 

0.99 

0.81 

0 

40 

134.56 

133.81 

1.37 

0.62 

0 

O 

152.86 

150.99 

2.29 

0.42 

o 

o 

180.22 

177.83 

3.34 

0.95 

70“ 

222.10 

231.74 

5.14 

14.77 

o 

O 

CO 

264.92 

275.36 

- 4.78 

5.66 
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TABLE XXXI Ii 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and One Restoration 

(Antenna = 12X horn, p = 26 feet, f = 10.69 GHz, T^ = 284 “k, S = 0 °/oo) 
^ HORIZONTAL POLARIZATION 


Incidence 

Angle 

L 

^bres 

L-Tb 

"•^bres “ ^b 

O 

. 0 

109.86 

108.42 

0.76 

- 0.68 

10° 

108.97 

109.18 

1.10 

1.31 

20° 

105.37 

105.56 

1.21 

1.40 

30° 

99.40 

99.63 

1.44 

1.67 

o 

O 

91.13- 

91.25 

1 .87 

1.99 

0 

50 

82.36 

79.33 

4.26 

1.24 

0 

O 

72.97 

65.68 

8.26 

0.97 

O 

0 

90.45 

54.91 

40.51 

4.97 

0 

o 

00 

135.40 

44.67 

97.43 

6.70 


VERTICAL POLARIZATION 


Incidence 

Angle 

L 

T 

bres 

' Ta-Tb 

^bres " ^b 



0° . 

109.89 

108.45 

0.79 

- 0.65 

o 

10 

111.50 

111,73 

1.10 

1.33 

o 

20 

115.65 

115.94 

1.21 

1.49 

30 

123.16 

123.68 

1.43 

1.96 

o 

40 

134.96 

135.99 

1.77 

2.80 

50° ■ 

153.65 

154.88 

3.08 

4.32 

O 

60 

181.24 

184.15 

4.36 

7.27 

o 

0 

228.78 

224.54 

11.82 

7.58 

CO 

o 

o 

279.24 

245.46 

9.54 

- 24.24 
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TABLE XXXIV 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and Three Restorations 

(Antenna = 12A horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 


HORIZONTAL POLARIZATION 


Incidence 

Angle 

T 

a 

T 

bres 

Ta-Tb 

^bres ” ^b 

0“ 

109.86 

108.10 

0.76 

HBBSI 

10° 

108.97 

109.39 

1.10 

1.52 

0 

20 

105.37 

105.51 

1.21 

1.34 

0 

30 

99.40 

99.36 

1.44 

1.40 

O 

o 

91.13 

90.82 

1.87 

1.56 

o 

O 

LO 

82.36 

77.12 

4.26 

- 0.98 

60° 

72.97 

62.85 

8.26 

- 1.86 

o 

70 

90.45 

56.27 

40.51 

6.33 

CO 

o 

a 

135.40 

54.16 

97.43 

16.19 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

^bres 

T. - T. 
a b 

T - T 

bres b 

o 

0 

109.89 

108.14 

0.79 

-■ 0.96 

10° 

111.50 

111.89 

1.10 

1.49 

20° 

115.65 

115.68 

1.21 

1.23 

30° 

123.16 

122.86 

1.43 

1.13 

40° 

134.96 

134.23 

1.77 

1.04 

0 

O 

153.65 

150.22 

3.08 

- 0.35 

o 

o 


178.93 

4.36 

2.05 

o 

O 

.228.78 

233.05 

11.82 

16.09 

CO 

o 

0 

279.24 

285.24 

9.54 

15.54 
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TABLE XXXV 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, No Interpolation, and One Restoration 

(Antenna = 8A horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

— 

Ta 

"*'bres 

■ — \ 
I 

—1 

cr 

"*^bres " 

0° 

114.58 

107.20 

5.48 

- 1.90 

10° 

111.64 

110,51 

3.77 

2.64 

tN3 

o 

0 

108.52 

106.92 

4.36 

2.76 

30° 

103.74 

100.66 

5.78 

2.70 

40° 

97.33 

91.56 

8.07 

2.30 

50° 

95.30 

81.09 

17.20 

2.99 

o 

60 

92.65 

68.25 

27.94 

3.54 

o 

70 

128.57 

54.50 

78.63 

4.56 

00 

o 

o 

174.54 

1 

36.31 

136.57 

- 1.66 


VERTICAL POLARIZATION 


Incidence 

Angle 


T 

bres 

T, - T. • 

cl b 

^bres " ^b 

0° 

114.61 

107.22 

5.51 

- 1.88 

10° 

114.13 

113.17 

3.73 

2.78 

20° 

118.61 

117.73 

4.16 

3.29 

O 

O 

CO 

126.90 

125.68 

5.18 

3.96 

o 

40 

139.70 - 

138.12 

6.51 

4.93 

O 

50 

161.73 

158.21 

11.17 

7.64 

o 

0 

190.37 

185.58 

13.49 

8.70 

o 

o 

240.58 

206.28 

23.62 

- 10.68 

o 

O 

CO 

280.25 

199.50 

10.54 

- 70.20 
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TABLE XXXVI 

Restored Antenna Temperatures for Finite Wave Tank with 


Random Error, No Interpolation, and Three Restorations 
(Antenna = 8X horn, p = 26 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

Ta 

IBHII 

\ -■'b 

^bres “ ^b 

O 

0 

114.58 

104.42 

5.48 

- 4.68 

o 

10 

111.64 

109.11 

3.77 

1.24 

o 

o 

CSJ 

108.52 

104.91 

4.36 

0.75 

o 

30 

103,74 

97.96 

5.78 

0.00 

O 

0 

97.33 

88.17 

8.07 

- 1.09 

50° 

95.30 

78.67 

17.20 

0.57 

O 

60 

92.65 

67.62 

27.94 

2.91 

. 70° 

128.57 

62.23 

78.63 

12,28 

O 

80 

174.54 

48.27 

136.57 

10.30 

V 

ERTICAL POLARIZATION 

Incidence 

Angle 

Ta 

T 

bres 

- T. 
a b 

T - T 

'bres 'b 

O 

0 

114.61 

104.44 

5.51 

“ 4.65 

0 

10 

114.13 

111.64 

3.73 

1.25 

• 20° 

118.61 

115.13 

4.16 

0.68 

o 

30 

126.90 

121.86 

5.18 

0.14 

0 

40 

139.70 

133.10 

6.51 

- 0.09 

o 

o 

161.73 

155,95 

11.17 

5.38 

60° 

190.37 

188.17 

13.49 

11.29 

0 

70 

240,58 

233.06 

' 23.62 

16.10 

00 

o 

o 

280.25 

259.57 

10.54 

- 10.13 


































8X horn antennas and 13 foot, 26 foot supporting booms are listed 
with one and three restorations. Figures 31 through 34- are graphs' 
using the B=0 points for each of the four cases and the optimum 
number of restorations. Tables XXXVII through XLIV and Figures 
35 through 38 show data analogous to Tables XXIX through XXXVI and 
Figures 31 through 34, respectively, but with interpolation and 
smoothing provided by the subroutine ICSSMU. 

With the 12X'horn and the 13 foot boom, the error that has 
been added is greater than the difference between the antenna 
temperatures and the brightness temperatures. Consequently, the 
restored brightness temperatures are not as good an approximation 
of the brightness temperatures as are the smoothed antenna 
temperatures. Multiple restoration makes the restored results 
inferior. For this antenna and boom length, these observations are 
valid for both the interpolated and uninterpolated data. 

For the 8X horn and the 13 foot boom data that contains the 
random error, multiple restorations are not desirable either with 
or without interpolation. When the antenna temperatures are not 
interpolated, some improvement at the larger incidence angles for 
horizontal polarization is achieved with one iteration. With 
interpolation, the smoothed antenna temperatures are more accurate 
than the restored results. 

Using the 12X horn 26 foot boom, and no interpolation,' the 
three-restoration results are better than those for one interation. 
With this antenna and boom combination, and no interpolation. 



TEMPERATURE { K) 


•280 



Fig. 31. Restoration of the Finite Wave Tank Data with 
Random Error and No Interpolation (Antenna = 
12A horn, p = 13 feet, one iteration). 
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Fig. 32, Restoration of the Finite Wave Tank Data with 
Random Error and No Interpolation (Antenna = 
8X horn, p = 13 feet, one iteration). 
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Fig. 34. Restonition of the Finite Wave Tank Data with 
Random Error and No Interpolation (Antenna = 
8X horn, p = 26 feet, three iterations). 
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TABLE XXXVII 


Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and One Restoration 

(Antenna = 12X horn, p = 13 feet, f = 10.69 GHz, T = 284°K, S = 

HORIZONTAL POLARIZATION 

Incidence I I I ^ 1 I 


Incidence 

Angle 


0“ 

109.42 

109.28 

0.32 

0.18 

o 

10 


108.90 

0.63 

1.04 

0 

20 

104.8*1 

105.18 

0.64 

1.01 

0 

30 

98.65 

99.07 

0.68 

1.11 

o 

40 

89.99 

90.57 

0.73 

1.31 

o 

50 

78.26 

78.46 

0.17 

0.37 

o 

60 

63.50 

62.34 

- 1 .21 

- 2.37 

0 

O 

61.06 

52.70 

T1.12 

2.75 

0 

o 

00 

77.78 

50.57 

39.81 

12.60 


VERTICAL POLARIZATION 





109.44 

109.30 

1 

0.34 

0.20 

111.03 

111.42 

0.63 

1.02 

• 115.13 

115.45 

0.69 

1,01 

1 122.51 

122.80 

0.79 

1 .08 

134.10 

134.38 

0.91 

1.18 

151.30 

151.52 

0.73 

0.95 

177.67 

178.80 

0.79 

1.92 

220.44 

226.19 

3.48 

9.22 

270.94 

275.13 

1.24 

5.43 
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TABLE XXXVIII 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and Three Restorations 


(Antenna = 12X horn, p = 13 feet, f = 10.69 GHz, T^^ = 284^, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


"'^bres 

"a -Tb 

^bres ” ^b 

0° 

109.42 

108.73 

0.32 

- 0.36 

o 

10 

108.50 

109.44 

0.63 

1.58 

20° 

104.81 

105.56 

0.64 

1.40 

30° 

98.65 

99.41 

0.68 

1.45 

o 

O 

89.99 

91.03 

0.73 

1.77 

o 

50 

78.26 

79.13 

0.17 

1.03 

60° 

63.50 

63.12 

- 1.21 

- 1 .59 

O 

0 

61.06 

47.93 

- 11.12 

- 2.01 

0 

CD 

CO 

77.78 

54.91 

39.81 

16.' 94 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

^bres 


^bres “ ^b 

0° 

109.44 

108.76 

0.34 

“ 0.34 

10° 

111.03 

111.94 

0.63 

1.54 

o 

O 

CM 

115.13 

115.84 

0.69 

1.40 

0 

30 

122.51 

123.14 

0.79 

1.42 

O 

o 

134.10 

134.72 

0.91 

1.53 

50° 

151.30 

151.10 

0.73 

0.53 

60° 

177.67 

176.30 

0.79 

- 0.58 

0 

70 

220.44 

221.48 

3.48 

4.52 

O 

80 

270.94 

287.72 

1.24 

18.02 
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TABLE XXXIX 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and One Restoration 

(Antenna = 8A horn, p = 13 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION ^ 


Incidence 

Angle 

Ta 

^bres 

"a-Tb 

Tbres - Tb 

0° 

109.52 

109.33 

0.42 

0.23 

10° 

108.55 

108.85 

0.68 

0.98 

20° 

104.85 

105.25 

0.68 

1.09 

30° 

98.65 

99.01 

0.69 

1.05 

0 

O 

• 89.84 

89.99 

0.58 

0.73 

0 

50 

79.00 

76.96 

0.90 

- 1.14 

o 

o 

69.13 

62.87 

4.42 

- 1 .84 

o 

O 

76.63 

• 58.06 

26.69 

8.12 

0 

o 

00 

98.74 

46.98 

60.77 

9.01 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

"'"bres 

Ta -Tb 

^bres “ ^b 

O 

0 

109.54 

109.36 

0.44 

0.26 

o 

10 

111.07 

111.36 

0.67 

0.96 

ro 

O 

0 

115.19 

115.55 

0.74 

1.10 

o 

30 

122.60 

122.91 

0.87 

1.18 

o 

40 

134.22 

134.65 

1.03 

1.46 

o 

50 

152.18 

153.11 

1.61 

2.55 

60° 

180.20 

182.76 

3.33 

5.88 

70° 

223.11 

227.37 

6.15 

10.40 

80° 

265.05 

253.82 

- 4.65 

- 15.88 
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TABLE XL 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and Three Restorations 


(Antenna = 8X horn, p = 13 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 


Incidence 

Angle 

Ta 

^bres 

Ta -Tb 

T - T 

bres b 

0° 

109.52 

108.98 

0.42 

- 0.12 

10° 

108.55 

109.52 

0.68 

1.65 

20° 

104.85 

105.91 

0.68 

1.74 

30° 

98.65 

99.68 

0.69 

1.71 

45. 

O 

o 

89.84 - 

90.60 

0.58 

1.34 

50° 

79.00 

75.89 

0.90 

- 2.20 

O 

O 

69.13 

59.59 

4.42 

- 5.12 

70°- 

76.63 

^ -- 

60.00 

26.69 

10.05 

O 

80 

98.74 

52.49 

60.77 

14.52 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

^bres 

^a-Tb 

^bres “ ^b 

o 

0 

109.54 

109.00 

0.44 

- 0.10 

0 

10 

111.07 

111.99 

0.67 

1.60 

20° 

115.19 

116.03 

0.74 

1.58 

0 

O 

CO 

122.60 

123.02 

0.87 

1.29 

40° 

134.22 

133.84 

1.03 

0.65 

50° 

152.18 

149.37 

1 .61 

- 1.20 

60° 

180.20 

177.19 

3.33 

0.32 

70° 

223.11 

233.63 

6.15 

16.67 

00 

o 

o 

265.05 

275.29 

- 4.65 

5.59 





TABLE XL I 


■ Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and One Restoration 
(Antenna = 12 \ horn, p = 26 feet, f = 10.69 GHz, T^^ = 284 “k, S = 0 °/oo) 


HORIZONTAL POLARIZATION 


Incidence 

Angle 


Incidence 

Angle 



Ta 

T 

bres 

T. - T, 
a b 


109.08 

107.33 

- 0.02 

- 1.77 

108.81 

109.54 

0.94 

1.68 

104.80 

105.18 

0.64 

1.01 

97.80 

97.67 

- 0.17 

- 0.29 

88.18 

87.22 

- 1 .08 

- 2.04 

80.84 

76.34 

2.74 

- 1.75 

74.02 

65.29 

9.32 

0.58 

97.44 

62.98 

47.49 

13.04 

140.57 

51.02 

102.60 

13.05 


a 


109.09 


111.32 


115.12 


121.84 


132.72 


!52.44 


181.29 


231.74 


283.17 


VERTICAL POLARIZATION 




"a'-'b 


107.34 


112.06 


115.58 


122.08 


132.90 


152.57 


183.41 


227.89 


250.76 


- 0.01 


0.92 


0.68 


0.11 


- 0.47 


1.87 


4.41 


14.77 


13.47 


r - T 
bres b 


- 1.76 


1.66 


1.13 


0.36 


- 0.29 


2.01 


6.53 


10.92 


- 18.94 
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TABLE XLII 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and Three Restorations 
(Antenna = 12X horn, p = 26 feet, f = 10.69 GHz, T^^ = 284 ’’k, S = 0 °/oo) 


HORIZONTAL POLARIZATION 


Incidence 

Angle 

^a 

^bres 

■ — 1 
Oi 

cr 

^bres “ 

o 

0 

109.08 

105.87 

- 0.02 

“ 3.23 

o 

10 

108.81 

110.29 

0.94 

2.42 

20° 

104.80 

105.70 

0.64 

1.54 

30° 

97.80 

98.19 

- 0.17 

0.23 

40° 

88.18 

87.47 

- 1 .08 

- 1.79 

o 

50 

80.84 

72.69 

2.74 

- 5.40 

o 

60 

74.02 

60.26 

9.32 

- 4.45 

O 

70 

97.44 

63.26 

47.49 

13.31 

80° 

140.57 

59.61 

102.60 

21 .64 


VERTICAL POLARIZATION 


Incidence 
Angle • 

Ta 

T 

bres 


^bres ~ "^h 

O 

0 

109.09 

105.88 

- 0.01 

- '3.22 

o 

10 

111.32 

112.72 

0.92 

2.32 

o 

20 

115.12 

115.85 

0.68 

1.40 

O 

30 

121.84 

121.86 

0.11 

0.14 

O 

40 

132.72 

131.49 

- 0.47 

- 1.70 

O 

50 

152.44 

146.93 

1.87 

- 3.63 

60° 

181.29 

177.35 

4.41 

0.47 

70° 

231..74 

236.06 

14.77 

19.10 

80° 

283.17 

290.88 

13.47 

21.17 
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TABLE XU II 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and One Restoration 


(Antenna = 8X horn, p = 26 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 

. 

T 

bres 

^a-'^b 

^bres " 

0° 

113.88 

107.11 

4.78 

- 1.99 

10° 

110.04 

108.22 

2.17 

0.36 

20° 

106.6-1 

103.86 

2.44 

- 0.30 

30° 

101.66 

97.04 

3.70 

- 0.92 

0 

40 

95.44 

87.88 

6.18 

• - 1.38 

0 

50 

95.99 

81.16 

17.89 • 

3.06 

o 

60 

95.64 

71.76 

30.93 

7.05 

0 

70 

135.46 

64.13 

85.52 

14.18 

o 

80 

180.96 

45.97 

142.99 

8.00 

VERT 

[CAL POLARIZAT 

ION 

Incidence 

Angle 

■Ta 

T 

bres 

1 

- T. 
a b 

- Tb 

O 

0 

113.94 

107.16 

4.84 

- 1.94 

10° 

112.58 

110.96 

2.19 

0.57 

o 

20 

116.84 

114.91 

2.39 

0.46. 

30° 

125.10 

122.56 

3.37 

0.84 

o 

40 

138.18 

135.18 

4.98 

1.98 

50° 

162.06 

158.00 

11.49 

7.44 

60“ 

191.71 

187.00 

14.83 

10.12 

o 

70 

243.91 

210.83 

26.95 

- 6.13 

80“ 

282.27 

202.52 

12.57 

- 67.19 
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TABLE XLIV 

Restored Antenna Temperatures for Finite Wave Tank with 
Random Error, Interpolation, and Three Restorations 

(Antenna = 8X horn, p = 26 feet, f = 10.69 GHz, T = 284°K, S = 0 °/oo) 

HORIZONTAL POLARIZATION 


Incidence 

Angle 


"'"bres 

Ta-Tb 

T - T 

' bres ' b 

0° 

113.88 

106.67 

4.78 . 

- 2.43 

10° ■ 

110.04 

106.82 

2.17 

- 1.04 

20° 

106.61 

101.23 

2.44 

- 2.93 

30° 

101.66 

93.00 

3.70 

- 4.97 

O 

40 

95.44 • 

82.45 

6.18 

- 6.81 

O 

• 50 

95.99 

77.42 

17.89 

- 0.67 

■ 60° 

95.64 

71.00 

30.93 ■ 

6.29 

0 

O 

135.46 

73.84 

85.52 

23.90 

CO 

o 

o 

180.96 

60.11 

142.99 

22.14 


VERTICAL POLARIZATION 


Incidence 

Angle 


"^bres 

T. - T. 
a b 

^bres “ ^b 

0° 

113.94 

106.72 

4.84 

- 2.38 

10° 

112.58 

109.45 

2.19 

- 0.95 

20° 

116.84 

111.74 

2.39 

- 2.71 

30° 

125.10 

117.59 

3.37 

- 4.13 

o 

0 

138.18 

128.54 

4.98 

4.65 

50° 

162.06 

154.77 

11.49 

4.21 

60° 

191.71 

189.37 

14.83 

12.49 

O 

O 

243.91 

238.35 

26.95 

21.39 

o 

o 

CO 

282.27 

263.17 

12.57 

- 6.54 
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Fig. 35. Restoration of the Finite Wave Tank Data 
with Random Error and Interpolation 
(Antenna = 12X horn» p = 13 feet, one 
iteration). 
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Fig. 36. Restoration of the Finite Wave Tank Data 
with Random Error and Interpolation 
(Antenna = 8X horn, p = 13 feet, one 
iteration). 
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Fig. 37. Restoration of the Finite Wave Tank Data 
with Random Error and Interpolation 
(Antenna = 12X horn, p = 26 feet, one 
iteration) . 
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Fig. 38. Restoration of the Finite Wave Tank Data 
with Random Error and Interpolation 
{Antenna = 8X horn, p *= 26 feet, one 
iteration). 
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there is enough difference between the antenna and brightness 
temperatures to prevent the error from dominating the process. 

Even at the incidence angles where the smoothed antenna tempera- 
tures are more accurate than the restored brightness temperatures, 
the restoration process does not yield results that indicate 
instability. When interpolation is used, there is some improvement 
in the horizontal polarization data but none in the vertical 
polarization with one restoration. Multiple iterations yield less 
accurate results. In addition to the random error, there is con- 
siderable interpolation error with the rapidly varying functions 
involved in the 12X horn and 26 foot boom case. 

For the 8X horn and the 26 foot boom, improved results are 
obtained by restoring the smoothed antenna temperatures with and 
without interpolation. With no interpolation, the results with 
three restorations are much superior than those with one restora- 
tion. The inversion is stable even with the added error and the 
results are improved significantly through the restoration process. 
With interpolation, the three-restoration results are inferior to 
those with one iteration. The results obtained with one restoration 
are, however, a definite improvement over the smoothed and inter- 
polated antenna temperatures. This is true for both polarizations 
and nearly all incidence angles. 

To partially summarize the parametric studies for the NASA 
finite wave tank, four tables will now be presented to show the 
recoimiended number of restorations for the various antenna, boom 



length, and data sampling combinations. Tables XLV, XLVI, 

XLVII, and XLVII summarize the 12X horn and 13 foot boom, 8X 
horn and 13 foot boom, 12X horn and 26 foot boom, and 8X ho»in 
and 26 foot boom cases. An "X", “V", or "H" indicate the 
recommended number of restorations for vertical -horizontal , 
vertical, or horizontal polarizations, respectively. Recommending 
no restorations indicates that the antenna temperature is a more 
accurate estimation of the true brightness temperature than the 
restored brightness temperature. These tables were based on the 

O O 

accuracy of the data from 0 to 60 incidence angle, since the 

O 

data above 60 is of little practical concern. With these tables, 
one should be able to use the most effective number- of restora- 
tions for the system under investigation. 

It should be noted that the two-restoration data- was investi- 
gated, but at no time did it yield the best results. 

Appendix II contains a listing of the Fortran program that 
performs the three-dimensional inversion. 

In addition to accounting for the non ideal pencil beam 
characteristics of the antennas, there is another major factor 
that needs to be compensated for in the restoration of measurements. 
This is the cross-polarization in the radiation characteristics of 
the antenna which was discussed in Section E of the theory. In 
practice, horns as well as other aperture antennas are not per- 
fectly polarized even in the principal planes, but do possess a 
smaller orthogonal component to the principal field. Since the cross- 
polarized term is orthogonal to the principal component, it responds 



TABLE XLV 


Optimum Restoration for the Finite Wave Tank 
with the 12X Horn Antenna and the 13 Foot Boom 


■’■ype of Data Sampling 

Recommended Number of Restorations 

0 

1 

2 


No Error 




X 

SPLINE Ihterpolation 


X 



Linear Interpolation 


X 



Random Error - No Interpolation 

X 




Random Error - Interpolation 

X 





cn 

"vi 


















TABLE XLVI 


Optimum Restoration for the Finite Wave Tank 
with the 8X Horn Antenna and the 13 Foot Boom 

Recommended Number of Restorations 


Type of Data Sampling 


o 


No Error 


SPLINE Interpolation 


Liiiear Interpolation 


Random Error - No Interpolation 


Random Error - Interpolation 



1 


X 


X 


CJ» 

03 


CO 


















TABLE XLVII 


Optimum Restoration for the Finite Wave Tank 
v/ith the 12X Horn Antenna and the 26 Foot Boom 


Type of Data Sampling 

Recommended Number of Restorations 

0 

1 

2 

3 

No Error 




X 

SPLINE Interpolation 




X 

linear Interpolation 


X 



Random Error - No Interpolation 




X 

Random Error - Interpolation 

V 

H 












TABLE XLVIII 


Ootimum Restoration for the Finite Wave Tank 
with the 8X Horn Antenna and the 26 Foot Boom 


Type of Data Sampling 

Recommended Number of Restorations 

0 

1 

c 

3 

No Error 




X 

SPLINE Interpolation 




X 

Linear Interpolation 


X 



Random Error - No Interpolation . 




X 

Random Error - Interpolation 


X 




o 


















to the polarization which is orthogonal to the principal wave for 
a given scan. As the incidence angle becomes larger and the 
difference between the horizontal and vertical emissions of the 
water get larger, the effect of this cross-polarization becomes 
more pronounced. To show the effect of cross-polarization on the 
measurements, antenna temperatures have been calculated for 
different assumed values of cross-polarization. The data listed 
in Table XLIX shows the effect of the different cross-polarizations 
for the 12X horn antenna and the 13 foot boom. Similar results are 
shown in Table L for the 8X corrugated horn and the 13 foot boom. 
Table LI shows the effect of cross-polarization for the 12X horn 
and the 26 foot boom and Table LII lists the results for the 8X 
horn and the 26 foot boom. From these tables, one can immediately 
conclude that cross-polarization becomes a very significant factor 
no matter how narrow the antenna pattern is or which boom is used. 
This is to be expected since the effect of cross -polarization is 
mainly a function of the difference between the orthogonal radiation 
characteristics of the environment. 

To show how well this cross-polarization phenomenon can be 
compensated for in the restoration process, antenna temperature 
profiles have been calculated, for various a's, assuming -20 dB 
cross-polarization and then restored to examine its importance. For 
the 12X antenna and the 13 foot boom, the results, with three resto- 
rations, are listed in Table LIII. By comparing these results with 
those in Table V, it becomes clear how well the restoration process 
compensates for the cross-polarization. The antenna temperatures are 



TABLE XLIX 


Antenna Temperatures for the Finite Wave Tank with Cross-Polarization 
(Antenna = 12X horn, p = 13 feet) 

(f = 10.69 GHz, T^ = 284 “k, S = 0 °/oo) 



Vh 

Tah 



"a» 

T 

av 

— 1 

< 

’’av ■ 


none 

-25 dB 

-20 dB 

-15 dB 

none 

-25 dB 

-20 dB 

-15 dB 

o 

a = 0 

109.15 

109.15 

109.15 

109.15 

109.15 

109.15 

109.15 

109.15 

a = 20° 

104.22 

104.25 

104.32 

104.54 

114,56 

114.53 

114.46 

114.24 

a = 40 

89.36 

89.50 

89.80 

90.71 

133.55 

133.41 

133.11 

132.20 

O 

a = 60 

65.69 

66.04 

66.80 

69.13 

177.97 

177.62 

176.86 

174.53 

o 

O 

CO 

II 

a 

69.54 

70.17 

71.50 

75.62 

267.99 

267.36 

266.03 

261.91 


























TABLE L 


Antenna Temperatures for the' Finite Wave Tank with Cross-Polarization 
{Antenna = 8A horn, p = 13 feet) 

(f = 10.69 GHz, T^ = 284°K, S = 0 °/oo) 




Tah 

Tah 

"^ah 

> 

1 — 

T 

av 

I'av 

^av 


none 

-25 dB 

-20 dB 

-15 dB 

none 

-25 dB 

-20 dB 

-15 dB 

o 

a = 0 

109.37 

109.37 

109.37 

109.37 

109.37 

109.37 

109.37 

109.37 


104.47 

104.50 

104.57 

104.79 

114.84 

114.81 

114.74 

114.52 

a = 40 

89.96 

90.10 

90.40 

91.32 

134.22 

134.08 

133.78 

132.86 

O 

a = 60 

70.48 

70.83 

71 .57 

73.84 

180.13 

179.78 

179.04 

176.77 

11 

00 

o 

o 


92.74 

93.9,0 

97.46 

263.84 

263.30 

262.14 

258.58 




















TABLE LI 


Antenna Temperatures for the Finite Wave Tank with Cross-Polarization (Antenna = 12A horn, p = 26 feet) 


(f = 10.69 GHz, T^ = 284°K, S = 0 °/oo) 




Tah 

— ] 


T 

av 

— 

< 

^av 


CROSS- 

POLARIZATION 

none - 

-25 dB 

-20 dB 

-15 dB 

none 

-25 dB 

-20 dB 

-15 dB 

O 

a = 0 

109.70 

109.70 

109.70 

109.70 

109.69 

109.69 

109.69 

109.69 

a = 20‘’ 

104.84 

r ' 

104.87 

104.94 

105.16 

115.13 

115.10 

115.03 

114.81 

a = 40 

90.68 

90.82 

91.12 

92.03 

134.52 

134.38 

134.08 

133.18 

a = 60 

73.50 

73.84 

74.57 

76.80 

181.23 

180.89 

180.16 

177.93 

a = 80 

129.16 

129.63 

130.63 

133.71 

277.43 

276.97 

275.97 

272.89 


a> 

45 » 



















TABLE LI I 


Antenna Temperatures for the Finite Wave Tank with Cross-Polarization (Antenna = 8Ahorn, p = 26 feet) 




(f 

= 10.69 GHz, T^ = 

284 K, S 

o 

o 

o 

o 

n 




^ah, 

Tah 

Tab ’ 

T u 
ah 

T 

av 

1 

■ T 

av 

Lv 

Lv 

CROSS - 
POLARIZATION 

none 

-25 dB 

-20 dB 

-15 dB 

none 

-25 dB 

-20 dB 

-15 dB 


112.92 

112.92 

112.92 

112.92 

112.91 

112.91 

112.91 

112.91 

HB 

108.41 

108.44 

108.51 

108.72 

118.48 

118.45 

118.38 

118.17 

a = 40 

97.85 

97.99 

98,27 

99.15 

140.03 

139.90 


138.74 

a = 60 

92.11 ■ 

92.42 

93.08 

95.11 

190.08 

189.77 

189.11 

187.08 

a = 80 

169,85 

170.19 

170.92 

173.17 

278.28 

277.94 

277.21 

274,-96 


cn 

cn 
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TABLE LIII 

Restored Antenna Temperatures for Finite Wave Tank with -20 dB 


Cross-Polarization and Three Restorations 


(Antenna = 12> 

horn, p = 13 

feet, f = 10.69 GHz, T = 

m 

284°K, S = 0 ° 


HORIZONTAL POLARIZATION 


Incidence 

Anqle 

Ta 

T 

bres 

T:, - T. 
a b 

^bres “ 

O 

0 

109.14 

109.08 

0.04 

- 0.02 

o 

10 . 

107.95 

107.88 

0.08 

0.01 

o 

20 

104.32 

104.17 

0.16 

0.01 

o 

30 

98.26 

97.98 

0.30 

0 n? 

0 

40 

89.80 

89.26 

0.54 

0.00 

O 

50 

79.16 

77.96 

1.06 

- 0.14 

O 

60 

66.80 

64.52 

2.09 

- 0.19 

O 

70 

59.74 

50.26 

9.80 

0.32 

0 

80 

• 71.50 

40.10 

33.53 

2.13 

VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

T 

bres 

Ta-Tb 

^bres " "^b 

0° 

109.18 

109.12 

0.08 

0.02 

10 ° 

110.44 

110.41 

0.04 

0.01 

0 

20 

114.46 

114.46 

0.02 

0.02 

o 

30 

121.69 

121.75 

- 0.03 

0.03 

o 

40 

133.11 

133.17 

- 0.08 

- 0.02 

o 

- 50 

150.52 

150.12 

- 0.05 

- 0.45 

o 

60 

176.86 

176.20 

- 0.02 

- 0.68 

o 

70 

217.48 

220.45 

0.52 

3.49 

CO 

o 

o 

266.02 

282.85 

- 3.68 

13.15 
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different for the two cases, but the restored brightness 
temperatures are almost identical. The effect of the. cross- 
polarization is therefore accurately compensated. In Table LIV, 
similar results are shown for the 8X antenna and 13 foot boom 
which can be compared with those shown in Table VII. Table LV 
shows the results of the restoration process of antenna temperatures 
with -20 dB cross-polarization for the 12X horn and the 26 foot 
boom. Comparing Table LV with Table IX, one can see that again 
the inversion process is able to remove the effect of the cross- 
polarization. The results of the inversion of the antenna tempera- 
tures with cross- polarization for the 8X horn and the 26 foot boom 
are shown in Table LVI. Comparing Table LVI with Table XI yields 
the same conclusion that the effect of the cross-polarization has 
been removed. It is concluded that the restoration process removes 
the effect of the cross-polarization for any antenna and boom 
length combination. 

Recently,' some preliminary measurements have been made on the 
wave tank system at NASA Langley Research Center, Hampton, 

Virginia. It would then be fruitful to examine the restoration of 
the data even though it may not be very accurate but it is 
representative of the response of the system. The measurements 
were taken at a frequency of 10.69 GHz using a 12X corrugated horn 

o o 

and the 26 foot boom. The values of a that were used were 0 , 10 , 

O O o o' 

20 , 30 , 40 , and 50 . For the first three values of a, the 
measurement had to be adjusted to remove a contribution attributed 
mainly to the standing wave pattern produced between the antenna 
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Incidence 


o 

0 • 


TABLE LIV 

Restored Antenna Temperatures for Finite Wave Tank with -20 dB 
Cross- Polarization and Three Restorations 

(Antenna = BX horn, p = 13 feet, f = 10.69 GHz, T^ = 284^, S = 0 °/oo) 

HORIZONTAL POLARIZATION 

nee T T - TT I 





a U 


109.34 

109.13 

0.24 

0.03 

108.17 

107.88 

0.30 

0.01 

104.57 

104.12 

0.41 

- 0.04 

. 98.64 

97.80 

0.68 

- 0.16 

90.39 

89.01 

1.13 

- 0.25 

81.03 

77.83 

2.93 

- 0.27 

71.60 

64.69 

6.89 

' - 0.02 

73.29 

51.66 

23.35 

1.72 

93,89 

39.27 

55.92 

1.30 


VERTICAL POLARIZATION 


— Angle- 

'a 

bres 

_ 

T - T 

bres b 

0° 

109.39 

109.17 

0.29 

0.07 

0 

10 

110.66 

110.39 ■ 

0.26 

- 0.01 

o 

20 

114.73 

114.37 

0.29 

- 0.08 

0 

30 

122.13 

121.39 

0.41 

- 0.33 

o 

40 

133.78 

132.58 

0.59 

- 0.61 

50 ° 

151,91 

150.04 

1.34 

- 0.53 

60 

■ 179.04 

177.90 

2.16 

1.02 

O 

70 

219.91 

229.94 

2.95 

12.98 

0 

o 

CO 

262.14 

272.08 

- 7.56 

2.38 





































































Restored Antenna 


• Cross - 

(Antenna = 12X horn, p = 
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TABLE LV 

Temperatures for Finite Wave Tank with -20 dB 
Polarization and Three Restorations- 

:6 feet, f = 10.69 GHz, T^^ = 284°K, S = 0 °/oo) 


HORIZONTAL POLARIZATION 



T 

. bres 

T;, - T. 
a b 

miQiiiii^mmiii 


109.04 

1.04 • 

- .06 


107.74 

.57 

- . *13 


103.88 

.78 

- .29 


mmm 

1.17 

- .41 


88.77 

1.86 

- .48 


77.74 

5.30 

- .35 


65.03 

9.86 

.32 


52.94 

39.51 

2.99 


41.10 

92.66 

3.13 

VERTICAL POLARIZATION 


^bres 




109.08 

1.08 



110.25 

.52 

- .15 


114.08 

.58 

- .37 


121.13 

.70 

.59 


132.45 

.89 

- -.74 


150.32 

2.43 

- .25 


179.33 

3.28 

2.45 - 

231.86 

9.50 

14.90 


179.33 

231.86 

280.98 


6.26 


11.27 
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TABLE LVI 

Restored Antenna Temperatures for Finite Wave Tank with - 20 dB 
Cross - Polarization and Three Restorations 


(Antenna = 8X 


Incidence 

Angle 


horn, p = 26 feet, f = 10.69 GHz, T^ = 284 K, S = 0 °/oo) 
HORIZONTAL POLARIZATION 




5.43 


111.38 


108.51 


104.13 


98.27 


96.18 


93.08 


125.92 


,92 


^bres 

Ta-^b 

^bres “ 

106.80 

6.33 

- 2.30 

1 

1 

3,52 

- .63 

103.58 

4,35 

- .59 

97.63 

6.16 

- .33 

89.39 

9.01 

.13 

80.24 

10.08 

2.14 

67.60 

28.37 

2.89 

54.05 

75.98 

4.10 

38.14 

132.95 

.17 


VERTICAL POLARIZATION 


Incidence 

Angle 

Ta 

^bres 

-Tb 

^bres ~ ^b 

O 

0 

115.45 

106.82 

6.35 

- 2.23 

o 

10 

113.82 

109.81 

3.42 

- .59 

o 

20 

118.38 

113.85 

3.93 

- .59 

o 

30 

126.75 

121.47 

5.03 

- .26 

0 

O 

139.61 

133,91 

6.42 

.72 

o 

50 

161.16 

156.82 

10.59 

6.25 

o 

60 


188.17 

12.23 

11.30 

0 

70 

237.73 

228.98 

20.76 

12.02 

O 

80 

277.21 

253.54 

7.50 

-16.16 

















































































and observed surface. The inverted data, using one restoration, 

o 

is shown in Figure 39. Since the 3-0 points are the more accurate 
-values on the restored brightness temperature profiles, curves were 
drawn through these points and are shown in the figure. For the 
larger incidence anglesi the restoration process shows a significant 
difference between the restored brightness temperatures and the 
measured antenna temperatures. 

B. Infinite Tank(Ocean) Data 

The developed programs can be used to predict and/or 
restore data from observations made at oceans or other large bodies 
of water. In these cases, the dimensions of the finite wave tank 
can be adjusted to fit the particular need. In Table LVII and 
LVIII are lists of data obtained by calculating the antenna tempera- 
ture profiles from the empirical brightness temperatures and then 
restoring them to recover the original brightness temperatures 
using the 12X and 8X horns, respectively. The restored 
brightness temperatures are almost exactly equal to the original 
profiles in these cases. The restoration process works better for 
the infinite tank case than in the finite tank case because the 
functions involved are smoother. All the finite tank cases involve 
a discontinuity in the water brightness temperature profile at the 
edge of the wave tank which causes high frequency content to be 
included in its spectrum. Since practical antennas can not detect 
this discontinuity (because of limited spectral resolution), the 
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Fig. 39. Measured Total Antenna Temperatures and Restored Water 
Brightness Temperatures for the NA5A laRC Wave Tank, 



TABLE LVII 


Restoration of Error-Free Infinite Tank 
Data (Antenna = 12X horn) 

(f = 10.69 GHz, T^ = ,284°K, S = 0 °/oo) 


B 

Tah 

Tbh 

T 

bresh 


Tbv 

T 

bresv 

0 

109.105 

109.099 

109.101 

109.105 

109.099 

109.097 

5 

108.712 

108.710 

108.701 

109.520 

109.508 

109.516 

10 

107.905 

107.905 

107.902 

110.380 

110,358 

110.362 

15 

106.145 

106.149 

106.149 

112.292 

112.249 

112.251 

20 

104.306 

104.316 

104.315 

114,338 

114.273 

114.274 

25 

101.164 

101.182 

101.182 

117.967 

117.862 

117.864 

30 

98.286 

98.3T1 

98.311 

121.442 

121.299 

121.299 

35 

■93.754 

93.788 

93.789 

127.237 

127.028 

127.029 

40 

89.836 

89.876 

89.879 

132.597 

132.324 

132.326 

45 

83.929 

83.977 

83.979 

141.369 

140.991 

140.990 

50 

77.269 

77.320 

77.320 

152.425 

151.909 

151,908 

55 

71.812 

71.861 

71.859 

162.580 

161.938 

161.940 

60 

64.009- 

64.037 

64.033 

179.234 

178.405 

178.407 

65 

57.876 

57.856 

57.867 

194.561 

193,613 

193.620 

70 

49.671 

49.485 

49.556 

219.304 

218.492 

218.451 

75 

44.174 

43.540 

43.839 

240.448 

240,648 

240.552 

80 

41.688 

37.871 

36.234 

263.427 

270,379 

270.476 




TABLE LVin 


Restoration of Error-Free Infinite 
Tank Data (Antenna = 8X horn) 

(f = 10.69 GHz, = 284°K, 5=0 °/oo) 


3 

T 

•ah 

Th 

^bresh 

Tv 

> 

1— 

T 

bresv 

0 

109.112 

109.099 

109.101 

109.112 

1 . I 1 

109.099 

109.099 

5 

108.719 

108.710 

108.701 

109.531 

109.508 

109.519 

10 

107.910 

107.905 

107.898 

110.401 

110.358 

110.367 

15 

106.145 

106.149 

106.146 

112.336 

112.249 

112.253 

20 

104.303 

104.316 

104.315 

114.408 

114.273 

114.276 

25 

101.154 

101.182 : 

101.184 

118.083 

117.862 

117.864 

30 

98.269 

98.311 

98.312 

121.603 

121.299 

121.300 

35 

93.728 

93.788 

93.787 

127.476 

127.028 

127.029 

40 

89.805 

89.876 

89.880 

132.910 

132.324' 

132.329 

.45 

83.890 

83.977 

83.975 

141.804 

140.991 

140.997 

50 

77.234 

77.320 

77.330 

153.013 

151 .909 

151.926 

55 

71.784 

71.861 

71.856 

163.296 

161.938 

161.949' 

60 

64.028 

64.037 

64.087 

180.098 

178.405 

178,370 

65 

57.981 

57.856 

57.951 

195.427 

193.613 

193.444 

70 

50.069 

49.485 

49.231 

219.419 

218.492 

217,892 

75 

45.496 

43.540 

42.715 

238.464 

240.648 

240.690 

80 

44.689 

37.871 

33.979 

252.473 

270.379 

272.900 
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exact brightness temperature profile can not be restored. With 
the smooth functions involved in the infinite tank case, one can 
obtain a very good approximation of the 

At Cape Cod Canal, Massachusetts, measurements were made by 
Swift [15] over a body of water which can be modeled as an 
infinite tank in one direction and finite in the other. The 
antenna used was a horn operating at 7.55 GHz and whose principal 
plane power patterns are shown in Figures 40 and 41. These two 
patterns were combined to construct the total three-dimensional 
pattern by using (122) and (123). In Figure 42 the measured ■ 
antenna temperatures have been plotted along with the restored 
(two- and three-dimensional) and the empirical brightness tempera- 
ture profiles. As can be seen, the restored three-dimensional 
and the empirical curves are very similar and different from the' 
measurements. The three-dimensional restoration works very well 
and it is more accurate than the two-dimensional, especially for 
larger incidence angles. 

As has been previously stated, a two-dimensional approximation 
of the wave tank system [7,8] has been used that takes advantage of 
the vector alignment in the ^ ~ 2 plane. In this plane 0 and h 

/S /N 

are aligned together as are 0 and v. This means that for the verti- 
cal scan only the vertical brightness temperature is received by 
the antenna and similarly the horizontal brightness temperature for 
the horizontal scan. However, the vector alignment in the other 
planes is not perfect and the opposite brightness temperature will 



RELATIVE POWER (dB) 


17B 



Fig. 40 . E-plane Power Pa 
Cape Cod Canal Ai 


: of the 7.55 GHz 
a. 



RELATIVE POWER (dB) 





TEMPERATURES { K) 


178 . 



Fig. 42. Measured Total Antenna Temperatures, Restored 
and' Empirical Water Brightness Temperatures for 
Cape Cod Canal Experiment. 



also contribute. The effect of this cross-coupling can only be 
accounted for by the three-dimensional modeling which calculates 
the vector alignment and integrates over all values of 9. 

The strength of the cross-coupling has been taken into account 
in the’ restoration process by the functions WF-|, WF^, WF^, and WF^ 

■in (144a)- (145b). In Figure 43 these functions have been plotted for 
the finite wave tank system and the 8X horn antenna and for the 
Cape Cod canal and the 7.55 GHz parabolic dish. From the data 
shown in these two figures, it is clear that cross-coupling is 
significant and the three-dimensional inversion is essential. 
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V. Conclusions . 

In the course of this investigation, the three-dimensional 
vector interaction between a microwave radiometer and a wave 
tank environment was modeled. With the computer programs developed 
one is able to predict the response of the radiometer to the known 
brightness temperature characteristics of the surroundings. More 
importantly, however, a computer program was developed that can 
invert (restore) the radiometer measurements. In other words, one 
can use this computer program to estimate the brightness temperatur 
profiles of a water surface from the radiometer response. 

The three-dimensional modeling of the problem was accomplished 
using two different coordinate system geometries. In one formu- 
lation, the z-axis was taken perpendicular to the radiometer 
antenna aperture and in the other formulation the x-axis was per- 
pendicular to the aperture. Computations were made to predict the 
radiometer response to the wave tank environment with both formu- 
lations and it was established that they both were' accurate models 
of the three-dimensional vector interaction. The three-dimensional 
models were also compared to a previously used two-dimensional 
scalar approximation of the problem. From this comparison, it was 
established that, unless the antenna used has a high main beam 
efficiency, the three-dimensional vector formulation is necessary 
to achieve an accurate result. 

With the x-axis formulation, it was shown that inversion 

(restoration), of the data was possible. Antenna temperature 

■\ 

profiles for the wave tank system were computed and brightness 
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temperatures were restored with a very good approximation. 

Errors were added to the computed antenna temperature profiles 
and the restoration process proved to be fairly stable. The 
effect of cross-polarization on the radiometer response was 
demonstrated as well as the capability of the restoration process 
to account for its presence. 

In addition to inverting (restoring) data for the wave 
tank system, it has been shown that the computer programs can be 
used to simulate. the viewing of large bodies of water. In this 
situation the restoration process is extremely accurate with the 
smooth functions involved. 

Preliminary measured data for the wave tank system, made 
available by NASA personnel, was restored taking into account the 
contributions from the surrounding earth and sky. Data taken at 
Cape Cod Can&l, Massachusetts, by NASA was also considered and 
resulted in a very successful restoration. 

With the restoration process and the future improved accuracy 
of the wave tank system, investigators should be able to experi- 
mentally verify the semi-empirical brightness temperature equations 
for various frequencies, salinities, incidence angles, and 
temperatures. The effect of surface roughness could then be 
experimentally measured with the controlled wave tank system. 

This knowledge could then be applied to analyze multiple frequency 
radiometer measurements received from satellites monitoring the 
ocean, in order to determine v/ind speed (surface roughness), water 
temperature, atmospheric conditions, and salinity. 
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Appendix I 

Transformation of Coordinates 


. A right-handed orthogonal coordinate system x, y, z can be 
transformed into any new right-handed orthogonal coordinate system 

II r II I ji t 

X , y ,2 , with the same origin, by three rotations about at 

least two different axes. An example of this- type of transformatior 
is shown in Figure I-l. this particular transform uses rotations 
about all three axes. The transformation of rectangular unit 
vectors for each rotation is described by the simultaneous equa- 
tions shown below in matrix form. For the rotation about the 


■x-axis {Figure I-l) 
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about. the y - axis (Figure 1-1} 
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and about the z-axis (Figure I-l) 
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-^11 -^ii r -<sii t 

The unit vectors x , y , and z can be found directly in terms 
of the original unit vectors (x,y,z) by combining (I*-l), (1-2), and 
(1-3) leading to 


/Nil ) 
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.•Ml 1 
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^11 1 
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(1-4) 


\t t U I 

Although 0 and 0 can be found directly from 6 and 0 by 

^11 i 

equating r to r , it will be more illustrative to show the trans- 
formation for each rotation. 

For the rotation about the x-axis we can solve (I-l) as 
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The radial vectors r and r are 

r = xcos0sin0 + ysin0sin0 + zcos0 

^ } I /V ) I I 

r'= XCOS0 sine + ysin0 sin9 + z cos0 
Using (1-5), (1-6) can be written as 


(1-5) 


(1-6) 

(1-7) 


/\ I 

r = X cos0sin0 

^ I 

+ y (cosAsin0sin0 + sinXcose) 

I 

+ z (-sinXsin0sin0 + cosAcos0) 
Equating (II-8) and (II-7) yields 

i r 

COS0 sine = cos0sinO 

I I 

sin0 sine = cosAsin0sin0 + sinAcose 


(1-8) 


(1-9) 

(I-IO) 
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cose -sinXsin0sin0 + cosXcose (L 

From (1-9), (I-IO), and (I-ll), we can find 9 and 0 
in terms of 0,0, and X as . 

0* tan"^ { *^osXsin0sin9 + sinXcose y (I- 

cos0sin0 

e. = cos {-sinXsin0sin6 + cosXcose} {I- 

For the rotation about the y-axis", (1-2) can be written as 


■^1 1 '-Nil 

y = [8J"- y 


cosy 0 siny x 
^1 

0 1 0 y 


-siny 0 cosy z 


We can express the radial vectors as 


'N|, II II ^11 „ „ H 

r = X cos0 sine + y sin0 sin0 + z cose 


/s I /VM 


r = X (cosycos0 sine - sinycose ) 


+ y.sin0 sine 


+ z (sinycos0 sine + cosycosB )' 


Equating (1-15) and (1-16) yields 


COS0 sine = cosycos0 sine - sinycose (] 

•I II II 

sin0 sine = sin0 sine {] 

>■ II I 

cose = sinycos0 sine + cosycosB (] 

From (1-17), (1-18), and (1-19), we can find the relationships 
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II li II 

between 0 and 0 and 6, 0» and u as 

o" = tan"'{ Aisg '^ iof' ^ ) (1-20) 

cosycos0 si no'- sinycos0 

II _ 1 II I 

6 = cos {sinycos0 s1n0 + cosycosO } (1-21) 

For the rotation about the z-axis, (1-3) can be rewritten as 
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The radial 

vectors are given by 
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All ^11 I 

If we set r = r , we find 


II 1 

COS0 

1 II 1 

sine 

II ii 11 II 

= COSVCOS0 sine + sinvsin0 sine 

(1-25) 
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1 II 1 
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sin0 

sine 

= -sinvcos0 sine + cosvsin0 sine 
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II 1 
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II 

= cose 

(1-27) 


If we multiply (1-25) by sinv, (1-26) by cosv, and then add we 
find that 

II I II I II ' II I 

sinvcos0 + cosvsinO = sin0 = sin(0 + v) 


(1-28) 
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from (1-27) and (1-28) we find 


U I II 


0 

= 6 

(1-29) 

II 1 

II 


0 

= 0 - V 

(1-30) 


As can -be seen, unless the rotation is only about the z-axis, 
the transformation of the spherical variables 9 and 0 involve 
transcendental equations. If the new coordinate system has a 
different z-axis, these transcendental equations cannot be avoided. 
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Restoration Computer Program 
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LATERAL EDGE OF THE VATPr-eartH BOUNDARY. 
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0BCR0S=THE AMOUNT DP C ROC S-PCLARl 2 ATT CN PRFSFNT IN CB . 
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« fJL?L9X«3wfl^0/ 

DATA NSTART/ 129, 12?, 12 7, 126, 126, 124, 126, 124 ,125 ,124,125,124,125, 

1 126*. 126’ ^^>125, 125, 125 ,125,126,125,127,125,127, 125,126, 

flDP = 2^ ’ ’ 

NOP21=NOP/2 .f 1.001 

NRPR=1 
IC0UNT=1 

(>EL = 360./FL0AT(NnP) 
or 6 ia=i,numalp 

ALP=ALPAI lAI 

’T**”V,TAAESH, TAAESV,TBPESH,T BRFSV ,NREST . 

1 ,NBRP,FFFF,LSTATF,HBCROS,PHa,W,TK,ALP,TC,INnATA) 

HPITF(6,7C0 ) alp 
DO 2 J= 100. 200 
A»G=t J-N0P21 l+OEL 

M( Jl .TERESHI J1 ,J,TAAHV I J) ,TAAESVIJ 1, 

1 T8P E SVI J ) ,ANG 
CONTINUE 

FORMAT (• ' ,3F10.3,I 10,4F10.3) 
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Ol'iliflil'813 ^'i|f|^^B^f”''Ct257),fM257 ),PYI256),WK(2S0 

THROUGHOUT ALL OF Tl**^ INTEG'^ATING SUBRCUTINES THE INTEGRATION WITI- RF^PFrr 
T? PEPFOPMEn with GAUSSIAN OUAORATupE numerical INTFGPAT^ 

WtlGHTWG FAC-CR SAPPLEO AND WX CONTAINS TFE 

5356 ,1.5067, 1,6707,1.A2B9f1-3827f 

? 'U?i9*"pp?7R ’9nnn')'’^•?knF’A »-T6341.. 67 lOA,. 57334,, A7385, 

nXT7*..5’rA^rt2Z?’ •?000^t, 12805, *07 035,. 029 02 ,*000549/ 
l”n^nAO?A2T?7^n 629792 9, 0. 02491 2439, 0.0326277B7, 

7 n"Si?oncJ >066,0. D4780S458 ,0.049598 040,0.049598049, 

O 039164137 rO. 032627787,0. 024912439^ 

3 0 .016297929,0. 0071 0843"^ ,0. 0142 1 69go ,0.032595869.0 049824879 

6 0:04ls24i79;0.i3259586?:0:9lla6Ut;°-°^ 

I F( N*>8R’.h E . I > GO TO ^17 

ji8"t-:tTf 

SriB^oTT?N^I."i'ANrfNV\\E'’^[.Sya='2s 

•MIll=BASF 2 LOG OF NQP 7PAU-6. 

. M(7)>=0 

«' 3) = 0 

frOP = ?56 
NDT=32 
M( 1 1=8 
M(21=0 
MI 31 = 0 
ALPlL=-5.0 

DBCP0S=10.>*7IPFCR0S/1 0. ) 

PFlNCF = l./{ l.+PBCPPS) 

CFPSSF=DPCROS/( l.+DBC PIS) 

DEGC0V=9C. 

NDP2=NDP/? + .001 
, NDP4=NDP/A + .001 
NOP21=NDP2-f I 
NDP41=NrP4i 1 
nn 8"i=i,NDr 
.TH^DISI I )=THFXXX( I ) 

WF( n = KX( I ! 
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IV G LfcVCL 


n 


INVFPT 


OATC = 7612? 


19/33/On 


P 1=3. 1 

H CCMTINUF 

L'tL=?.*pi/rLnAT(Nrpi 

• ■ Pfto=Pi/ieo. 

OELTHE = OEGCCV*r AO/PLOAT(N 
.... Rf-AC(5,lCri (RANCnPU),J = 

100 format (PF10.41 
717 Cf’NTINUF 

ALP = ALPO'=°AD 

PHI l = ATAN2IH/2, ^OHC;*SI N( ALP) ■ RHO*COS ( ALPU-ALP 
PMI2=ATAN2t W/2. -PHn-51 N( ALP) , RHO*CGS I ALP ) ) 4-ALP 

SCANNING^PLANE^^'’^" WATFP-FARTh BaUNDAPIES IN THE 

NPHI?=PHI2/0EI 

NPHl l=PHl l/DEL 

DFL THF=0E GCO V«R AO /FLOAT! NDT) 

on 6 1 = 1, NOT 


NOFAO=.2»(ALPn-ALPlL)4-. 00) 

«876 FORMAT(8FI0.3) 

PFAO(5,OP76) i TAPPVI )),J=1,\DP) 
R FAO( 5, 9876) ( TAPPHI J) ,J = l ,'IOP) 
IF IlMDATA.EO.l) GO ro. 1233 
IF f If DATA.FO.?) GO rr 1301 
GO TO 1302 

1301 continue 

1 = 1 

on 319 J=1,N0P,A 
TINTPVI I )=TAPPV( J) 

TINTRHI I ) = TflppH( j) 

I = I -4 1 

319 CONTINUE 

TINTRVI I )=TAPPV( 1) 

T INTPHI I ) = TAPPH( 1) 

no 110 J=l,65 

110 - 

ANGLF=0EL*(J-1) 

CALL TFPPA (ANGLE, YYl 
TAPPV(J1=YY 
in CONTINUE 

CALL SPLINA IAMG,T!N.'PH,65) 

0<) 112 J = 1,N3P 
AMGLE=OFL»( J-! ) 

CALL ’’ERPA t angle, YY) 

TAPPH(.JI = YY 
112 rOMTlMUF 
GO TO 1??T 

1302 CONTINUE 

IP ( 1NOATA.FO.3I Gn -^C 1303 
GO TO 1304 

1303 CONTIMUF 

1 = 1 

on 9 J=1,N0P,4 
TtNTRV( I 1 = TAPPV( J) 

TINTFHt I ) = tapph(.J) 

1 = I 4 1 

9 CONTINUF 

TINTDVd 1=TAPPV(1) 

TINTPHt I | = tapp.h( l) 
on 14 J=l,64 
I I = (J-1 )w4.vl.001 
TApPVI II )=TINTF V( J) 
TAPPH(II)=TINTPH(J) 

11=11+1 
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IV G LEVEL ■ ?l 


JKVER T 


DftTE = 7G122 


59/33/78 


TAPPV( in = .7 8«’TI\'TRV( J) + .25’’TINTPV(J*-I1 
TAPPHt ni = .7E<-TINTKH( J) + .25^TINTPH( J+1) 

I I = 1 1 + * 

TAPPvni 5=.50*TINTP VIJ)4-.50»TtMTRV( Jf i) . 

TAPPHC^I I )-.50»TINTPHI Jl4-. 50*TINTRHlJ+n 

TAPPV( in = .25*TIN'TkV( J) + .75 'TINTFVCJUl 
TAPPH( in=:.2'5»-TiMTPH( J)*-. 79 fTIVTPHU*-!! 

14 COMTIr.iif- - 
GO Tn 1323 
1304 CnNTI'lUE 

IF ( riDATA.ro. 4) GO TO 130? 

GO TO 1306 
130 5 Cr.NTT>'iOE 

SSSSS=278. 
on 10 J = 1',N0P ■ 

ANO( JlsDFl J-ll 
-DYIJI-1 . 

10 CONT.IVur 

C Al L ir sr.c- U ^ U NO , T/ p p V ,D 7 , S 3 S SS , NOP , A , R , c , D , WK 1 
tappvijT=a'( j + n 

11 COM TIN OF 

call ICS?GU(ANC, TAPPH,DY,SS3SS.NDP,A,B,C ,0,WK) 

DO 12 J=1,N0P 
TAPPHt J)=A(J + 1) 

12 CONTINUE 
on TO 133’ 

1306 CONTINUE 

IF t IMDATA .FO. 51 GO 70 1307 
CO Tfi 133^ 

1307 CONTINUE 
SSS':5 = ?7R. 

I- 1 

no 20U J=1,N0P,4 
TINTPVtn=TAPPV(J) 

T-INTPIIt I )=tapphI ,1) 

■ r=i + i 

209 CONTINUF 

T INT"V( I l = TAPPVt 1 ) 

TINTRHn)=TAPPHIl) 

1)0 210 J=1,NDP 
AMGt J) = 4.*-0EL’r( J-I 1 
OYt Jl=l. 

210 CIINTIniiF 

call irSSGII (ANG.TItlTPV.OY, ;SSSSt65 ,A,B,C.O,WK) 
or 2111-1.64 
on 211 11=1,4 
tU = 4«( I-IHII 
M=t II-l 1*-04L 

... . JAPPVt ! M ) = MD( ! ♦!! *h+C( I + I • 1 *m-«( 1 + 1) l»w+At I + l) 

211 CONTINUF 

’ <''NG.TIN''PH,DY,'.SSSS,65,A,B,C,0,WK) 

01' 212 1-1,64 
DO 212 n = l,4 
1 I ! = 4i«t I-l ) + I! 

H=( IT-U-OEl 

TAPPHt I II ) = l (111 I + l l*H+Ct I + ll ) '‘H + BII+D 1 ♦H+A ( I +1 ) 

,212 COMTIMUE 

1 333 CONTINUE 

IFILSTA TK.eo .5) 

1 CALL.NEUSKY (NPHIl ,'!PHI2 ,T',ESV,TAESH,TK,TF,bH0,FF( , 

IF ( L STATE.E0.6 1 

.T.ESV.TAESH ,TR BV.TBPH , LST AT E, FFFF, RHQ) 

IP (LSTA+E.FO.T) 

'T CALL THEAPT ( TRB V, TOPIt , PEP' I P , PE RHC , Pf RV P, PEPVC 1 
IF I LSTATE.ro. 7) 

I CALI EARSKY (NPHIl ,'lPHI2,T\ESV,TA£SH,TReV.TRRH,LSTATC,FFFF,RH01 
: TfUS loop SUP TRACTS ti'E ciijTii AND sky ccntributi cns fpcm tff nfasurements. 


DO 134 J+l.NDP 
TAPP VU 1 = TAPP VI Jl-TAt SVI Jl 
T.''PPM( J ) = TAPPH( J)-TAFSm J I 
CAMVt J1=TAPPV( J)' 

TAMHI J l^TAPPHtO) 
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iv 0 LEVPL ?1 


TNV=PT 


PATe = 7A122 


19/33/08 


13^ CONTINUr 

on 109 1!=1,MPFST 

CALL WATER < WHIl t NPH I •» , tai, y, , l STA TE , FFFP, OHC 1 

• no 7 j=j 

, TAV( J )=TAWVi J) 

TAIII J )=TAWM[ J) 

EFP V( J,l = TAMV( J )-TAV( i) 

£A“H( JI=7AHH( J)-TAH( J) 

TAPPVl J) = TAPPV( J)+FR>- V( Ji 
TAPPH( J l = T4pPH(J )+cpPH( J) 

TAPPh( J I =TAFPHt J I tfPPH ( J ) *p7p HP ( J) 4-FpRV ( lISDCDurf tl 
TAPPV( J1-TAPPV( JH-FRP V Up^RVP J ippPH I *pfepvr f‘i 

TAP (;SV( Jl-TAPPVI Jl 
7 CONTiFiUE 
109 CONTIMUE 

ALP.IL’ALPD' 

' NRPPsjMPRP 4- V 
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^ooroonononon 


il; G LEVFL 21 


PATE = 7hl22 


SUwnyrWIE- PEECEN ( NPHll,NPHI2,FFFF,pPEHP,p£RHCt PPnVP.PFPVC, 

ILSTt*T»frm) 

alpha and beta will he COHPnSED OF 


SUBPOUTINF. , ^ 

PEPHP=FPACTf CNAL f, CKTR [ RijTI PN CF TEH TO TAWH 
PFOMC^FRACTinNAL CONTRIBUTIDN OF TBV TO TAWH 
oCBVP-FlUCTinMAL (.CNTP 10‘JTI PN OF TBV TO TAWV 
pcRvC^FoACTj OMA L C. CNTP I BUT! ON OF TPH TC TAWV 

l‘''':?pT?^ (cVTs^F:?H?i;f IVu^ . ^'>T . ndp .del . dfltkf. alp 

CO'IMON/BG/GTHEH.CPHIH.GTHEV.G mi V .GCRCH , GCROV 
COMPLFX T1V(,PC(256) ,TlHCRQt256) 

,STlVr.(256) .CMPLX.CP 

i?S U’’'£*?561 .T1HC1256) ,T1VP(256) ,T1VC(256) 

CCHPLEX GTHEHI 256) ,GPHIH(256) ,GTHFV(2S6) jGPHIVI 256) 

CPHPLFX GCBOH<256),GrpOV(256) 

DIMENSION PERHPI256) ,PERHC( 256) ,PE«VP (256) .PFFVC 
DIHENSMTN T1HI 256) ,T1V(25A) 

DIMENSION INV(128),S( IZ8),M(3) 

DIMENSION TAPPVt256),.TAPPH(>56) ,STHE(65) 
dimension THFDIS(32) ,HF( 32) 

PI=^.1A159 

HAPyPI/IBO-. 

NALP=ALP/DFL+.5 

NOP2*Nnp/?+.0Pl 

NDPAyNDP/Af.001 

NDP21-NPP2+1 

NDPAI=NOPA-Fl 

N1=MD»21»NPHI1 

N2-MPP21-NPH12 

CP=CMPLX( 1,,0. ) 

CALP=COS(ALP) 
or 3 J = 1,,MPP 
riH( j )=i,o 
Tl V( J )=1.G 

stihp(j)=cmplx< 0. ,0. ) 

5T1HC,( J)=CMPLX( 0. .0. ) 

ST1VP(J)=C.MPLX(C.,C.) 

STIVCI J)=CMPLX< 0. ,0.) 

CDNTINIIF 

OP 10 1=1. MOT 

RFAOm (GTHEHI J) .jYl.NDP) 

REM)!)) (GPHIHI J 1,J = ) ,N0P> 

PFAD(l) (GTHFVI J),J=l,NOP) 

°FA0(1) IG)‘HIV( Jl,J=l,NDP) 

READII) (GCPOHIJ 1,J = 1 ,NDP) 

PEAO(l) IGCROVI J) ,J=1 ,N0P) 

OP 13 J=1.NDP 
T 1HP( J )=CMPLX( C.,C. 1 
TIHC(J)=rMPl X(C.,C.) 

T1 VP( J)=f MPLX( 0..0, ) 

T!VC( J)=CMPLX(0.,C.) 

T1VCPO(J)=CMPLX(0.,0. I 
TIHCPOI J)=CKPLX< 0, ,0. ) 

)3 CONTINUE 

00 A J=N2,II1 

PH1*FLPAT( IABSI J-LDP21 + NALP) ) «DEL+.0001 
THF = P1-THE0IS( I ) 

TPHI=t4n(phi ) 

THEFW=ATAM( FFFF /( PhO*C ALP*S ORT (1 . pni*TPHn ) > 

THI;Ew= THEp W+P I /2- 
STF1=SIN(thF) 

CTHF^cns( the ) 

SPMI=SI N( PHI ) 

FN(3MySPHI<-spHI/ISTH1*STH1 '»SPHI*SpMIFCTHE*CTHE) 

GAmMA=APSIN( ( SINtPHn*SIN(PIII I*SI NITHE) 4-CnS (THE) ‘ 

1 /SOB T{ SIN( PHI )*SIN( PHI ) + C’’TAKt THE ) *CCTAN (THE ) ) 1 


aEPROBUClBIlbl^^ W Tip' 
original^ IS' POOR 
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IV G LrVtL J1 


PERCSN 


DATE = 76122 


] 2 
10 


100 

200 


ll 


NGAM.MA=GAMMA /DE L+ . ^ 

NCAwrA=NOP?l-NALP + l'GAMMA 
lFnt«E .CT .THEFW) GL' TP 5 
T1HP( Jl = Tin(NQAt»H^ )*F6’CW+CP 
T 1HC( J)=T1H<MGA«HA)»( l.-FMGV) *CR 
T1VP( J>=TlV(NGArtMA )*FW0W*Co 
TIVC( J I’sTlVC NCAPHA I * ( 1 . -F NC»J 1 *CK 

T ivraot J J-2. + T1 VI MGAMMA) *S0P T! FNflls* 1 1 . -F NOW 1 )*CP 
T niCPni J ) = ?.*TlH(Nr,AMHA)«SOPT(FNnH'‘( l.-FNOW) I*CP 
CONTINUE 
CONTINUF 

CALL HAPP mKP,P,INV,'?,l,I‘=EPPl 
CALL HAPP ( TIHC ,H,I NV,S,l ,I "EPPl 
CALL PARK ( T1VP,M,INV,S,1 ,t=ERRl 
call HAPP ITl VC.HrINV.S, I , I=6PP1 
CALL HAPP( TlHCP0tHiINV,S,2,'FEPR) 

CALL HAPPITl VCBn,P,INV,St2 .IFERR) 

[!0 12 .1 = ] ,NOP 

$T1HP(J) = ST1HP( Jl + t T’.hpI J1 AITHEKU) +T1 HC [ J )«GPH IHIJ 1 
1 -tihCROIJ I»'GCP0H( J) 1 
I »STHEI 1 1 

‘ ST1HC(J)=STIHC(J)f(T 1VP( J1*.'.PHIH(J) ♦TlVCtJ)*GTMF.MJ) 
1 ♦Tivcpnu)»ccPOH( ji I 
1 ♦^Tuei 1 1 

5T1 VP( J) = ST1 VP( ,() + ( TlVPIJl+GPHIVI J) *T1 VC tJ ) «GT PEV ( J ) 
1 -TIVCPOIJ )»GCP0VIJ1 ) 

1 ■*";the( 1 1 

snvr ij)=sTivc( Ji*-( tip pi j)*gthevij) ftihcij)*gphiv(j) 

1 ♦TlHrPCI(J)*GCRr»V(J) 1 
1 «STHE( I 1 
CONTINUE 
CONTINUF 

acting 1 

CALL HAFP I STlH0,M,INV,S,-2 ,IFEPO) 

CALL HARM f STIPC tP.INV.S 1-2 ,IFEPP) 

CALL HAPP ( STlVP,M,IIIV,Si-2 tIFF.RPI 
CALL HATH (STlVCtM, 1NV,S,-2,IFERP) 

WK ITEI 6,1001 

FOPPATI • >, >PFPHP,FEPHC , PEP VP ,PEPVC ,I NOEX ' 1 
FOPMATC ',AF13.6,16) 

I'C 11 J = I,NOP 

DcNPOP^REALI STIHRUKSTIHCI JJ 1 
UENVFR = PFAL1ST]VPU)+ST1VCI J1 I 
. PERHPI JJ*|PEALISTll>PIJl) 1 VDENHOR 
PEPFCIJ )=1PFALI STII-CI J1 ) ) /OENHOR 
P1RVPIJ)«IPEALIST1VP1 J)) I/PENVER 
PFPVC IJ)=(RE-AL( STl VCI Jl) l/O^NVER 
PX1=PFRHPIJ1 
PX2=P60HC(J) 

PX’ = P=P VP( J) 

PXA=PFRVCI J) 

■ Pf-HPU l = PXl*P»INCi +PX2*CPn5SF 
PFRHC! M=PX2*PPINCFH>X1*CR05SF 
PERVPI J ) = PX3*pPINCf •■PXA«CRO?SF 
PEPVC IJ I = PXA*po INCr »-PX3*CRnSSF 
CONTINUE 
I 0 6 J = l,26e 

WRITE (6, 2001 PERHPU) ,PEFUr( J1 , PEP VP ( J) , PERVC ( J ) , J 

CONTINUE 

FETUPM 

FMD 


10/33/08 
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V C LtVfcL 2! 


M 4 TEP 


OATf = 7612? 


19/33/09 


C 

C 

X 

C 

C 

C 

C 


SiJIiRnuriM- WATf^(N('la 1 ,NPMI ; ,takV,T&WH,L5TATE, FFFF,DH01 

THIS SUElFnUTlN!; .CALCUl.&TFS THE 4NTENMA TEMPERATURE CCMTP I HUT [ON FROM THE 
WATER AS A FUFCTIUN OF BETA '-nP EACH PCLARIZATION USING AN ESTIMATE OF 
THF hfightness TEMPEPATUPE. 

TAl>V*VrP TICAL KA TER ANT-NKA TI=MPERATURE CONTRIBUTION 
'^AWH=HNP UONTAL water ANTENNA TFHp£PATl)RE CONT R I BUT 1 ON ‘ ' 


coMMC'/Bi /TAPP v.tapph.sthe t ndt.mcp.oel.ofuthe, alp 

1 ,PF inch »r.PrSSF,THFO[ S »KF 

OMMON/BLOCl /TPHP ,TH1IC , TRVP , TBVC ,$T AWH , STAHV ,T1 VCPC.T IHCRO 
C0MM0N/DLnC2/S,M, INV 

COMMON /B6/0THEH ,GPH! il ,GTHE V ,G FH I V ,TXRCH , GCPCV 
CDPOl F)f TlVCPnl 266) , TIHCPOI 296) 

COMPLEX GTHFHI 256) ,GPHIH(2S-5) .GTHEVI256) .GPH!V(256) 

COMP) EX GCPnH(?56),GrPnv(25S) 

complex TGHP(256I ,TBHC(256) ,THVP (256 1 jTBVC (256 ) 

COMPLEX STAWHI 256) ,$TAWV( 25b) ,CMPLX,CP 
niMENSION TBWH(256) .TBWVtaS-i) ,TAWV(2 56) ,TAWH(?56) 

DIMENSION IMV< 128) , St 128) ,M(3) 

DIMENSION TaPPV( 256) , TAPPHt 256) ,STHE( 65) 

DIMENSION THEDlSnai fWF(32) 

P 1-3. lAl?" 

RA0=PI/1PC. 

N6LP = ALP/D(-L+.5 
NL'P?=NOP/? + .COI 
NOPA=NDP/A+. CCl 
NOP21=NOP2+1 
NrPAl=NDPA+l 
' Nl=Nr)P?lFNPHIl 
N2-M01>21-MPHI? 

CR^EMPl X( l.,C.) 
rALP=cns(Ai p ) 

DO 3 J=1,NDE 
. TPWM( J)=TAPPMt Jl 
TCWVt J)=TAPPV(J) 

I Ft ( J .LT.NDPZl -llPhU ;l -OR. ( J.GT, N0P21 +)IPH[ I ) ) GOTO I 
GO ro , 2 
1 TPKHJ)=0.0 
TPWVtJ)=0.0 
7 COMTINUF 

. STAHIHJl sCMPlxIU. ,0. ) 

•STAWVIJ) =fMPLX(0. ,0. ) 

3 CCNTPIUF 

on 10 1 = 1 , NOT 

READin (GTHEHt J) ,J = ’ ,NDP) 

RtAOtl) (GPHlKt J) ,J=1 ,NDP) 

RFADtl) [GTMEVt J) , J=l,KDP) 

READd) (GPHI V( J) ,J=l ,NOR) 

RFAOt 1 )• ( GCROHt J) ,J = 1 ,NDP) 

KFAntll (CCPOVt J) ,J = ’.,NDP) 

UP 1? J=1,N0P 

T8HP( J) =CMPLXtu. , 0 . I , 

TPHC(J) =CMPLXtO. ,0- ) 

TBVPl J) =CMPLXIO.,P. ) 

IRVCIJ) =CHPLX(0. ,0. ) 

TlVCFOt J)=C«PLX( J. ,0. ) 

TlHCBOt J)=CMPlX(u,,0.) 

13 COMTINUF 

on A J=tl2,Nl 
TME=P !-THEOISt ! ) 

PH!=FLOAT( lABSt J-W/R > H-NA LP) ) *rjFLF. 0001 
TPH1 = TAN(PHI > 

rMEFW=ATANt FTFF/tP hO^CALP*S!)RT( l.+TPHI»TPHI ) 1 ) 

ThFTWs rHFFU<-PI /2. 

STH1=SIN( THF ) 

CTHF=:nStTHF> 

SPHI^SIMIPHl ) 

rNill, = SP)II*SPIII /f STIil »ST|-il ‘S’HdSPMl <-C THE ♦CTHF) 
T7D|tI=TAN)PHI lATANtPHl 1 
CT2THF=1./(TANI THE)*rANtTHF) ) 

GAMHA = ATAM SOP T< T2PH!+CT2TltE4-CT2THF*T2PH! I 1 

GAMMA=GAMMA/OFl 

NGaMMa-gAhma 

oelgam=gamma-ngahha 

NGAHMA=NDP21-NALP+NGAMHA 
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iV G LEVei 


n/33/OF 


T NCAMM/s )+tt LGA TP WH ( KG AMHA H ) -TBHH ( NGAFHA ) I 

IFiKGA'-iNiA.GE.ra ) TUV=TBWV(NH 
IF(THE.GT.THGEW) THH=0.0 
IFIThF.GT. THEEVl TPV=D.O 
TP HP ( J J= TPH*FHnw*CP 
T(ihC< J)=TbH»( l.-FKUH)»rp 

T 0 VP ( J ) = Ti) V*FMnw»G F 

TRvc( j)=Tov*n.-rNOMi tcp 

llHC«Ot JI=?.»TPH»S(JKT(F('!CH*( 1 .-FNCW) I *CR 

TI VCPFK JI=?,*TPV*it'RT<FNnW*tl,-FKCWI l+CR 

CON T I NUE 

COHTIMUF 

CAUI HftPit (TBHP,M,iNV,S.l ,I -FRR) 

CAM. HAFM ( TBhC.H,! NV,S,?,rFERR) 

CAUL H40M (TBVP.M.JNV.SjEtIFEPPI 
CALL HAFM (TbvC.M,: KV.S.2 ,I <=ERP) 

CALL HAPrU Tl VCROjH, INV,S ,2 , [FEPRl 
CALL HAPMtTlHCF n,K. I NV, S 1 2 . I F EPP 1 
on 12 J=l,Nfjp- 

STAWMC J) = STAWH( J) + ( THHP( J) >^TTHFH( J1 +TBHC ( J) *6 PH IH IJ I +T6V C IJ I 
J VsTHFl'n TT1VCR0(J )*GCPnHtJ ) ) 

, J > < ■'2 VC( J1 *3the vt J1 + TB VP ( J1 *CPH IV ( J ) ftbhp U ( 

I *STHn < J)»GPHI V( J) fTIHCROC J] *GCPnV { J)-T1VCR0(J 1«GCRnV(J 1 I 

COtlTINUF 
CO'ITINoe 
PEwIfO 1 

CAl L HAPK 1 STAWH,H,IKV,Si-2 ,1FF0R) 
call harp 1$TAWV,K, INV,S,-2,IFERP1 
no 11 j=i,NOP 
TAWVI J )=PFAL( STAWVI JJ) 

TAWIH J)=PEAL(STAWH( jn 
TX1=TAHV( J) 

TX2=rAWH( J) 

TAHVl J l=TXl*PP.INCF+TX2*CROSSF 
TAKH{ J)=TX?”PP INCH- T,<1*C PGR 3F 
comtihut 
FETUPM 

ENG 
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ononnoonrs 


IV r. LfVFL 


-rVl-L NfWSKY UATF.s 761?? n/33/08 

SimPTIiriNF NFWSKY ( N-’MII ,NP-II? AFSV .TAESH, TK, T F , F Hfl, FFFF ) 

' CALCUI./ TFS T'JF AN'T^mma TEMPERATURE CUNTR 1 RUT IPNS FROM THF 

taesv=tmf vertical antenva tempfrature contribution from the farth 

TAi:SH=THl- HORIZuNtAL ant-nNA TF^PEPATUPF CCNTRIPUTION FROM THE EARTH- 

C'iMMONZRE/CTHFH.OPHIM.GTHEV ,G TH I V ,GC RCH .OCRCV 

1 .GPt'IH t2SS) tGTHEVI256l ,GPHIV(256) 

COHPLFX GCROHI 256) .GCPOVI 2E'j) ,TBEXT(256 

i'^AAAiTa ^?|5<256),LMPIX,CR, 161256), StIesh(25GI,STA = SV(256), 

DI'trNSICN'’lNV( 12n),sn28),M(3) 

DIMENSION TAESVI256) ,TAESM( ’56) 

•^*^”V(256) .TAPPHI 156) ,STHE(65) 

OIMFNSION^^IFPISOi) ,WF(32) 

NAL'’^ALP/DEL + -5 
CALP=CnSlAL'P) 

N0pA=NDP/A.4-.001 

NDPHA=NOP-N0PA+l 

- N0P21=MriP/?»l.C0i 
NDPAl =NDP/*, + lA oor 
N l-^NOP21+NPHI 1 
N2=N0P21-NPHI2 
CK=rMPLX( l.,0.) 

■_ DO T J=1,NDP 

TnU) = CMPLX(TF,0.O) 

TRESUKMPLXIP.,0. ) 

STAESVI J )=CMPLXI 0. , C. ) 

STA'^SHI J I =CK.PL X( u. , ”. ) 

T rUNTINUF 

’’G IS THT GAUTCTIC !>P)CHTNESS temperature. 

TG=’>. 

TtFF=l.l?»TK-50. 

°Tn=J./TEFF 
TAI|P=-ALCG( 1 .-PTC) 

DU A J=l,Nt'PAl 
ETI = nS[DtL'^l J-U ) 

IF ( C Tt ,L T. O.OOOOOl ) CC Tr ?? 

I' ( J.Ea .NDPAl ) Gt. TO 22 
GO TO 33 

22 TIM J)=TeFF"CP 
GO TO A A 
33 COMTINDE 

SrCTHP^l ./f TI 

AA CON^INUf'"'^''** ' •-t-XP<-TAUO*SrCTHP) ) + TG *FXP( -T AUO*SFCTHP) )«C« 

II ( J.EO.) ) GO TC 55 
tu(nD->+ 2- J ) = Tp( j) 

55 CflNTlAilJF 
A COfITINUF 

ir (MALP.EO.O) GC TO 21 
on 13 .l*),NPP 
TG^XTI J) = T6( J)' 

13 CO'ITINlIf 

AAt.PUNAlPF) • 

'Of lA J -MAI PI, NOP 
-7|H J“NALP)=TBt XTIJ) 
lA CONTI nuf 

DO 15 J = 1,NALP 
TPINOP-NALPF.))=TliEXT( J) 

15 CONTInuF 

21 continue 

or 10 1=1, NOT - 

- THE = Pl-THFfI-Sn 1 
‘’fao(i) (gthfhi ji ,j=i ,Nnp) 
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IV G LtVEL zr 


MFWSKY 


DATE = 7M?2 


19 / 33/08 


12 


10 


II 


kFADIll tGr>Mrn(J),J = l,NDP) 

PFAD(I) (GTIIf V( J) , J = 1 ,NDP) 

PPAD(l) (GPHIVIJ),J = 1,N0P) 

PlAOn) (GCROHI J) ,J=l,NDP) 

“FAOdl- (CTRPVC J1 ,J=i,NDPl 
rn> 12 J=1,N0P 

Wil 7i J.C ? ^ 1 + LP 1 1 *DF L . 000 1 
I i*f-i i ~ TAN C P*-l I J 

THEFvliTMFPH+PI/2* ^HO”CALP »SOP Til .+TPHI>‘TPhI ) ) 1 
TB'^Sl J)=T5(J) 

CONTINUE 

CALL HAPP( Tpf C,M, imv,£,2 ,IF PP) 

L*o |o j=ljl*OP 

I'.M J> »1GTHEH( J1 +GFHIH J 1 I 1*STHF( tJ 

CDMTInuf SViJl H TP.ESI J) ♦IGTUFVl JI+GPHI VU) ) ) wSTHEC I 1 

PfklMO 1 

CALL HARP ( STAFSH.M, INVt S,- ’ ,IFEPP) 

^ ""7=^ ■ N V. 8 ; IFFPR ! 

IAF8H1 J l=RfcAL( STAfcSHlJ I J 
TACSVI Jl = PEAL ( STAGS VI J ) ) 

TX1 = taESH(J) 

TX2=TAFSV( J) 

TA^SMIJ ) = TXl<-PRti'iCF+ :x2*CPn ISF 
J/ '^Syl J 1 = 'X?*PPIiNCF+' Xl*CRO>SF 
CON T! 'HIT 
P FtuI'N 
FNI) 
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nnoooon 


IV G LEVEL 21 TPEAPT OATF = 7ftl2? 

suproutjnf tkfapk TPnv,TBPH ,pephp,pfbhc, PERVP,p=Rvr i 


1 9/33/Ofi 


10 


11 

12 


this subpputine ta\es the measured amenna temperatures ant rt-stopes them 

DONE TO OBTAIN A BETTER APPROXIMATION OF THE BRIGHTNESS 
TEMPEPATUPES OF THE -EARTH ANT SKY TO BE USED I N 'SUBPOUTINE EAPSKY. 

T8BV«THE PESTdKED VFPTISAL MEASUREMENTS 
T(1BM = THH FFSTURTI) HORIZONTAL PE ASUPEMF.NTS 

CriMHrj'l/BI /TAPP V.TAPPH.STHE, NOT, NOP, DEL.OEUTHE, ALP 
1 ,PP !NCF,CP0$SF,Tri6DIS,RF 
CPMMnN/JUOr.l /TBESH, tbeSV, STACSH ,STAESV,AAA1 ,FPRV,£PRH 
CO'(HDN/HIOC2/S,H, INV 

CnMHON/BA/GTHFH,G‘>H I H ,CTHE V ,G PHI V ,GCROH , GCPOV 

DIMFMSION TBBV(256) ,TI’BH( ?Sj) ,PFBHP{ 2S6I ,PERHCC>S6) ,PFRVP(756I 
1 ,1’FPVCI ?SA) 

COMPL'-X GTMFHI 256) ,GPHIH(25SI ,GTHFV(256I ,GPHIV(?56I 
C.'IMPIEX r-rPOM(256I,GCPnV(25'>) 

complex THESH(256) , CMPL X ,CR , TRES VI 256 1 , ST AESHI256 ) , S T AESV (256 ) , 

1 AAAI ( A ,256 I 

DIMENSION INVdZa) ,S<128I,M(5I 

TAESVt2561 ,TAESH( 256) ,ERRV(256) ,EPPH(256) 

DIMENSION TAPPVI2561 ,TAPPH(’56) ,STHE ( 65 ) 
dimension THEniS('J2) ,KF( 32) 

’JO 3 J=1,MPP 
TB.1V( J)=TAPPV( J) 

TBBH( J)=-^APPH( J) 

TBPSH( JJ=C MPLX(TAPPH( J) ,0.0) 

TRTSV( J) =CHPLX( TAPP V(Jl ,0.0' 

STAESH( J)=r,MPLXIO.O,O.0f 
STA6SV( J)=C MPlX( 0.0,0. 0) 

CON T I NUF 

on 12 11 = 1 , A 

CALL HARM' ( 7BFSH,M, I NV ,S ,2 , 'FFRP ) 

CALL HARM ( rpFSV.N, INV,S,2 ,IFEPP) 

or 10 1=1, NOT 

P.FA0( 1 1 .( GTHEHl J) , J = l ,NDP) 

PFAO(l) (GPHIH(J),J=l,NPP) 

RFAD(1> (GTHEV( J) ,J=1,NDP) 

pead(1) 1 gphiv( J) ,j=i ,ndp) 

PCAD(l) (GCPnH( J) ,J=1,NDP) 

READ(i) (CCPOV( J) ,J»] ,NTP) 

0-0 10 J = 1 ,NPP 

J)*(GTHEH( J)-fGPHIH( J) ) )*STHE( I ) 

STAESV( J)=STAF SV( J) + ( TBE.SV( J) * ( GTHFV ( J) -fGPH I V ( J ) ) )*STHE( I ) 

CONTINUF 
REWIND 1 

CALL HARM ( STAFSH,M,!NV,S,-2,IFEPP) 

CALL HARM ( STAE SV,M ,IMV,S ,-P tlFERR) 

DO 11 J=l,NOP 

TAFSH( J)*orAL( STAESHt J) ) 

TAFSV( J )=PFAL( STAcSV(J) ) 

TX1=TAESH( J) 

TX>=TAESV( J) 

TAESH( j )=TX1»PP! ')CFA7X2*CPn 'SF 
TAFSVI J 1=T/2*PF INCF»-TXlwC00S5F 
ERR V( J)=TAPPV( J)-TAF SV( JJ 
ERPHI J )=TArpH( J)- TAF SH( J) 

Tj)BH( J )“TpRH( J ) TERPtU J )»PERHP( J) FFRPV( J) *PERKC( J ) 

TBBVJ J) = TBBV( J)AEPRV( J)*PFPVP( J) *-ERRH ( J ) wPERVC t J ) 

TBFSM(J)=CMPLX(Tbbh( J) ,0.0) 

TBRSV(J)=CHPLX(T00V(J) ,0.0) 

CONTIMUF 

CONTINUr 
^ '•’’UPN 
E'jrt 
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OOUOOOCUUOOO 


IV G LEVFL 


21 


EA1SKY 


DATK = 7A122 


lS/33/Oe 


SllBRnUTINF EftF<’KY (NPMl ,NPH I? , TiF S V, TAESH, THBV , T BflUt L$T atF, FFFF , 

'-1 PHT) 

THIS SUBROUTIMF CALCULATES thf aMTCNNA TEFPEPATURE CONiRIBUTinN DUF Tn THE 
FAPTH AND SKV FUR BOTH PO LAR ' Z ATI ONS. IF LSTATF=6, THE HEABURED ANTENNA 
TLMPEPATUPE-S ARF USFD as THF BRIGH’^NESS TEFFEPATUPFS of TFE EAPTF AND SKY. 
IF LSTATE^"', the results of JUBPOUTTNE tbEAPT ape usfo as THF BPIGHTNFSS 
TFMPEPATUf cS. 

TAESV=THE VFR riCAL ANTE'l-JA TEMPEPATURE CONTP 1 8l)TI CN FPCM THE EARTH 
AND SKY 

■TAESH=THE HOPUUNTAL ANTtNKA TEMPEPATURE CONTbIBUTIQN FROM THE FAPTH 
ANOSKY-^' 

CO«HnN/Bl/TAPPV,TApPH ,STHL- , \DT ,NDP ,DE L . DELTHE , AL P 
1 ,PR INCF.CFnSSF.THE'M S,RF 

COMMON/PLOC I /TBESH, TBFSV.STAESH, stars V, AAAI 
CPMMO'l/BLOC ?/S ,N,1 NV 

CPWMON/06/GTHEH,GPHni,OTHEV,GfHI VtGCRCH.GCROV 
COHPLFX OTHEH(25o) ,GPHIHI2S>il ,GTHEVI2S61 ,GPHIV(256J 
COMPLFX GCPOHt 255) .r.CRCV(2S-.) 

CORPLEX ICE SHI 2551 .CRPLX ,CR , TBE S V ( 255 J , STAE5H ( 256 1 , ST AESVI 255 ) , 
t AAAi: A,?56l - 

DnENRIOM TBBH(25o)’.7B8Vt25’>) ,TBBBH (255) .TBBBVI256 1 
DIMENSION INVI128) ,Sil2P) ,M(3) 

DIMENSION TAESV1 256) ,TALSH( 256) 

DIMENSION TAPPVI256) ,TAPPH(256I ,STHE(65) 

DIMENSICM THFDISI32I ,WF(3?)' 

NALP = ALP/0EL*-5 
NPPZl=NDP/2+1.00l 
■Nl=NOP2I+NPHII 
N2=NDP2I-NPHI2 
CAlP=CnS(ALP| 

P I-3.1A1F9 

CR^CMPLXt 1 . ,0. 1 

IF ILSTA-RF.E0.7) CO TO 1 

DO 2'J = T,NDP 

TBB8IU J ) = TAPPH( J1 

TBBBVt J)=TAPPV(J» 

CONTINUF 
GO ■^'1 A 

continue 

DO B J=1,N0P 
TliBBHI J 1 = TRBHI J) 

■TPBBV! J )=TPBV( J ) 

CONTINUE 
CONTINUE 
00 3 J=1,N0P 
STAESH(J)=CMPLX(0. ,0.) 

STAESVI Jl=CHPLX(0. ,0. ) 

continue 

DO 10- 1=1, NOT 

REAntU' I GTHP.HI J1 ,J = 1 ,NDP) 

KEADCI) ( GPHIHI J) ,J = l ,NDP) 

RLAO(l) (GthEVI J) tJ=l,NDP) 

PFADUl IGPHIV! J) ,J = l-,NDP) 

PEAOni ( GCROHIJl ,J = : ,NDP) 

RFAOU) (CCPOVI J)-,J = l,NPP) 

T(,F=P I-THED] S( 1 1 
00 20 J=1,NDP 

PHI=FL0ATI iAUSI J-N 0P21FNALP) ) ADFL+.OOOl 
TPIIl=TANIPHn 

TH=FH=ATAN(rFFF/(!>H0*CALP»S3PT(l . -ftphI *TPHI ) ) ) 

TllFEW=rHf FMf-PI /2. 

TBrSHI J)=CMPLX( TBBBH(J) ,0.0) 

TBLSVI J ) = CMP LX( TUBB Vi J) ,0. 0) 

IFt U.LE.NDP21 ).AND.( J.CT.N2-3 ) ) GO TO 31 
GO TO 37 

31 TBf SHI J »=CHP.LX( TBBBHf N2-3) , 3, 0) 

TBESyiJ)=CMPLX( T»BBV(N2-3) ,0.0) 

32 CONTINUE 

in(J.GT.NDP21).ANn.l J.LT.NI + 3) ) GO TO 
GO TO 34 • 

33 TPFSH(,Jl=fvpLX(T8BBHINH-3) ,3.0) 

TBFSVI J)=CMPLX( TUB0V(N! + 3) ,0. 0) 

34 CONTINUE 
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Itf G LEVEL ?1 


SAB’SKY 


DATE - 761?2 




7b 


36 

20 


30 

10 


11 


lEt ( J .C.N?! .AND. ( J.LE-.Nl 1 .AND. ( THr.LT.THEEWl ) GO TO 35 
GO Tri 3 6 

TRESHI J l^CMPlXt G., 0. I 
TflE3V(Jl»CMPlX(0.,.3. 1 
lontinuf 

CONTI NOT 

tail HAFK (TBESH,M,INV,5,2,irEPB) 

CALI HARM I TRESV,M,INV,S ,IfFRB) 

on 30 j=i,Nr.p 

STAESHt J1 = STAF SM( J) + (TU£SH( J) *( GTHEH ( J) +GPH I H (J ) 1 )*SThEU I 

J} *fGTHEVUl +GPHIVU 1 ) 1*STHE< n 

COMTiNUE 
CONTFNOE 
PE WIND I 

CALL MAPI* ( STflEStitM,(NV,St'2TlFEPfi) 

TALL HAPP ( STAFSV.N ,’NV, S.-2 , IPEPPl 

on 11 J=1,NPP 

Tt<^SHI J 1=PFPL( STAESHI Jl ) 

TA^SV; J l=PFALt STAESVI JH 
TXi=TAFSH( Jl 
TX2=TACSV( J1 

TAESHI J ) = TX1*PPINCF + TX2*CR3SSP 
TAFSVt J>=YX7*PPINCP<^ rxi*CRe55F 
CONTIMIJC 

prTuRAj 

EMO 
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!V G Ll'VEL 


1 


SPLl NA 


DATE = 7612? 


19/T3/C8 


5()BBnUTlNr- SPLINAl XI , YI ,M 
CSUUPQUTIMF rOF- INTFFPulATIMG .OIPPPRENTI ATING AND 
CINTFGRATIMG A SFT OF DATA POINTS 

CX 1= IMDEPEMOFN-t vapiable df input data, it fust bf in 

CASCENOING ORDER 

cyi=oepfndcnt variable or input oata 
CN=NUNBF 0S or iriPOT data 

DI'^ENSIOM W( 3001,01 3001 ,P tOOOl ,A(300 1 ,C(300 ) ,5(3001 , Z I 300 ) 
DIMENSION X! (N1 ,Y1 ( N) ,X( 3001 , Y(300) 

DATA A(] ),C(1 l,/.{ n/-l. ,0. 0,C.O/ 

OfT 11 I»I,J" 

Xt 1 1 = X1 ( I I 
YUl^YKI) 

11 CONTINUF 
DO 1 J=?,N 

ui ji=xijt-x( j-ii 

1 CONTINIIF 

NN=U-1 
fT ? J=2,NM 
WJ = Wl I ! 

WP=WI J+ll, 

WS*WJ+WP • 

OJ=V'J/WS 
U( J J=OJ 

QA = 1.-.s*oj*a( J-II 
A( J) = .5*( l.-CJI/QA 

B( J) = 3.»(V,P*Y! J-ll-WSTYI Jl (-WJOYIJ+in ZWP/WJ/W 
riJ) = (B(JI-.6*(}JYC( J-11 l/OA 

2 continue 

S(M) = ClNKl/( l.«-AINNIl 
S(NM! = S(N1 
. NHM=N-7 

no 3 JJ=1,NMM 
J=NNM-JJ+1 
S( JI=C(J)-A( 

3 CONTINUE 

rn-A 3=2, N 

Z ( Jl =ZU-1 i+. still Jl +1 Y( JI + Y( J-1 1 -.OB25»W( Jl-t*?* (S ( J l+S I J-1 1) ) 

4 CONTINUr 
RETURN 

ENTRY TrpPAI XV ,YVI 
C INTF.PPOLAT ION-HUNT FOR SUBINTERVAL 
DP 7 JJ = ?,U 
J = JJ 

IF (XV.GT.ytJl) GO T07 
GO TO R 

7 CONTINUE 
ccalculatf function 

8 WJ=WIJ1 

ni = (xv-xij-i)i/wj 

D2 = ( X( J 1-XVI /RJ , 

D?=HJ*WJ/6. 

YV=DI »{ Y( J 1+03*101 to 1-1 .) tSl J) 1 

YV= YV+P2*( Y( J-1 l + 03“f 02 ♦02-1. 1*S 1 J-i I 

RTTUPN 

FNTP Y OEP I VA( XV,YVI 
CDIFFEPFNTIATIPN-HUNT FOP SUBI NT = PVAL 
DO S JJ = 7,N - 
J = JJ 

IF( XV.CT.Xl J) 1 GO TO S 
GO TO 6 

5 CONTINUF 
CCALCULATE ofpivativf 

6 V.J=W(J) 

01-(XV-X(J-111/«J 

D2=( XI Jl-XVl /WJ 

VV=( Y( J 1-Y( J-11 1 /HJ-.S*KJ*((ni«n!-.3333333)*SUl 
!-(0?''02-.3S3333?)*S(J-ll 1 
PFTUPH 

ENTRY GPFATAI XV, YVI 
C integration-hunt FOR SUBINTERVAL 
DO 9 JJ=Z,N 
J = JJ -• 

IF (XV.CT.XJ) 1 GO TO 0 
GO TO 10 

9 rnNTtMiic. 
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IV G LtVU ?l 


SPLINA 


PATfc = 75122 


CCALCULATE FMOFFINITE INTFftr-oAl, 

10 wj=w( j I 

Dl = ( ( XV-X{ J-l) )/HJ)**2 
P2=( I xt J 1- XVI /hJl *■*? 

03^ ,0«25*WJ*WJ 

XV=Z( J-1 1+0. S+WJ*(01+( Y1 J1 + 03*(D1-?, )«S t J) I + U.-D2I+ (Y(J-1 
]+r3*(02-l.J*£( J-Dl 1 
B FTOPM 
cNn 


19 / 33/08 
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LEVEL 21 


HARK 


DATE = 76122 


lo/ 33 /Ofi 


SllBPnOTirF KAf K ( A, M, l NV, S, IF 1 1 FERR I 


HVU15000 


this PRIIGFAF. CALCULATES THE FOtJPIFR tranSFCRK OP THE SAKPLEII FUNCTION A. 


• ■ - ]p 


ni'lENSlQM 

t 

10 IF( lAaS(IFSET) - 1) cco.RO 

12 HTT=HAXO<(M I ) ,M( 21 ,M(3) ) -2 
Pf’OT2=S0PTI 2-OFC) 

IF IMTt-kt I IA.14.13 

13 IFEPE=1 
Rf-TUPN 

14 IfOR-o 

tM = H( 1 > 

H7*Mt 2) 

H3=MI 3) 

N 1-2**131 
tl2 = 2*«M2 
N3=2**M3 

16 IF, I IFS6T) 18,18,20 
la NX= N1*N2*N3 
FN = MX 

Pn lo I = i,NX 
AI 7*1-1 1 = A{2»!-11 /EN 
I” A( 2*1.1 = -A( Z»I l/FN 
20 NP( 

- NCI 21= NP(-l ) *N2 
MP( 3I=NPI 2)*N3 
DO 25P IP=1,3 
IL = Np( 3)-NP(I0) 

IL I = IL + 1 
HI * «( ID) 

IF (HII2SO,25C,3G 
30 IDIF=NP( ID) 

KBIT=NP( ID) 

KEV = 2*<MI/2) 

IF (HI - HRV 160,60,40 

H IS ODD. on 1=1 CAST 
40 K.O 1T=KI<IT/2 
KL=KB!T-2 
DO 50 I=1,!L1,I0IF 
KLA.ST=KL*l 
no 50 K=I,KLA$T,2 
KP=K4KI1IT ■ 

00 ONE STEP WITH L*1,J=0 
AIK > = A(K)+AI KOI 
A(KDI=A(k)-A(Knl 

• TsAIKDI 
A( KO)=A(K l-T 
AIK1=A(K) + T 
T=A(Kn+l> 

AIKD + 1 l*A(K + l)-T 
50 A(K+11=AIK411 +t . 

IF (*I - 1)250,250,52 

S: LFIRST =3 
OFF - JLAST 
JLAST=1 
or TCI 70 


2**1 L-2) -1 


M IS FVFN 
60 LFIPST = 2 
JLA5T=0 

70 on 24C L=LF1PST,HI ,2 
JJOIF=KFIT 
KBlT=KniT/4 
KL*KBIT-2 

DO FDR J=n 

DO no 1=1 ,UI,InIF 

kla,st»i+kl 

on 80 K=I ,KLAST,2 
K 1=K4KBI T 


CHAP 9 30 
CHAR 94C 
CHAR 05C 

CHAP 97C 
DHAO 98C 
OHAB 9SC 
DHAP lOOC 
OHARIOIC 
DHARI02C 
OHAR 1030 
0HAR104G 
DHAP lose 
CHAR 1060 
CHAR 107C 
0HAR1O8C 
DHAP 1080 
CHAR HOC 
OHARlllC 
CHAR112C 
DHARH30 
OHAR 1 14C 
0HARH50 
CHAR 1160 
CHAR 1 17C 
CHAO 1180 
0HARH90 
CHAR 1200 
DHAR1210 
0HAR1220 
DHAR1230 
DHAP 1240 
CHAR 12 50 
CHAR 126C 
OHAR 1270 
CHAR 1280 
CHAR 12SC 
CHAP 1300 
OHAR 13 10 
CHAR 13 20 
OHAR 1330 
OHAR 1340 
DHAR1350 
CHAR 1 360 
OHAR 1370 
OHAR 1380 
CHAR 1390 
OHAR1400 
0HAR1410 
CHAP 1420 
DHAP 1430 
DHAP 1440 
DHAP 1460 
CHAR 1450 
DHAR1470 
OHAR 1480 
CHAR1490 
OHAR 1500 
OHAR 15 1C 
OHAR 1520 
CHAP 1530 
DHAP 1540 
CHAR 15 50 
CHAR 1560 
OHAR 1570 
rHAR15R0 
OHAR I 590 
DHAR16CC 
CHAR 1610 
OHAR 1620 
DHAR1630 
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onrvooooonf'on r-n o o n noonoono'^c-nn 


IV G LFVEL 71 


QATE = 76122 


19/33/Oa 


K2»K1+KB!T 

K’-K2«KBIT 

no TWI STEPS WITH J=0 
A(K-) = AtKH-A(K2) 

AtK2l-A(l' )-A{K2) 

AIK 1 J=AIK1)+A( K3l 
A(«3 I=AIK1)-A(K3) 

A(KJ=A( KI+AI ><1 ) 

A(KI)=A(K )-A(Kl 1 
A{K2>^A(K2)+AIK3)«I 
A(K3I=A(K7)“A( K3l»i 

■ T-=A(K2) 

AIK2I=A(K)-T 
A(K 1 = A(K H-T 
T=A(K2+1) 

A(K2tl)=A(K<-l)-T 
AlK + 1 I=AIK<-1 UT 

T=A(K3) 

A(K31-AIK1)-T 
AtKll=i(Kl)+T 
T=AIK 3 + 1 1 
AlK3+n=AIKl + VI-T 
AIKl+l I=AIK1 + U + T 

T=A(k11 

AI>U1=A(K)-,T 

A(K)+AIK1+T 

T=AIKl+ll 

AIKl+ll=A(K+U-T 

A(K + 1)=A(K + 11 + T 

R=-AtK3 + : 1 
T = A(K31 
A(K3)=A(K2)-P 
A(K2)=A(K2)+P 
AIK3+11 = A(K2 + 1 l-T 
SO AIK2 + 1) =A( K2 + 1 ) + T 

IF ULAST) 235,235, 82 
82 JJ=JJOrf +1 

on FOR J=1 - ’ 

ILAST= IL +JJ 

on 85 I = JJ,IlAST,IOIF 

KLAST = KL + I 

no 85 K=[,KI AST, 2 

K! = K+KBIT 

K7 KI+KI'IT 

K3 = K2+KB1T 

LPTTIMG W=( t + ! )/ftOOT2 ,R3=(-'i+I I /«roT2 ,1.2 = 1 , 
A(IO = A(K)+A(K2)*I 
A(K?)=A(Kl-A(K2)*i ■ 

A (Kn = A(Kl)*V>+ AIK 31*1.3 
A(K3)=A<K1 I*U-A(K3J*W3 

A(K»=A(KI+AtKl) 

A(K1)=A(K1-AIK1) 

A(K2)*A(K21+A(K.31*I 

A(K3)=AlK7)-A(K3)*i 

R =-A|K2+n 
T = A(K?) 

AIK2J = Al K)-P 
AIKl = AIKl+P 
A{K2+1J=A(K + 1)-T 
A( K + 1 )=A( K + n + T 

AWn = AIK!>-A{Kl + H 
AWI = AIKl + l I+A(K1I 
r=-A(k31-A|K3+1) 

T=A{k3)-AI K3 + 1) 

A(K3)»(AKP-R I /POOT2 


DHAR16A0 
CHAR1650 
DHAR1660 
0HAB1670 
CHAR 1680 
CHAP 1650 
CHAR 1700 
CHAR 17 10 
CHAP 1720 
DHAR1730 
0HAR17A0 
DHAR1750 
DHAR1760 
DHAR177C 
CHAP 1780 
0HAR1790 
CHARI SCO 
OHAR lalO 
CHAR 18 20 
DHAR183C 
CHARISAO 
DHAR1850 
OHAR I860 
DHAR187C 
CHAP 1880 
DHAR1850 
CHAR 1900 
OHAR 19 10 
OHAR 1920 
OHAR 1930 
0HAR1940 
CHAR 1950 
CHAR 1960 
DHAR1970 
DHAR 1980 
CHAR 1990 
0HAR2000 
DHAR2Q10 
DHAR 20 20 
CHAR 20 30 
OH AP 2040 
CHAR 2050 
CHAR2060 
CHAP 2 070 
CHAR 2080 
OHAR 2090 
OHAR2100 
OHAR 2 1 10 
OHAR 2 120 
DHAP2130 
0HAR2140 
OHAR 2 1 50 
OHAR 2 160 
CHAP217C 
DHAR2180 
DHAR2190 
CHAR 2200 
CHAR2210 
DHAR2220 
CHAR2230 
0HAR2240 
0HAR2250 
DHAR2260 
OHAR 2270 
-CHAR22R0 
DHAR229C 
0HAR23O0 
DHAR2310 
CHAR2320 
DHAR 2 3 30 
CHAR2340 
DKAR235C 
DHAR 23 60 
0HAR2370 
DHAR2380 
DHAR 2 390 
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rt(K3»I ) = (/lWI-T) /R00T2 


DHAR 2400 


A(K I) = ( AWP+R )/F00T2 


DHAR241F 


A(K1*-1)=(AWH-T) /R00T2 


0HAR242C 


T= Al K1 1 


0HAR2430 


A( K1 ) = A(K )-T 


0HAR244C 


A(K.)=A1KH-T 


DHAR245C 


TsA(Kl^l) 


CHAR 2460 


A(Kl4-l) = A(K+] l-T 


CHAR247C 


A{R»-ll=A(K*l ) + T 


DHAP248C 


P=-A(K3+1 ) 


CHAR2490 


T = A(K3) 


DHAR 2500 


A(K31-AIK2I-R 


CHAR251C 


AtK2)=A(JC2)+R 


DHAR 2520 


a(k3H1 = A(K2 + 1I-t 


DHAR253C 

8S 

A(K2An=A(,K2 + l) + T 


0HAR254C 


IF(JtAST-l) 235,235,00 


DHAR255C 

<)0 

JJ* JJ + JJOIF 


0HAR2S60 

C 



CHAR2570 

c 

NOW on THE RSHAININO J’S 


OHAR256C 


on 230 J-?,JLAST 


DHAR 2 59 C 

c 



DHAR 2600 

c 

FETCH M's” 


DH AR 2 610 

c 

Opc- W=K*»INV( J) , W?=W**-9, W3=V 


CHAP262C 

06 

I*INV(JH) 


DHAR 26 10 

9fl 

U=MT-I 


DHAR 2640 


w(ii=sun 


DHAR 2650 


K( 7>=S( I 1 


DHAR 2660 


1?=2*I- 


DHAR267C 


12r=t|T-12 


DHAR 2630 


IF( 120120,110,100 


CHAR2690 

C 



DHAR 2 700 

C 

2*1 IS IN FIRST OUADRANT 


DHAR271C 

100 

HZ( U = SII2CI 


CHAR 2 7 20 


W2I21 = S[ 12) 


OHAR273C 


00 Til 130 


DHAR 2740 

uc 

W2( 1) = 0. 


0HAR275C 


W?I2)=1-. 


DHAR2760 


CO TCI 130 


0H4R2770 

c 



DHAR 27EC 

c 

2» I IS If; SECOND QUA3FAN'' 


CHAR2790 

120 

I2CC = I2C+NT 


DHARZnCG 


I2S=- I?C 


GHAR201C 


WZm=-SI17C) 


DHAR 2820 


H2(2)=S( I2X.C 1 


DHAR2B30 

130 

1 7*1 + 12 


DHAR 2840 


I3C.=MT-I3 


CHAR 2 a 50 


I FI I 3C I16C, 15C,140 - 


0HAR2860 

c 



OKAR 287C 

c 

13 IN f IF ST QUADRANT 


DHAR 2880 

140 

WT( 11=S( I3C ) 


DKAR2890 


W?I 21=Sl 13) 


DHAR29C0 


GO TO 200 , 


DHAR 2910 

• 150 

W3I 1)*0. 


DHAR 2920 

• 

W312)=1. 


DHAR 2930 


(,n *n 500 


DHAR 2940 

c 



CHAR2950 

160 

I3i;C=I3C +NT 


DHAR 2960 


|F( I3CC )1°0,1PC,170 


CHAR 29 70 

c 



0HAR29B0 

c • 

I'l IN SFCOMO OUAORAKT 


DHAR2990 

no 

I 3C =- I 3C 


DHAR 3000 


W3I 1)=-S(I3C) 


CHAR3010 


K3( 2)=SI I3CC ) 


CHAR 3020 


on TO 200 


CHAP, 30 10 

180 

W3I11*-1. 


0HAP3040 


H3(2)=0. 


CHAR 3050 


C TO 200 


0HAR3C6C 

c 



CHAR3070 

c 

3*1 IN THIRD OUADRAN'' 


DHAR3Q80 

190 

I 3CCr.=NTtl 3CC 


DHAR 3090 


l?Cf = -I3CC 


CHAR 3100 


H3tll=-S( I'3CCC) 


0HAR3110 


W3(2)=-S( T3CD 


DHAR312C 

700 

ILAST*IL+JJ 


CHAR 3 130 


DO 220 I=JJ,ILAST,IDIF 


DHAR314C 


KLAST*KL + I 


0HAR315C 
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UP 220 K = I *KLAST,2 


DHAR316C 



KI’sK^KBlT 


DHAR31TC 



K2=Kl+KniT 


DHAR3180 



K3=K2+K0I T 


0HAR31SC 

C 




DHAR 3200 

C 


00 TWO STFP.S with J WET 0 


DHAR3210 

c 


A'('<) = A(K)+A{1<21*W2 


DHAR 3 22 C 

C 


a(K2I=a(K)-A[K2)*H2 


DHAP3230 

c 


A(Kl) = A(Kn-WfAtK3I*W3 


0HAP32AC 

c 


A(K31=A(K1 |»k-AI K3I*W3 


DHAR 3 2 sc 

c 




DHAR 3 260 

c 


a(k)=;A(k)+aiki ) 


DHAR327C 

c 


A(K11»AIKI-A1K1) 


DHAR3280 

c 


A(K2J = MK2)4A(K3)*I 


DHAR3290 

c 


A(W3I=A(!C2)-MK3)*I 


DHAR3300 

c 


* 


DHAR331C 



K=AIK21*V2( l)-A<K2<-ll»W2(31 


0H4R332C 



T=A(K2)»-W2< riwAI K2+ 11 *H2 III 


0HAR3330 



A(K2)=AIK1-R' 


DHAR334C 



A(K) = AIKlwP 


DHAR335C 



A(K2+l)-AfK*ll-T 


DHAR336C 



A(KH 1=AI K + l l+T 


DHAR337C 

c 




OMAR 328C 



Pa.A(K ?)*W3(4 l-A(K3+ 1 J *W3<2I 


CHAR33'>0 



T = A(K3)*W3(2)+A(K3-m»W3(ll 


OKAR3400 



AWW^AIKl )*WI1)-A(K1 + 11 »H(21 


CHAR34U 



AW!=AIKl)“W(2l+A(Kl+ll’»Wm 


DHAR3A2G 



AIK3I=AWP-P 


0H4R343C 



A(K3+1I=AVI-T 


CHAR344C 



A(K1|=AHP+R 


DHAR 3450 



A(K1»U=AWI*-T 


DHAR3460 



T=A»k1) 


DHAR 3A7C 



A(K1)=A(K)-T 


CHAR 3480 



AIK >=!A(k)wT 


DHAR3410 



T = A|KH-1 J 


DHAR 3500 



A(K1+1)=A(KH l-T 


CHAR3S10 



AIK + I >=A( K + 1 l + T 


0HAR352O 



R=-A(K3+1) 


GHAR3530 



1=A(K3I 


DHAR 3540 



A(K3l=AIKrl-P 


DHAR 3550 



A(K21=A(K2)«P 


CHAP 3560 



A( K iH l=AI K2M l-T 


DHAR3570 


220 

AIK2W1 ) = A(K2+1 H-T 


DHAR35F0 

c 


ENO or I AND K LOOPS 


CHAR 3 590 

c 




DHAR 3600 


230 

JJ=JJOIF*J'J 


OHAR3610 

c 


‘fcNO PF J-LHOP 


DHAR 3620 

c 




DHAR 3630 


23S> 

JLAST = 4'=JLASTw3 


OHAR3640 


2^0 

CO^^T^^UJF 


OHAR36SO 

c 


EMO OF L LOOP 


CHAR 3660 

c 




0HAR367C 


250 

CPN71NUF 


DHAR3680 

c 


EMO OF ID LOOP 


CHAR36<5C 

c 




CHAR 3 700 

c 


HE 'low HAVE ’•HE COMPLEX rOU’IFF SU''S 

BUT THFIR AOORCSSES 

ARE 0HAR3710 

c 


lilT-RFVEPSEO. THE FOLLOWING ROUTINE 

PUTS THFM Ih ORDER 

0HAR372C 



NTSC=NT<-Nr 


CHAR 37 30 





DHAR 3740 


350 

IFIkl^PTI 37C,360,36C 


DHAR375C 

c 




CHAR 3760 

c 


M3 GP. OR EO. MT 


0HAR37T0 


360 

100 3=1 


DHAR378C 



N3VN^=N3/MT 


DHAR 3790 





DHAR 3800 



c,n Tp 3180 


0HAR3810 

c 




CHAP 3820 

c 


*13 LESS THAM MT 


CHAR3830 


370 

ion ^=2 


0HAR384C 



N 3VNT^l 


char 38 so 



nTVN‘»^NT/N3 


DHAR 3850 



.'1IMN3=N3- 


0HAR387C 


3B0 

JJ03 ^ NTSO/N? 


CHAR 3880 



M2'^TsM2-MT 


0HAR3890 


450 

IF (M2MTI47C, 400,450 


OH AR 3900 

c 




CHAP 39 10 
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GATE = T',122. 


iR/is/oe 


c 


c 

c 


c 

c 


c 

c 


c 

c 

c 

c 

c 


H2 GP. OF EO. >*1 
•460 1 GO 2=1 

N2VNTs=N2/NT 

KINN7=NT 
GO TO 480 


M? L-rSS THAN MT 
470 IG02 = 2 
N2VNT = ] 

NTVN?=-NT/M2 
HINN2=N2 
460 JJD?=NTbO/N? 
H1MT=N]-HT 

550 IFIM1MT157Q, 560,560 

Ml Gf- . OP FO. MT 
560 I GO 1=1 

N1VNT-=N1/NT 
MINNI =NT 
GO TO 580 


Ml I.FSS THAN MT 
570 I GO 1=2 
N1VNT=1 
NTVN1=NT/NI 
HINN1=N1 
5B0 JJOI=NTSO/M1 
600 JJ2=1 
J=l 

00 880 JPP’=1 ,NTVNT 
IPPT=INV( JJ3) 
on 870 JP?=1,MINN3 
GO TO I 610,620), ICOT 
610 IP3=INV<JP3)»N3VNT - 
■ GO TO 630 
620 IP3=ri'iVlJP3)/NTVN3 
630 I 3=( 1?P3+!P3)*H? 

700 JJ2=1 

DC 870 JPPP=1,N2VNT • 
IPP2=!NV( JJ21 + I3 
DO 86C JP2=1,I'!NN2 
GO TO (710,720,1002 
710 1P2=INV< JPPl^NPVNT 
CO TO 7^0 

720 IP2=1NV(JP?1 /NTVN2 
730 I2=( IPP2 + IP2)*N1 
800 JJl=l 

DO 860 4PP1=1,N1VNT 
IPP 1=INV( JJl )♦! 2 - 
on 850 JP1=1 ,KIMH1 
GO Tn t eio, 820) , IGOl 
810 IP1 = INV(JP11*MIVNT 
GO TO 8 30 

820 IP!=rNV(JPl)/NTVNl 
830 ! = ?*( IPPl-HPD + l 

IF (4-1 ) PAC.PEC.eSC 
840 T=A(( ) 

A< 1 )=A( 4) 

A( 4 1=T 
T=A( I H ) 

A( I+1)=A(4<-1) 

A(4H) = T 
850 4 = 4 >-r 
860,441-=441 + 4401 

P70 442=442+4402 

bKD OF 4PP? AND 4P3 LOOP.' 


880 443 = 443+4403 
END OF 4PP3 LOOP 

890' 1F( IFSET)R91,P.95,89‘' 

891 DO 392 I-= 1 ,NX 

892 A( 2*1 ) = -A(2»I ) 

895 prTU“N 


OHAR30 2O 

0HAR3930 

0HAR3940 

CHAR 39 50 

DHAR3960 

0HAR3970 

OHAP 3980 

0HAR3990 

DHAP4000 

CHAR 4010 

CHAP 4020 

DHAR4030 

CHAP 4040 

0HAR4050- 

DHAR406Q 

GHAR4070 

0HAR4069 

nHAR4090 

GHAR4100 

CHAR 4 1 10 

0HAR4120 

0HAP4130 

CHAR4140 

CHAR 4 150 

0HAR414C 

DHAR4170 

0HAR418O 

0HAR4150 

0HAR4200 

0HAR4210 

CHAR 4 2 20 

CHAR 42 30 

0HAR4240 

OHAR42 50 

CHAO4260 

DHAR4270 

OHAP 4 2 80 

CHAR 4290 

OHAR4300 

0HAR431O 

DHAR4320 

CHAR4330 

□HAR4340 

CHAP 4350 

0HAR4360 

DHAR4370 

OHAP 4380 

DHAR4390 

0HAR440C 

CHAP 4410 

DHAR4420 

DHAP.4430 

CHAR 4440 

DHAR445Q 

DHAR4460 

CHAR4470 

OHAR4480 

OHAP 4690 

CHAR4500 

OHAP4510 

0HAR4520 

CHAR 4 5 30 

CHAR 4 5 50 

DHAR4560 

OHAP 4580 

0HAR459O 

DHAR4600 

OHAP 4 6 10 

CHAR4620 

DHAR4630 

OHAR464 0 

DHAR4650 

DHAR4660 

0HAR467C 

CHAR4680 

DHAP4690 


214 



(•so no oon oo on r>oo no 


IV C LEVEL 21 


HARM 


DATE = 76122 


TH= rOLLOWlNG PROGRAM ClIMPUTES THE SIN AND I NV TABLES. 

900 MT=MAXOfH( 1 ) ,Kf 2»,M( 31) -2 
•AT = HAX0I?,MT) 

«0A IE (HT-J81 9C6.906,13 
906 IPER2=0 
M T= 2*#MT 
NTV2=NT/2 


SET UP SIN TABl 
THETA=PIE/2**(L«-i) FCR L = 
910 thctA=.78539P1 6339794 

= 2**IMT-L»-l) FDR L= 
=NT 


jstfp 

JSTEP 


jotr=2*«(HT-L) Fn° L=i 
J0!F-N’'V2 ■ 

S1J0IF)=SIN1THFTA) 

00 9R0 L=2,PT 

theta* thDTA/2. 000 

j$TcP2=JSTEr 
J STEP=jriF 
J01F = J5TFP/2 
S(.ID1FI = SIN(THE’^41 
JC1=NT-JDIF 
S( JCl l=C0$I THETA) 
JLAST=NT-JSTEP2 
!F(JL4ST - JSTEP) 959,920,920 
920 00 940 J*JSTEr,JLAST, JSTEP 
JC=NT-J 
J0=J+jniF 

940 SIJD1 = SIJ)*S( JC1)«-S( J0!F)«SI JO 
950 continue 

S'T lip IMV(J) TABLE 
960 MTt E)<P=NTV2 


fnLEXP=2A'»lMT-Ll . TOR L=1 
LMIEXP*] 


lmifxp=2*«(l-i). fop L*l 

INVt 1 1»0 
on 990 L=1,NT 
IflV(lM)FXP + n = HTLFXP 
on 970 J = ?,LM1FXP 
JJ=J+L«1EXP 

970 INV( .IJ )=IMV( J)4MTLF XP 

mtlExp^mtlexp/7 

980 LM1EXP*L«1EXP»2 
982 IF( IFSFT112, £95,12 
END 


19/33/08 


DHAR470O 
0HAH471C 
DHAR4720 
0HAP4730 
OHARA 7 40 
DHAR4770 
0HAR4780 
DHAR4790 
CHAR4800 
DHAR4810 
DHAR4820 
CHAR 48 30 
0HAR4B4C 
DHAR4850 
DHAR4860 
0HAR4870 
□HAR4880 
CHAP 4890 

0HAP4910 
0HAR4920 
CHAR 49 30 
DHAR4940 
CHAP 4 9 SC 

DKAP4970 

CHAR4990 
OHAR5000 
DHAR 5010 
CHAR5020 
OHAR5030 
0HAR5O40 
DHAR50S0 
DHA95060 
OHAP5070 
OHAR508O 
DHAR5090 
OMAR 5 100 
DHAR 5 110 
OHAR5120 
DHAR 5130 
CHAP 5 140 
DHAR5150 
DHAR 5160 
OHAR5170 
OHAR5180 
DHAR 5 190 
0HARS20O 
0HAR521O 
0HAR5220 
0HAR5230 
DHAR5240 
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I G LEVEL 21 


DATE = 76122 


ig/SJ/oa 


SURAOIJTJHC ICSSGU ( X, V,Dy ,S,N ,A ,B,C,D,HKI 
•ICSSGU s library 1 


FlINCTinN 

USAGE 

PARAMFTCPS X 


■RPFCISION 

LANGUAGE 

LATEST REVISION 


- CUBIC SPLINE DATA SNCOTHING 

- CALL ICSS.''U(X,Y,CY,S,N,A,B,C,0,WK) 

- VECTDP OF N ABSCISSA xm.xizl, 

• ..,XIN1. WC ASSUME XII 1 .GT.Xtl-n 
FOR I =2 ,3,. , A 

- VECTOR OF N FUNCTIONAL VAUJFS Ytll,Y(21, 

.^.,Y(N). Y(I) IS THE FUNCTIONAL VALUE 

- VECTOR OF INPUT PARAMETFPS DY ( 1 ) , GY I 2 I , 

,OYtN) where OVm IS AN ESTIMATE OF 

THR STANDARD DEVIATION OF THE QRCINATF Y( I I 

- A NON“N=GATI VE PARAMETER WHICH CONTROLS THF 

EXTENT OF SMOOTHING 

- NUMBER OF DATA POINTS 

' OUTPUT VECTOR OF SALINE COEFFICIENTS 
A( 1) ,A( 2) .. .. .AIN+l) 

- OUTPUT vector of SPLINE COEFFICIENTS 

0(11 ,0121 ,B1 N+11 

- OUTPUT VEC^DP OF SPLINE COEFFICIENTS 

Ctl 1 ,C<2) ,C(N+1) 

- OUTPUT VECTOR OF SPLINE COEFFICIENTS 

0( 1 1 ,0(21 ,D(M+1) 

THF CUBIC SPLINE FtXX) IS CONSTRUCTED IN 
TH= I-TH ( I «l ,2 . ,N-11 INTERVAL BY 

F( XXI =( t0( IFI) YH+C (1+11 J*H+B(lFl) )*H+At 1+11 
WHERE H = XX-xm 

- WORK AREA OF OIHENSICN .GE. 7*N+9 

- SINGLE 

- FORTRAN 


- SBP.TFMRRR ?.q, 1972 


SUBROUTIMF ICSSGU ( X,.Y . DY, S ,N ,A ,B , C ,C . HR I 


ni MEN SION 


X(ll,Yni,OY{ll,A(ll,H(ll,C(ll,0(ll,WK(l) 
SET UP WORKING AREAS 


IBl = NPl ■ 

IP2 = IBl+NPl 
IB? = IB2+NPI+1 
lOA = IB3+NP1 
IBS = IBA+NPl 
IBS = IB‘^+NPl + 1 
WKI 11 = O. 

WK ( 2 1 = 0. 

WKIIB2I = 0. 

WKUB3) = 0. 

IJK2 ^ IB2+NP1 
WKI IJK21 = 0. 

IJKS = IB5 + 1 
WK( IJKS)=D. 

UKS = IBS + 2 
WKI IJK5)-0. 

WKI IB61 = -0. 

IJKS = IB5+NP1 
WKI IJKS) = 0. 

•H=X<*2 l-XI 1) 

F2»-S 

F=( Y( 21-YI 1) 1/H 
IF (M.LT.-3I GO T1 25 
DO 5 l=?,N 

G=H . 

H=X( II- XI I- 11 

E = F ■ 

F=I Y( I l-YI I-in/H 

A( I l=F-r 

IJK1 = IB3+1 

WKI IJKA 1=1 G+H)-». 6666567 

IJKA = lRA+1 

WKI IJKAI=H/3. 

IJK2 = IB2+I 


ICMU0020 
-ICKU003C 
ICMUOOAO 
ICHU0030 
ICMU0060 
ICMU0070 
ICHUOOeC 
ICMUOOSC 
ICHUOIOD 
ICMUOllC 
ICMU012C 
ICMU013C 
ICMU016C 
ICHUOISC 
ICHU016Q 
ICMU017C 
ICMUOIBC 
ICHU0190 
ICMU020C 
ICMU021C 
IC MUD 2 20 
ICMU023C 
1CHU02AC 
ICHU025D 
ICHU026C 
ICMU027C 
ICNU02B0 
ICMU029G 
1CMU030C 
ICMU0310 
IC MU0320 
ICHU0330 
-ICMUO3A0 
ICNUO3S0 
ICMU036C 

ICHU038C 

ICMU039C 

1CMUD6D0 

ICMU041C 

ICMU0420 

ICHU0430 

ICMUOAAD 

ICHUOASC 

ICMU046C 

ICMU0A7C 

ICMU04SC 

ICHUOASC 

ICHU050Q 

IC.MU051C 

ICHU052C 

ICMU0530 

ICHU05A0 

ICMU055C 

ICHU0560 

1CMU057Q 

ICMU058C 

ICMU0590 

ICMU0600 

ICMU061C 

ICMU0620 

ICMU0630 

1CM0064C 

1CHU0650 

ICMU0650 

ICMU067C 

ICMU0630 

1CMU0690 

ICMUQ70C 

ICHU0710 

ICMU0720 

ICHU0730 

1CMU0740 

ICMU0750 

ICMU076Q 
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DATP = 76122 


l9/33/0<i 


WKl IJK21=nV( 1-2) /G 
WK( [1 = 0V( I )/H 
TJKl = lEl+I 

WK( IJKl l = -rY( 1-1 ) /G-DY(T -1 )/K 
5 COMTINUF 
nn 7 1^=2, N 

iJK2*in2+i 

H(!> = WKl n*WKU )+WKl lJKl)*WKn JK1)+WK! IJK?)»WK( !JK2) 

CM) = WKII )»WK{ t ,JK1 + : )f/(K (IJKl)*WK( IJK2 + 1) 

RI 1) = WKI I )*MK( I JK2 + r,) 

7 CONTINUE 

NEXT ITEPATION 

10 ir (N.lT,3) go to 25 
DO IE 1 = 3, N 
IJKl = 

IJKD = 1-1 
NKdJKl )=F7NK(1JK0) 

1JK2 = 1B2+I-2 

IJKO = 1-2 

WK( IJK?J=G=WKn JKO) 

IJKO = I 
IJK3 = 1R3+I 

WK( 1JK01 = 1. /(P*0(I )fWK(I JK3)-P*WK( r JKl )-G*WKUJK2) ) 

UK5 = IBEfl 
UKN = 1JK5-1 
IJKO = IJKN-1 

WKIUKE) = Ain-HKl IJK1)*HK(1JKK1-WK(UK2 )*WKUJKCi) 

IJK-A = IPA+I 

F=P«C( I )»-V,Kt IJKA)-H*WK(I JKl) 

G=i| 

H=0( I )*P 
15 CPtITINUr 

no 20. 1 = 3, N 

J=N- 1+3 
1JK5 = IB5+J 
IJK6 = IJK5+1 
1JK7 = IJK6+1 
IJKl = 161 + J 
1JK2 = IR7+J 

NK(1JK5) = WK( J) #WK(1 JK51-WKU JKl )*WK( 1JK6 )-WKl IJK7)*WK( IJK71 
20 rONTiNUF 
25 E = C 
H=0 

' CCMFUTE U AND ACCUMULATE E 

no 30 I=2,K 

G=H 

IJK5 = IB5+I 

H = ( KKl I JK5+1)-WK(I JK5) 1 / (X( n-Xtl-1 ) ) 

IJK6 = 1R6+1 

WKt IJKf ) = (n-G)*OY(l-l ) *I)Y( I-I ) 

E=E+WK( I JK6)*(H-G) 

30 CONTINUF 

G=-H*DY(N)*0Y(N) 

IJK6-S IP6+NP> 

HK( 1JK6)=0 

G=F2 

j. 2=j:+n*p 

.IFIF2.GF.S .OP. F2.LE.G1 GO TO AS 

r = o. • 

IJK6 = 1B6+2 

H = [WKnjK6+3)-wKn jK6n/(x(2)-x(in 
IF(N.LT.3) GO TO 40 
no 35 1 = 3, N 

G=ll 

1JK6 = TB6+! 

« = ( WKt I JK6+1)-WK(! JK61) /tX( I )-XI I-l) 

IJKl = lBl+I-1 
IJK2 = IB2+1-2 

■ r, = U-G-WKt ! JKl )+WK(I-l)-WKn JK2 I *WK( 1-2) 

F = F+G*WK( I )»0 
WK m = G 
35 CONTimiF 
40 >I=E-P»F 

IF(H.LE.O) GO TO 45 

LPDATf ruF LAGPANGE PULTIPLIEH p 


ICMU077C 
ICMUC78C 
ICMU079C 
ICMU080G 
ICHUOaiC 
1CMU082C 
1CMU083C 
ICMU084C 
1CMU08 5C 
ICMU086C 
1CHU087C 
ICHU088C 
1CHU089C 
IC“UO90C 
ICMUOOlO 
ICKU0920 
1CHU093C 
ICMU0940 
1CMU0950 
1CMU096C 
ICHU0970 
ICMU098C 
ICHU059C 
ICMUIOOO 
ICMUlOlO 
1CMU1020 
1CMU1030 
1CMU1040 
ICMU1050 
ICMU1060 

!CM'J1070 

ICMU1080 

1CMU1090 

ICMUllOO 

ICMUlllO 

1CMU1120 

ICHUU30 

ICMUU40 

1CMU1150 

ICMU1160 

ICMU1170 

ICMU118C 

ICMU1190 

ICKU1200 

ICMU1210 

ICMU1220 

ICHU1230 

ICMU124G 

1CMUL250 

1CMU126C 

1CMU127C 

ICM01280 

1CMU129C 

ICMU13CC 

ICMU1310 

ICMU1320 

1CMU133C 

1CMU1340 

1CMU1350 

ICMU136C 

ICHU1370 

ICMU1380 

ICMU1390 

1CKU1400 

ICMU1410 

ICMU1420 

ICMU1430 

1CMUL440 

1CMU1450 

FCHU1460 

ICMU1470 

1CMU148C 

ICMU1490 

1CM01500 

ICMUI510 

ICMU1520 
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on 


rcssGu 


1 9/33/08 


C I EVfL 21 


55 

9005 


P = 
C(1 


+ CS-F2)/((SQRT( S/EI+D+H) 


RATE = 76 122 
FOR THE NEXT ITE°ATI0N 


10 


45 on 50 !=2,NP1 

IJK6 - 

At I )=Yn-n-P»WK( IJK5 I 
1JK5 = IR5+I 
Cl 1 ) = V.K( ; JK5) 

50 rONTINUE 

on 55 1=2, N 

H=X( I )-Xt i-n 
01 I l=tr t I*U-C( I ) ) /( 3.=HI 
CnNI NIJF n/H-(H*OtI)fC(in'‘H 


IF E LESS THAN OR OCIIAL TO S, 

COMPUTE THE COEFCICIFNTS AND RETURN. 


R F ruR N 
END 


ICKU1530 

ICHU1540 

1CKU1550 

ICHU156C 

ICMU1570 

1CHU1580 

ICMU159C 

ICHU1600 

ICNU1610 

ICNU162C 

ICHU1630 

ICNU164C 

ICMU16EC 

ic^uioeo 

KMU1670 

ICHU168C 

ICM01690 
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